Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,487,225 T→C D126D (GAT→GAC)  unknown → Flp pilus assembly membrane protein TadE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,487,2250TC97.9% 129.2 / NA 48D126D (GAT→GAC) unknownFlp pilus assembly membrane protein TadE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/39);  minor base A (1/0);  total (9/39)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCTATCACTGAAAGTGCAAGAATTGCTAAAAATAGAACTGCAGAGGGTAATGATTATGCGGTAGAATTTGAGAAAGCTCTAAAACAACAGCTTGTTTATCAAGAATCTTCAACTTTAGGTTATCTCGCCCGTTTGGATAAAAATGGCGGTTACAAGATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGTTCTTGTGTTAAGTCTTTATTGGATGGTAGGTTTAGACGGCCTACTAAAGATAGAGACGGTAAGATTATACCCCCC‑AATGGTCGGTTAGCAACATTGGCTGTATATTCATTGACCTACAATTATGAGTTTCTGGTTTCACTGCCATTTTTACCTAAGGATAGTGTTGGCGGTTTATTGTCACGTAAGTTTGTAGCCGTACAGGAATTTAATAGATCTAAATTCCAATATCCCGGCAGGTCAAATAAATAATCAGGCAATTT  >  1/1486980‑1487434
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                  
ttagccataaCTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACagaaa                                                                                                                                                                                                               <  1:602452/242‑3 (MQ=255)
        aaCTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTTTTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTaata                                                                                                                                                                                                       >  1:560789/2‑249 (MQ=255)
         aCTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTaatat                                                                                                                                                                                                      <  1:8590/251‑4 (MQ=255)
              aaGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATAc                                                                                                                                                                                                 <  1:332321/251‑1 (MQ=255)
              aaGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATAc                                                                                                                                                                                                 <  1:303701/251‑1 (MQ=255)
              aaGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATAc                                                                                                                                                                                                 <  1:482367/251‑1 (MQ=255)
              aaGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATAc                                                                                                                                                                                                 <  1:606470/251‑1 (MQ=255)
                gTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTTATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTTTTAGATGGTAAGTTTAGGCAGCCTACAAAAGAAAGAAACGGTAATATTATAc                                                                                                                                                                                                 >  1:485284/1‑249 (MQ=255)
                gTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTaatat                                                                                                                                                                                                      <  1:120730/244‑4 (MQ=255)
                        atcgCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATgg                                                                                                                                                                                      <  1:237468/248‑1 (MQ=255)
                         tcgCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGt                                                                                                                                                                                     <  1:223428/249‑1 (MQ=255)
                             tAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGt                                                                                                                                                                                 <  1:151568/251‑1 (MQ=255)
                                      aacggcggaaggTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAAc                                                                                                                                                                          <  1:323401/239‑1 (MQ=255)
                                             gaaggTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTg                                                                                                                                                                  <  1:640058/247‑1 (MQ=255)
                                             gaaggTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTg                                                                                                                                                                  <  1:110985/247‑1 (MQ=255)
                                                     cgATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCAt                                                                                                                                                         <  1:480908/250‑1 (MQ=255)
                                                     cgATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCAt                                                                                                                                                         <  1:399960/250‑1 (MQ=255)
                                                       aTTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCTCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCAtt                                                                                                                                                        <  1:634449/250‑1 (MQ=255)
                                                       aTTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCAtt                                                                                                                                                        <  1:784424/250‑1 (MQ=255)
                                                         tatgctgcaGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCAttaac                                                                                                                                                     <  1:841116/243‑4 (MQ=255)
                                                                  gAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTaatat                                                                                                                                                                                                      <  1:168314/194‑4 (MQ=255)
                                                                  gAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCa                                                                                                                                                          <  1:247496/237‑1 (MQ=255)
                                                                    aTTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCt                                                                                                                                                   <  1:724858/242‑1 (MQ=255)
                                                                                 ctAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATgg                                                                                                                                                                                      <  1:507792/194‑1 (MQ=255)
                                                                                            acttatttATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCCAAATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCAttt                                                                                                                   <  1:160468/246‑1 (MQ=255)
                                                                                                     tCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATTAAATTGATGTTTAATATGTTGATTGTGGGAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGGAGGTTAGCAACATTGGCTGTGTATTTATTAACCTACAAGTATTAATTTTTTGTTTTGCTGCCATTTTTACCTaggg                                                                                                         >  1:492464/1‑248 (MQ=255)
                                                                                                      cAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCa                                                                                                                                                          <  1:871881/201‑1 (MQ=255)
                                                                                                            acctCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCt                                                                                                                                                   <  1:854501/199‑1 (MQ=255)
                                                                                                             cctCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAAATGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGGTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGGATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATaacgtt                                                                                                 >  1:162447/3‑246 (MQ=255)
                                                                                                                  actatAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AAt                                                                                                                                                                                        <  1:856712/155‑1 (MQ=255)
                                                                                                                  actatAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGtagg                                                                                                                                                                                  <  1:450764/161‑4 (MQ=255)
                                                                                                                  actatAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTAc                                                                                                               <  1:163409/228‑1 (MQ=255)
                                                                                                                     atAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATaacgttagcag                                                                                            <  1:483397/247‑11 (MQ=255)
                                                                                                                            aTCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATTTTTATTAGATGGGAAGGTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATAACCCCC‑AATGGTAGGTTAGCAAAATTGGCTGGGTATTCATTAACCTACAAGTATGAAATTTTGGTTTCGCTGCCATTTTTACCTAGGGATaacgtttgcaggttgttgtta                                                                                  >  1:649348/1‑231 (MQ=255)
                                                                                                                                               aaaTGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGt                                                                           <  1:392231/239‑1 (MQ=255)
                                                                                                                                                     cggttatgaaaTTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGaa                                                          <  1:164474/240‑1 (MQ=255)
                                                                                                                                                        ttatgaaaTTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAAc                                                                                                                                                                          <  1:872591/128‑1 (MQ=255)
                                                                                                                                                          atgaaaTTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCTTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAAt                                                    <  1:424603/246‑1 (MQ=255)
                                                                                                                                                          atgaaaTTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATaacg                                                                                                   >  1:875529/6‑201 (MQ=255)
                                                                                                                                                                  gATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATta                                                                                                                                                                                                    <  1:425287/100‑1 (MQ=255)
                                                                                                                                                                     gTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATaacgttagc                                                                                              <  1:446746/198‑9 (MQ=255)
                                                                                                                                                                       tAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATaacg                                                                                                   >  1:839709/1‑188 (MQ=255)
                                                                                                                                                                                    tgtgGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAATAGATCTAAAtttcagcatccaggtag                          <  1:462618/251‑16 (MQ=255)
                                                                                                                                                                                     gtgGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAATAGATCTAAAtttcagcatccaggtagg                         >  1:355380/1‑235 (MQ=255)
                                                                                                                                                                                     gtgGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAATAGATCTAAAtttcagcatccaggtagg                         >  1:38067/1‑235 (MQ=255)
                                                                                                                                                                                                                                gTTTAGGCAGCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAATAGATCTAAATTTCAGCATCCAGGTAGGTCAAACAaataattaattcatga  <  1:768933/231‑11 (MQ=255)
                                                                                                                                                                                                                                         gCCTACAAAAGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAATAGATCTAAAt                                          <  1:129261/182‑1 (MQ=255)
                                                                                                                                                                                                                                             acaaaaGACAGAAACGGTAATATTATACCCCCC‑AATGGTAGGTTAGCAACATTGGCTGTGTATTCATTAACCTACAAGTATGAATTTTTGGTTTCGCTGCCATTTTTACCTAGGGATAACGTTAGCAGTTTGTTGTCACGTAAGTTTGTAGCCGTACAAGAATTTAATAGATCTAAATTTCAGCATCCAGGTAGGTCaaaca                  <  1:502462/199‑3 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                  
TTGGCTATCACTGAAAGTGCAAGAATTGCTAAAAATAGAACTGCAGAGGGTAATGATTATGCGGTAGAATTTGAGAAAGCTCTAAAACAACAGCTTGTTTATCAAGAATCTTCAACTTTAGGTTATCTCGCCCGTTTGGATAAAAATGGCGGTTACAAGATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGTTCTTGTGTTAAGTCTTTATTGGATGGTAGGTTTAGACGGCCTACTAAAGATAGAGACGGTAAGATTATACCCCCC‑AATGGTCGGTTAGCAACATTGGCTGTATATTCATTGACCTACAATTATGAGTTTCTGGTTTCACTGCCATTTTTACCTAAGGATAGTGTTGGCGGTTTATTGTCACGTAAGTTTGTAGCCGTACAGGAATTTAATAGATCTAAATTCCAATATCCCGGCAGGTCAAATAAATAATCAGGCAATTT  >  1/1486980‑1487434

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑