Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1050AT21.2% 205.5 / 23.6 80intergenic (‑310/+80)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (49/14);  new base T (0/17);  total (49/31)
Fisher's exact test for biased strand distribution p-value = 2.61e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGA  >  NC_000913/752971‑753239
                                                                                                                                      |                                                                                                                                      
gaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAttttt                                                                                                                                   >  1:94248/1‑140 (MQ=255)
 aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  2:1575035/1‑140 (MQ=255)
 aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  2:366225/1‑140 (MQ=255)
 aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  2:699742/1‑140 (MQ=255)
 aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  1:1638706/1‑140 (MQ=255)
 aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                                    >  1:598300/1‑138 (MQ=255)
 aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                                    <  2:598300/138‑1 (MQ=255)
  cgccATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGGATTGGGTTGCTTATAATTTGAGCTGGTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                 >  2:522732/4‑140 (MQ=255)
  cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                 >  2:1185973/1‑140 (MQ=255)
  cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                 >  1:1962915/1‑140 (MQ=255)
  cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                 >  1:1441656/1‑140 (MQ=255)
      aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  <  2:1728997/135‑1 (MQ=255)
      aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  1:1728997/1‑135 (MQ=255)
      aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGtt                                                                                                                             <  1:1575035/140‑1 (MQ=255)
          tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                                         <  2:1393148/140‑1 (MQ=255)
             tAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCat                                                                                                                      >  2:585163/1‑140 (MQ=255)
                ggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                   <  2:463636/140‑1 (MQ=255)
                  ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                 >  1:1230638/1‑140 (MQ=255)
                    ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                               >  2:1733062/1‑140 (MQ=255)
                    ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                               >  1:1340279/1‑140 (MQ=255)
                          gCGACTACCATGAATTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGTTTTGTTAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAaaccaacc                                                                                                         <  2:2043394/140‑1 (MQ=255)
                          gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                         >  1:1393587/1‑140 (MQ=255)
                           cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                        >  2:854428/1‑140 (MQ=255)
                            gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                       >  2:738380/1‑140 (MQ=255)
                            gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGCTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                       >  1:1911453/1‑140 (MQ=255)
                               tACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCa                                                                                                    >  1:994628/1‑140 (MQ=255)
                                 ccATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCata                                                                                                  >  2:1808739/1‑140 (MQ=255)
                                 ccATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCata                                                                                                  >  2:1399061/1‑140 (MQ=255)
                                    tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGa                                                                                               >  2:1892147/1‑140 (MQ=255)
                                    tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGa                                                                                               >  1:1836323/1‑140 (MQ=255)
                                        ttttttAATTGAAAGTTTTGGGTTGCTGATAATTTGAGTTTTTTTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAAcgg                                                                                           <  2:1783477/139‑1 (MQ=255)
                                         tttttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggc                                                                                          <  1:1198272/140‑1 (MQ=255)
                                         tttttAATTAAAATTTTTGGTTTGCTGATATTTTGAGTTTTTTTTTTTTTTTTAAATATCTATATAGGTCTGTAAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAAcggc                                                                                          <  2:1393587/140‑1 (MQ=255)
                                           tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                                        >  2:178952/1‑140 (MQ=255)
                                           tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAT‑AAAATCAACCCGCCATATGAACGGCggg                                                                                       >  2:95442/1‑140 (MQ=255)
                                                ttGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAT‑AAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                  <  1:95442/140‑1 (MQ=255)
                                                 tGAAAGTTTGGGGTTGTTGATATTTTGAGCTTTTTTTTTTTTTTAAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                  <  2:1962915/140‑1 (MQ=255)
                                                   aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                >  1:978002/1‑140 (MQ=255)
                                                        tttGGGTTGCTGATAATTTGAGTTTTTCTTTTTTTTTTAAATATCTATATAGGTTTTTTAATGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                           <  2:1836323/139‑1 (MQ=255)
                                                                   aaTAATTTGAGTTTTTTTTTTTTTTTTAAATATCTATATAGTTTTGTTAATGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGCTTAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                <  2:1638706/139‑1 (MQ=255)
                                                                      attttGAGTTGTTCTTTTTTTTTTAAATATTTATATAGTTTTTTTAATGTTTTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                             <  2:994628/138‑1 (MQ=255)
                                                                                       tttttttAAATATCTATATAGGTCTGTTAATGTTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGata                                            <  1:585163/138‑1 (MQ=255)
                                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:395004/1‑140 (MQ=255)
                                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:186332/1‑140 (MQ=255)
                                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:1621803/1‑140 (MQ=255)
                                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:1279496/1‑140 (MQ=255)
                                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:111768/1‑140 (MQ=255)
                                                                                                   tCTATATAGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTTTTTAGGCATATAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                <  1:1733062/140‑1 (MQ=255)
                                                                                                   tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                >  2:1094672/1‑140 (MQ=255)
                                                                                                         taGGTCTGTTAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGAATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                          <  1:854428/140‑1 (MQ=255)
                                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  2:1809533/1‑140 (MQ=255)
                                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  2:871769/1‑140 (MQ=255)
                                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  1:1108078/1‑140 (MQ=255)
                                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  2:1549510/1‑140 (MQ=255)
                                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  1:848137/1‑140 (MQ=255)
                                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  1:84002/1‑140 (MQ=255)
                                                                                                              ttGTAAATGGTTTTTTTTTTTAAATTTTTTTTTTTTTAGGAATAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                     <  2:972955/139‑1 (MQ=255)
                                                                                                              cTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                     <  1:1094672/140‑1 (MQ=255)
                                                                                                              cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCCTGATATCGCTTTTAAAGTCGCGTTTTTc                     >  1:588810/1‑140 (MQ=255)
                                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCa                    >  1:616629/1‑140 (MQ=255)
                                                                                                                tttAGTGGTTTTTATTTTTATATTTTTTTTTTTTTAGGAATAAAAAGCTCAGTCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:1673612/139‑1 (MQ=255)
                                                                                                                tttAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:738380/139‑1 (MQ=255)
                                                                                                                tttAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:1279496/139‑1 (MQ=255)
                                                                                                                tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:1167718/139‑1 (MQ=255)
                                                                                                                tttAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:84002/139‑1 (MQ=255)
                                                                                                                tttAATGGATTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:1673601/139‑1 (MQ=255)
                                                                                                                tttAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:1738762/139‑1 (MQ=255)
                                                                                                                ttaaaTGTATTTTTTTTTTACATTTTTTTGTTTTTGGGCAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:1108078/137‑1 (MQ=255)
                                                                                                                gTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:1340607/140‑1 (MQ=255)
                                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTAAACTTTAGAATACAACATTAACGCTTGATATCGCTTTTCAAGACGCGTTTTTAata                  >  2:764027/1‑140 (MQ=255)
                                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                           <  1:750537/131‑1 (MQ=255)
                                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                           >  2:750537/1‑131 (MQ=255)
                                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                  >  1:1167718/1‑140 (MQ=255)
                                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                  >  2:872375/1‑140 (MQ=255)
                                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                  >  1:1180740/1‑140 (MQ=255)
                                                                                                                   aaTGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                <  2:588810/140‑1 (MQ=255)
                                                                                                                   aaTGGATTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                <  1:1239157/140‑1 (MQ=255)
                                                                                                                   aaTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                <  1:1809533/140‑1 (MQ=255)
                                                                                                                      ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg             >  1:972955/1‑140 (MQ=255)
                                                                                                                        tttttattttTACATTTTTTTTTGTTTAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATTTCGCTTTTAAAGTCGCGTTTTTCATATCCTGta           <  2:134047/139‑1 (MQ=255)
                                                                                                                          ttttttttTAAATTTTTTTGTTTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata         <  2:848137/140‑1 (MQ=255)
                                                                                                                           ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAATTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc        >  2:529047/1‑140 (MQ=255)
                                                                                                                                 ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATTAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTCTATACAGCTGa  >  2:571153/1‑140 (MQ=255)
                                                                                                                                 ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAcctgtct  >  1:166389/1‑134 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGA  >  NC_000913/752971‑753239

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