Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA27.7% 134.3 / 23.7 65intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (40/7);  new base A (0/18);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 9.65e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752986‑753264
                                                                                                                                           |                                                                                                                                           
gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                                                             >  1:186975/1‑140 (MQ=255)
   ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                                          >  2:296909/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  1:1448129/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                  >  2:2267/1‑140 (MQ=255)
            cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                                 >  1:1620579/1‑140 (MQ=255)
                 aCCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCat                                                                                                                            >  2:1089722/1‑140 (MQ=255)
                        attttttAATTGAAATTTTTGGGTTGCTGATAATTTGATTTGTTCTTTTTTTTTTAAATATTTATATAGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGAAAATAAAAATCAACCCGCCATATGAAcg                                                                                                                     <  2:694637/138‑1 (MQ=255)
                          tttttAATTAAAAGTATTGGTTTGCTGATAATTTGAGTTTTTTTTTTTTTTTAAAATATCTATATAGGTCTTTTAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAAcggc                                                                                                                   <  2:1620579/140‑1 (MQ=255)
                            tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                                                                 >  1:100254/1‑140 (MQ=255)
                                    aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                                         >  1:638421/1‑140 (MQ=255)
                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                              <  2:248208/140‑1 (MQ=255)
                                                 gCTGATAATTTGAGCTGTTCTATTTTTTTAAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                            <  1:1453791/140‑1 (MQ=255)
                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGc                                                                                           >  2:1514187/1‑140 (MQ=255)
                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTTCAACTTAGc                                                                                           >  2:1615138/1‑140 (MQ=255)
                                                               cTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                              >  1:32371/1‑140 (MQ=255)
                                                                  ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGc                                                                           >  2:1317228/1‑140 (MQ=255)
                                                                  ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGc                                                                           >  1:345331/1‑140 (MQ=255)
                                                                       ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                                      >  1:1381443/1‑140 (MQ=255)
                                                                        ttttttaAAATATCTATATAGTTCTTTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGata                                                                     <  2:129994/138‑1 (MQ=255)
                                                                         ttttttAAATATTTATATGGTTTTGTTAATGGATTTTTTTTTTACAATTTTTTTTGTTTAGGAAAATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                    <  2:1421840/139‑1 (MQ=255)
                                                                         ttttttAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                    <  2:371875/139‑1 (MQ=255)
                                                                         ctctttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                    >  1:217638/4‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTGGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  2:360047/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCCTTAACGCTTGATATTGc                                                                 >  1:777363/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:920094/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:1626613/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:803166/1‑140 (MQ=255)
                                                                                      tttataGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTc                                                       <  2:345331/138‑1 (MQ=255)
                                                                                         ataGGTGTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                    <  1:1264723/140‑1 (MQ=255)
                                                                                          taGTTCTTTTAATGTTTTTTTTTTTTACAATTTTTTGTTTTTAGGAATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   <  1:1068209/140‑1 (MQ=255)
                                                                                          taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTGGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   <  1:360047/140‑1 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                >  2:578375/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                >  2:631789/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                >  2:687614/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                >  2:1520453/1‑140 (MQ=255)
                                                                                               tttttAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                              <  1:1153577/137‑1 (MQ=255)
                                                                                               ttGTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                              <  2:638421/139‑1 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGgtttttt                                             >  2:396340/1‑138 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  2:453202/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATAACGCTTTTAAAGTCGCGTTTTTCa                                             >  2:1297675/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttttta                                             >  1:651649/1‑138 (MQ=255)
                                                                                                 tttATTGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:819202/139‑1 (MQ=255)
                                                                                                 tttAATGTATTTTATTTTTACAATTTTTTGTGTTTAGACATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:186975/139‑1 (MQ=255)
                                                                                                 tttAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:508386/139‑1 (MQ=255)
                                                                                                 tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:990281/139‑1 (MQ=255)
                                                                                                 tttAATGGATTTTTTTTTTAAAATTTTTTGTGTTTAGGCAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:32371/139‑1 (MQ=255)
                                                                                                 tttAATGATTTTTTTTTTTACATTTTTTTTTTTTAAGAAAAAAAAAAATAACACCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:1106289/139‑1 (MQ=255)
                                                                                                 tttAAGGGTTTTTATTTTTAAAATTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:688632/139‑1 (MQ=255)
                                                                                                 gTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:1297675/140‑1 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTtata                                           >  2:27587/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:191573/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:935479/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:840306/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:508386/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:638620/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGTAATCAACCATTAAAGCTTGATATCGCTTTTAAAATCGCGGTTTTcttt                                           >  2:819202/1‑137 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAAATCGCGTTTTTCata                                           >  2:1286878/1‑140 (MQ=255)
                                                                                                  taaaTGGTTTTTTTTTTTACATTTTTTTGTTTTTAGAAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           <  2:1409144/138‑1 (MQ=255)
                                                                                                    aaTGTATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  2:803166/140‑1 (MQ=255)
                                                                                                    aaTGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGCAAAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  2:217638/140‑1 (MQ=255)
                                                                                                     aTGGTTTTTTTTTTTAAAATTTTTTTTATTTAGGCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                        <  2:777363/140‑1 (MQ=255)
                                                                                                      tGGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAACAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                       <  1:233258/140‑1 (MQ=255)
                                                                                                            ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                                 >  1:1467442/1‑140 (MQ=255)
                                                                                                                  ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCATATCCTGTATACAGCTGa                           >  2:1437185/1‑140 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  1:453202/140‑1 (MQ=255)
                                                                                                                              ttttttAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa               <  1:296909/137‑1 (MQ=255)
                                                                                                                              gtgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa               <  1:1520453/140‑1 (MQ=255)
                                                                                                                                   tAGGCATATAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt          <  1:631789/140‑1 (MQ=255)
                                                                                                                                    aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc         >  1:568357/1‑140 (MQ=255)
                                                                                                                                    aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc         >  1:409714/1‑140 (MQ=255)
                                                                                                                                    aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc         >  2:982719/1‑140 (MQ=255)
                                                                                                                                         tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTCATATCCTGTATTCAGCTGACGCGGACGGGCAATCTTCATAcc    >  2:743894/1‑140 (MQ=255)
                                                                                                                                           taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt  >  1:837859/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752986‑753264

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: