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

GGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACC  >  NC_000913/752987‑753262
                                                                                                                                          |                                                                                                                                         
ggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                          >  1:1325577/1‑140 (MQ=255)
 gTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                         >  1:115450/1‑140 (MQ=255)
   tataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaa                                                                                                                                       >  2:1105808/1‑140 (MQ=255)
    ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGGTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                      >  2:190429/1‑140 (MQ=255)
    ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                      >  2:419179/1‑140 (MQ=255)
    ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTAAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                      >  1:150743/1‑140 (MQ=255)
          gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                >  1:939723/1‑140 (MQ=255)
          gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                >  2:162649/1‑140 (MQ=255)
          gCGACTACAATAAAGTTTTTAATTGAAAGTATTGGTTTGCTGATAATTTGAGTTTTTTTTTTTTTTTTAAATATCTATATAGGTTTGTAAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAAcc                                                                                                                                <  2:847917/140‑1 (MQ=255)
                         tttttAATTAAAAGTATTGGGTTGCTGATAATTTGAGTTTTTCTTTTTTTTTTAAATATTTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAAcggc                                                                                                                 <  1:1424842/140‑1 (MQ=255)
                                                tcTGATAATTGGAGTTTTTCTTTTTTTTTTAAATTTCTATATAGGTTTTTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTAGGGCATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                          <  2:1235426/139‑1 (MQ=255)
                                                 cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                         >  1:193580/1‑140 (MQ=255)
                                                 cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                         >  1:1461165/1‑140 (MQ=255)
                                                  ttaTAATTTGAGCTGTTCTTTTTTTTTAAAATATCTATATAGGTCTGTTAATGTATTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCa                                                                                        <  2:1461165/138‑1 (MQ=255)
                                                      aaTTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                                    <  1:190429/140‑1 (MQ=255)
                                                       ttttGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                   <  1:1105808/139‑1 (MQ=255)
                                                       aTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                   >  2:1248619/1‑140 (MQ=255)
                                                           gAGCTTTTTTTTTTTTTTAAAATTTTTATATGGGTCTGTAAAGGGATTTTTTTTTTAAAATTTTTTTTTTTTAGTAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                               <  2:230193/140‑1 (MQ=255)
                                                           gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                               >  2:861252/1‑140 (MQ=255)
                                                            aGCTTTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGATTTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                              <  1:1210190/140‑1 (MQ=255)
                                                                gTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                          >  1:272295/1‑140 (MQ=255)
                                                                      ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                                    >  2:488050/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACCACTTAGCAATCCACCATTAAAGCTTTAAATCgt                                                               >  2:820815/1‑139 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTTAAGCTTGATATCGc                                                               >  1:1397075/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  2:374639/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  2:460752/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  1:204044/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACCGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                               >  2:132297/1‑140 (MQ=255)
                                                                           tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTAAAATTAGCAAACAACCATTAACCCTTGATATCGc                                                               >  2:1318551/1‑140 (MQ=255)
                                                                                   tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                                       >  1:281545/1‑140 (MQ=255)
                                                                                        ataGGTTTGTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                  <  2:939723/140‑1 (MQ=255)
                                                                                         taGGTCTTTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                 <  1:162649/140‑1 (MQ=255)
                                                                                         taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                 <  1:1248619/140‑1 (MQ=255)
                                                                                         taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                 <  1:419179/140‑1 (MQ=255)
                                                                                         taGGTCTGTTAATGGATTTTATTTTTAAAATTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                 <  1:1341875/140‑1 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCTACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  2:555083/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  2:247987/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  1:904060/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  1:1313096/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  1:186244/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  1:197464/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  1:402619/1‑140 (MQ=255)
                                                                                            gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAAAGCTTGATATCGCTTTTAAAGTCGCGtttt                                              >  1:744743/1‑140 (MQ=255)
                                                                                              tttttAATGTATTTTTTTTTTACATTTTTTTTTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                            <  2:1192570/137‑1 (MQ=255)
                                                                                              ttttaaaTGGATTTTATTTTAAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                            <  1:275153/135‑1 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  2:436250/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  1:1134567/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  2:565369/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  1:439110/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                           >  2:1125885/1‑140 (MQ=255)
                                                                                               tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAACAACCATTAAAGCTTGATATCGCTTTTTAAGTCGCGTTTTTct                                           >  2:275153/1‑139 (MQ=255)
                                                                                                tttAATTGTTTTTTTTTTTACAATTTTTTGTTTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  2:1129451/139‑1 (MQ=255)
                                                                                                tttAATGGTTTTTTTTTTTACATTTTTTTTTTTTAAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  2:204044/139‑1 (MQ=255)
                                                                                                tttAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  2:1325577/139‑1 (MQ=255)
                                                                                                tttAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                          <  1:76348/139‑1 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                         >  2:1379649/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                         >  2:719213/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                         >  1:922567/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                         >  2:113985/1‑140 (MQ=255)
                                                                                                 ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                         >  1:1192570/1‑140 (MQ=255)
                                                                                                   aaTGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                       <  1:362034/140‑1 (MQ=255)
                                                                                                       gATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTTCAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGGTTTTCATATTCTGt                                   >  1:383064/1‑140 (MQ=255)
                                                                                                        tttttattttTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                  <  1:1427435/139‑1 (MQ=255)
                                                                                                         ttttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCCTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGtat                                 >  1:1129451/1‑140 (MQ=255)
                                                                                                           ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtataa                               >  1:913880/1‑139 (MQ=255)
                                                                                                           ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                               >  1:342628/1‑140 (MQ=255)
                                                                                                                   ataaTTTTTTGTTTTTACGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                       <  2:400249/138‑1 (MQ=255)
                                                                                                                      tttttttGTGTTGAGCAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                    <  2:383064/139‑1 (MQ=255)
                                                                                                                        ttttttttttTAGACATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                  <  1:374536/140‑1 (MQ=255)
                                                                                                                        ttttttTGTTAAGGAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                  <  1:820815/140‑1 (MQ=255)
                                                                                                                            tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACAGCTGACGCGGACGGGc              >  2:40437/1‑140 (MQ=255)
                                                                                                                            tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc              <  1:565369/140‑1 (MQ=255)
                                                                                                                            tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc              <  2:922567/140‑1 (MQ=255)
                                                                                                                                tttAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATc          <  1:374639/140‑1 (MQ=255)
                                                                                                                                  aagcaataaAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt        <  2:934674/135‑1 (MQ=255)
                                                                                                                                        tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc  >  2:519748/1‑140 (MQ=255)
                                                                                                                                        tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc  >  2:1165853/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACC  >  NC_000913/752987‑753262

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: