Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA21.4% 181.8 / 16.7 70intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (47/8);  new base A (0/15);  total (47/23)
Fisher's exact test for biased strand distribution p-value = 6.80e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTT  >  NC_000913/752986‑753256
                                                                                                                                           |                                                                                                                                   
gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                                                     <  1:1238510/140‑1 (MQ=255)
      taAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAAt                                                                                                                               >  2:680904/1‑140 (MQ=255)
         aTGCGACTCCCATGAATTTTTTAATTGAAATTTTTGGTTTGTTAATATTTGGAGTTTTTTTATTCTTTTTAAATATCTATATAGTTCTTTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCaa                                                                                                                            <  2:109615/140‑1 (MQ=255)
          tGCGACTCCCATGAATTTTTTAATTAAAATTTTTGGTTTGCTAATATTTTGAGTTTTTCTTTTCTTTTAAAATATCTATATAGGTCTTTTAATGAATTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAAc                                                                                                                           <  2:1222116/140‑1 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGCGTTTAGGCATATAAAAATCAAcc                                                                                                                          >  1:140397/1‑140 (MQ=255)
             gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                        >  2:919938/1‑140 (MQ=255)
              aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                       >  1:109937/1‑140 (MQ=255)
              aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                       >  1:212285/1‑140 (MQ=255)
                tACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCa                                                                                                                     >  1:142939/1‑140 (MQ=255)
                  ccATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCata                                                                                                                   >  2:1336130/1‑140 (MQ=255)
                          tttttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTTTTTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAACCAACCCGCCATATGAAcggc                                                                                                           <  1:1280705/140‑1 (MQ=255)
                           ttttAATTAAAATTTTTGGGTTGCTGATATTTTGAGTTTTTTTTTTTTTTTAAAATATTTATATGGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAAcggcg                                                                                                          <  2:692066/140‑1 (MQ=255)
                                    aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTaaaaa                                                                                                 >  2:564565/1‑139 (MQ=255)
                                      aGTATTGGGTTGCTGGTAATTTGAGCTGTTCTCTTCTTTTTAAATATCTATATAGGTCGGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGGGGGGTAAAatat                                                                                               >  2:202999/1‑140 (MQ=255)
                                      aGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAatat                                                                                               >  2:360940/1‑140 (MQ=255)
                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                      >  1:1336645/1‑140 (MQ=255)
                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                   >  2:912486/1‑140 (MQ=255)
                                                     aTAATTTGAGCTTTTTTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                <  2:543976/140‑1 (MQ=255)
                                                             aGCTTTTCTATTTTTTTTAATTATCTATATAGTTCTGTTAATGGTTTTTATTTTTACATTTTTTTTTGTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                        <  2:1134655/140‑1 (MQ=255)
                                                                    cTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                 >  1:1134655/1‑140 (MQ=255)
                                                                       ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                              >  2:470916/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTGTTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                         >  2:258096/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTATCGCTTGATATCGc                                                         >  1:345492/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                         >  2:231079/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                         >  1:776602/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                         >  2:203465/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                         >  2:1366138/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGc                                                         >  2:896848/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAACAATCAACCATTAACGCTTGATATCGc                                                         >  1:1300002/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTTAAATATTTACAACATAGCAATCAACCATTAAAGCTTGttatacc                                                         >  2:927432/1‑137 (MQ=255)
                                                                                          tgggTTTGTTAATGGATTTTATTTTAACATTTTTTTTTTTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                           <  2:1297455/138‑1 (MQ=255)
                                                                                          taGTTCTGTAAAGGATTTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                           <  2:1201170/140‑1 (MQ=255)
                                                                                          taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                           <  1:964182/140‑1 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                        >  2:982185/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                        >  1:75514/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                        >  2:269474/1‑140 (MQ=255)
                                                                                              tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCCCTTTTTAAGTCGCGgtttt                                       >  1:1316028/1‑140 (MQ=255)
                                                                                              tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTTAACTTAGCAATCAACCCTTAACGCCTGATATCGCTTTTAAAGTCGCGttttt                                       >  1:1432092/1‑140 (MQ=255)
                                                                                               ttGTTAATGGATTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                      <  2:714131/139‑1 (MQ=255)
                                                                                               cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAAATTAGCAATCAACCCTTAACGCTTGATATCGCCTTTAAAGCCGCGtttttt                                      >  1:1172814/1‑139 (MQ=255)
                                                                                               cTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                      >  1:450055/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTct                                     >  1:710084/1‑139 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                     >  2:558698/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                     >  2:1159068/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                     >  2:238751/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                     >  1:1170402/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAAGACGCGTTTTTCa                                     >  2:1441699/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                     >  1:1201170/1‑140 (MQ=255)
                                                                                                 tttAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                    <  1:844100/139‑1 (MQ=255)
                                                                                                 tttAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                    <  2:450055/139‑1 (MQ=255)
                                                                                                 tttAATGGATTTTTTTTTTAAAATTTTTTGTGTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                    <  1:680904/139‑1 (MQ=255)
                                                                                                 ttaaaTGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                    <  2:212285/137‑1 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCATCCCGCCATATGAACGGCGGGTTACAGTATTTATAACTTCCCAATTCACCATTAAAGCTTTGTATCGCCTTTAAAGTCGGGTGTTTCata                                   >  2:171803/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                   >  1:1096248/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                   >  1:1049455/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                   >  2:1189019/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                   >  1:1151437/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                   >  1:829290/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                   >  1:1421333/1‑140 (MQ=255)
                                                                                                   tAATGGATTTTTTTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat                                  <  1:258096/140‑1 (MQ=255)
                                                                                                    aacccatttttttttAACATTTTATGGGTTTTCGGCATATAAAAATCAGCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTCACGCTCGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                 <  2:1432092/135‑1 (MQ=255)
                                                                                                    aaTGTATTTTTTTTTTAAATTTTTTTGTGTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                 <  2:1096248/140‑1 (MQ=255)
                                                                                                     atgttttttattttTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                <  1:360940/135‑1 (MQ=255)
                                                                                                       ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                              >  2:246443/1‑140 (MQ=255)
                                                                                                          ttttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCCTTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                           >  1:236719/1‑140 (MQ=255)
                                                                                                            ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                         >  1:1178727/1‑140 (MQ=255)
                                                                                                                        ttttttGTTTTAAGAAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa             <  2:1085881/140‑1 (MQ=255)
                                                                                                                         tttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc            <  1:1366138/140‑1 (MQ=255)
                                                                                                                           tttGTTTTTAGACATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACgg          <  1:269172/140‑1 (MQ=255)
                                                                                                                           tttGTGTTAAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACgg          <  1:564565/140‑1 (MQ=255)
                                                                                                                             tttttttAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCTTATCCTGTATACAGCTGACGCGGACGGGc        <  2:710084/136‑1 (MQ=255)
                                                                                                                             tgttTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc        <  2:971874/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc        <  1:246443/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc        <  1:470916/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc        <  1:912486/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc        <  1:1159068/140‑1 (MQ=255)
                                                                                                                             tgcttttAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc        <  2:140397/136‑1 (MQ=255)
                                                                                                                              gtgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa       <  1:238751/140‑1 (MQ=255)
                                                                                                                               tgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCaa      >  1:971874/1‑140 (MQ=255)
                                                                                                                                   tAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt  <  1:231079/140‑1 (MQ=255)
                                                                                                                                   aagaaaaataAAAAAAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt  <  2:1174066/133‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTT  >  NC_000913/752986‑753256

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