Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1050AT31.2% 143.6 / 31.5 80intergenic (‑310/+80)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (45/10);  new base T (1/24);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 1.67e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACG  >  NC_000913/752970‑753241
                                                                                                                                       |                                                                                                                                        
agaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                                      <  1:73253/140‑1 (MQ=255)
agaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                                      >  1:1560604/1‑140 (MQ=255)
 gaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAttttt                                                                                                                                     >  1:1409612/1‑140 (MQ=255)
  aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                    >  1:358520/1‑140 (MQ=255)
  aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                    >  1:343653/1‑140 (MQ=255)
   cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATGTAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                   >  2:573131/1‑140 (MQ=255)
   cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                   >  1:631959/1‑140 (MQ=255)
   cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGAATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACCATTTTTtg                                                                                                                                   >  2:1410179/1‑140 (MQ=255)
    gTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtgt                                                                                                                                  >  2:668172/1‑140 (MQ=255)
     tCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtgtg                                                                                                                                 >  1:258627/1‑140 (MQ=255)
                gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                                      <  2:168977/140‑1 (MQ=255)
                   ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtatata                                                                                                                   >  2:160760/1‑138 (MQ=255)
                   ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                   >  2:314709/1‑140 (MQ=255)
                   ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                   >  2:360599/1‑140 (MQ=255)
                           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                           >  1:469740/1‑140 (MQ=255)
                           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                           >  1:7536/1‑140 (MQ=255)
                           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                           >  1:454179/1‑140 (MQ=255)
                           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTAGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                           >  1:1351668/1‑140 (MQ=255)
                            cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                          >  2:119588/1‑140 (MQ=255)
                            cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                          >  2:1490843/1‑140 (MQ=255)
                             gACTACCATGAAGTTTTTAATTGAAAGTTTTGGGTTGCTGATAATTTGAGCTGTTCTATTCCTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                         >  1:1313963/1‑140 (MQ=255)
                             gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                         >  1:285913/1‑140 (MQ=255)
                             gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                         >  1:767216/1‑140 (MQ=255)
                             gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                         >  1:1103680/1‑140 (MQ=255)
                                          tttttAATTAAAAGTTTTGGGTTGCTAATATTTTGAGCTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAAcggc                                                                                            <  2:1531412/140‑1 (MQ=255)
                                                  taaaaGTTTTGGGTAGTTGATATTTTGAGTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                    <  2:1351668/138‑1 (MQ=255)
                                                  taaaaGTATTGGTTTGTTAATAATTTGAGTTTTTTTTTTTTTTTTAAATTTCTATATAGGTCTGTTAATGTATTTTTTTTTTACATTTTTTTGTGTTTACGATTAAAAAAATCATCCCGCCATATGAACGGCGGGTTaaa                                                                                    <  2:1560604/138‑1 (MQ=255)
                                                  taaaaGTATTGGTTTGCTGATAATTTGAGCTTTTTTTTTCTTTTTAAATATTTATATAGGTCTGTTAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                    <  2:469740/138‑1 (MQ=255)
                                                  tGAAAGTATTGGGTTGCTGATAATTTGAGCTTTTTTTTTTTTTTAAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                    <  2:1046100/140‑1 (MQ=255)
                                                  tGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                    >  1:785390/1‑140 (MQ=255)
                                                   aaaaGTTTTGGTTTGCTGATAATTTGAGCTTTTTTATTTTTTTTAAATATTTATATAGGTCTTTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaaa                                                                                   <  2:851264/139‑1 (MQ=255)
                                                       tttttGGGTTGTTAATATTTGGGGTTTTTTTTTTTTTTTAAAATTTTTATATAGGTTTTTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATAtt                                                                               <  2:785390/137‑1 (MQ=255)
                                                       tttttGGGTTGCTAATATTTTGAGCTTTTTTTTTTTTTTAAAATATCTATATAGTTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATAtt                                                                               <  2:84978/137‑1 (MQ=255)
                                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                       >  2:433594/1‑140 (MQ=255)
                                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                    <  1:160760/140‑1 (MQ=255)
                                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                    >  1:1153303/1‑140 (MQ=255)
                                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTATTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTGAAAAATTTACAACTTAga                                                                    >  2:1487896/1‑139 (MQ=255)
                                                                     aTAATTTGAGTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTGTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                 <  2:1360145/140‑1 (MQ=255)
                                                                     aTAATTTGAGCTTTTTTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGTTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                 <  2:454179/140‑1 (MQ=255)
                                                                     aTAATTTGAGCTTTTTTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                 <  2:1191510/140‑1 (MQ=255)
                                                                       aaTTTGAGTTGTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                               <  2:285913/140‑1 (MQ=255)
                                                                       aaTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                               >  1:7230/1‑140 (MQ=255)
                                                                        aTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                              >  1:1191510/1‑140 (MQ=255)
                                                                            ggtcttttCTTTTTTTTTTAAATTTTTATAAAGTTATTTAAATGGATTTTTTTTTTAAATTTTTTTGTGTTGAGTAATAAAAAAAAAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                          <  2:7536/134‑1 (MQ=255)
                                                                             ggtttttCTTTTTTTTTTAAATATCTAAATAGGTTTGTTAATGATTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                         <  2:358520/135‑1 (MQ=255)
                                                                              gCTTTTTTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                        <  2:254170/140‑1 (MQ=255)
                                                                                 gTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                     >  1:1003813/1‑140 (MQ=255)
                                                                                       ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                               >  2:416326/1‑140 (MQ=255)
                                                                                         ttttttAAATATTTATATAGTTTTTTTAATGTATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGTGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                             <  2:426552/139‑1 (MQ=255)
                                                                                          tttttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                            <  2:1103680/140‑1 (MQ=255)
                                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                          >  1:921191/1‑140 (MQ=255)
                                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                          >  1:1134265/1‑140 (MQ=255)
                                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                          >  1:570257/1‑140 (MQ=255)
                                                                                            tttAAATATCTATATAGGTCTGATAATGGATTTTATTTTTACAATTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGCTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGatatatct                                         >  1:1074962/1‑136 (MQ=255)
                                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                         >  2:581308/1‑140 (MQ=255)
                                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                         >  1:158370/1‑140 (MQ=255)
                                                                                                              tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                        >  1:608315/1‑140 (MQ=255)
                                                                                                               tttttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                       <  2:218395/137‑1 (MQ=255)
                                                                                                               ttGTTATTGGATTTTTTTTTTACATTTTTTTGTGTTTAGCCAAAAAAAATGCAGCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                       <  1:1487896/139‑1 (MQ=255)
                                                                                                               ttGTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                       <  2:1003813/139‑1 (MQ=255)
                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                      >  1:1221337/1‑140 (MQ=255)
                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                      >  2:1342964/1‑140 (MQ=255)
                                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                      >  1:1033288/1‑140 (MQ=255)
                                                                                                                 tttAATGTATTTTATTTTTACAATTTTTTGTGTTTGGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  1:360599/139‑1 (MQ=255)
                                                                                                                 tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  2:1033288/139‑1 (MQ=255)
                                                                                                                 tttAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  2:1289729/139‑1 (MQ=255)
                                                                                                                 tttAATGGTTTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  1:1053488/139‑1 (MQ=255)
                                                                                                                 tttAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  2:951212/139‑1 (MQ=255)
                                                                                                                 ttaaaTGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  1:573131/137‑1 (MQ=255)
                                                                                                                 tatAATGAATTTTATTTTTAAAATTTTTTGTTTTAAGAAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  2:1382453/138‑1 (MQ=255)
                                                                                                                 gTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  2:608315/140‑1 (MQ=255)
                                                                                                                 gTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                     <  1:314709/140‑1 (MQ=255)
                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                    >  2:489891/1‑140 (MQ=255)
                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                    >  2:660896/1‑140 (MQ=255)
                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                    >  1:252894/1‑140 (MQ=255)
                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                    >  1:185639/1‑140 (MQ=255)
                                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTCata                    >  2:1045859/1‑140 (MQ=255)
                                                                                                                    aatgttttttttttttAAATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                  <  2:1029545/134‑1 (MQ=255)
                                                                                                                    aagggtttttttttttAAATTTTTTTGTTTTTAGGATTAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                  <  1:118645/134‑1 (MQ=255)
                                                                                                                    aaTGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                  <  2:158370/140‑1 (MQ=255)
                                                                                                                    aaTGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                  <  2:1391412/140‑1 (MQ=255)
                                                                                                                    aaTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                  <  1:660896/140‑1 (MQ=255)
                                                                                                                         tttttttttttACATTTTTTTGTTTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta             <  1:1093722/139‑1 (MQ=255)
                                                                                                                            ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc          >  1:1472247/1‑140 (MQ=255)
                                                                                                                                    aCAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg  >  1:1156864/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                        
AGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACG  >  NC_000913/752970‑753241

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