Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,745,509	0	A	T	40.3%	74.9 / 35.7	67	T236S (ACC→TCC)	yiaN	2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base A (13/27);  new base T (26/1);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.11e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCAC  >  NC_000913/3745372‑3745659                                                                                                                                           |                                                                                                                                                           ggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTG‑GTTT‑CCCCTcccca                                                                                                                                                 >  2:587935/1‑146 (MQ=255)        ggCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGcc                                                                                                                                      <  1:920873/151‑1 (MQ=255)                gCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTG‑TTTT‑CCCCTACCCCAG‑CCgggg                                                                                                                                        >  2:704255/1‑141 (MQ=255)                     aaTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCTG‑TTTT‑TCCCCCCCGCAG‑CCgggg                                                                                                                                        >  2:180762/1‑136 (MQ=255)                     aaTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCCG‑TTTT‑CCCCTCCCCAAA‑CCGGGGGCGGTGCTTCttttt                                                                                                                         >  2:1050570/1‑149 (MQ=255)                       tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTCCTCCGGCCTG‑TTTT‑CCCCCACCCCAG‑ACGGGGGCGGTGCTGCttttttt                                                                                                                       >  2:1165716/1‑147 (MQ=255)                       tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTa                                                                                                                        <  1:953429/150‑1 (MQ=255)                        ttcccccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTAt                                                                                                                       <  2:709673/147‑1 (MQ=255)                        ttcccccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTAt                                                                                                                       <  1:1184285/147‑1 (MQ=255)                           ccccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATgcg                                                                                                                    <  1:833654/150‑1 (MQ=255)                           ccccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATgc                                                                                                                     <  1:762865/149‑1 (MQ=255)                            cccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑Ct                                                                                                                 <  2:1139649/151‑1 (MQ=255)                            cccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCCCAG‑CCGGGGCCGGTGCTGCTTTCTTTtcg‑c                                                                                                                  >  2:33935/1‑150 (MQ=255)                              cGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGTTCCCGCCTG‑TTTT‑CCCCCACCGAAG‑ACGGGGGCGTTGCTTGTTTCtttt                                                                                                                      >  2:272447/1‑142 (MQ=255)                              cGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGt                                                                                                               <  2:831722/151‑1 (MQ=255)                                    cAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTc                                                                                                          <  1:686081/150‑1 (MQ=255)                                    cAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCCCCCCGGCCTG‑TTTT‑CCCCCTCCCAAG‑CCCGGGCCGTTGCTGCTTTCTATTCG‑CTGTTTGtt                                                                                                          >  2:903339/1‑149 (MQ=255)                                         aGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCCACCGAAA‑CCGGGGGCGTTGCTGCTTTTTATTCG‑CTGTTTTTCGCcac                                                                                                     >  1:863008/1‑150 (MQ=255)                                          gCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAg                                                                                                   <  1:284963/151‑1 (MQ=255)                                           cAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTT‑CCCCTACCCAAG‑CCCGGGGCGTTGCTGCTTTCTATGCC‑CTTTTTTTTGcccccgt                                                                                                  >  2:524118/1‑146 (MQ=255)                                           cAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGt                                                                                                  <  2:1103166/151‑1 (MQ=255)                                           cAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCCCCGGCCCG‑TTTT‑TCCCTACCCAAG‑CCGGGGGCGGTGGTGCTTTCTTTGCG‑CTGTTTGTCGcccc                                                                                                     >  2:1073545/1‑146 (MQ=255)                                            aaCGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCCG‑TTTT‑TCCCCTCCCAAG‑CCCGGGCCGGTTCTGCTTTTTTTTCG‑CTTTTTTTCGCCCCCGtt                                                                                                 >  2:87609/1‑151 (MQ=255)                                            aaCGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCCCCCGAAG‑CCGGGGCCGTTGCTGCTTTCTTTTCG‑CTGTTTTTTGCCACAGtt                                                                                                 >  1:431998/1‑151 (MQ=255)                                            aaCGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCCTCCGCCCG‑TTTT‑TCCCCACCCAAA‑CCGGGGGCGTTTTTGCTTTCTATGCG‑CTGTTTGTTGCCCCAGtt                                                                                                 >  1:552955/1‑151 (MQ=255)                                             aCGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTTCTCCCGGCTG‑TTTT‑TCCCTACCCAAG‑CCCGGGCCGTTGCTGCTTTCTATGCG‑CCGTTTGt                                                                                                           >  1:1066077/1‑140 (MQ=255)                                                 tGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCCCCGGCCTG‑TTTT‑ACCCTACCGAAG‑CCGGGGCCGGTGCTGCTTTCTTTGCG‑CTGTTTGTCGCCCCAGTTAttt                                                                                             >  1:843/1‑150 (MQ=255)                                                       aaaTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCTCTCCCGCCCG‑TTTT‑TCCCCCCCCCAG‑ACGGGGGCGTTGTTGGTTTCTTTTAC‑CTGTTTGTTGCCCCCGttttttttc                                                                                          >  2:182141/1‑144 (MQ=255)                                                          tCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTTCTCCCGCCCGTTTTT‑CCCCTCCCAAAC‑CCGGGG‑CGGTGCTGCTTTTTTTGCG‑CTGTTTTTTGCCACCGTTTTTT‑TTCGTGgaa                                                                                   >  2:843254/1‑150 (MQ=255)                                                           cTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGGCTC‑TTTT‑CCCCCACCCAAG‑CCCGGGGCGCTGCTGCCTTCTTTGCC‑CCGTTTGTCTCCCCAGTTATTT‑TCCACGaaa                                                                                   >  2:859019/1‑149 (MQ=255)                                                             ggCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTTGTCCGGCCTG‑TTTT‑CCCCTACCGGAA‑ACGGGGCCGGTGTTGCTTTTTCTGCG‑CTGTTTTTGGCCCCGGTTGTTT‑ACCGTGGAATGa                                                                                >  1:132998/1‑150 (MQ=255)                                                             ggCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAAt                                                                                  <  2:579837/148‑1 (MQ=255)                                                                 ctctTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCtt                                                                            <  2:1054555/150‑1 (MQ=255)                                                                     tttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCTCTCCCGGCCG‑TTTA‑CCCCCACCCAAG‑CCGGGGCCGTTGCTGCTTTTTTTGCG‑CTGTTTTTCGCCCCCGGTATTT‑ACCGGGAAAAGAACTTcccc                                                                        >  1:128042/1‑150 (MQ=255)                                                                      ttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTTCTCCCGCCTA‑GTTT‑CCCCCACCGCAA‑ACGGGGGCGTTGCTGCTTTTTTTTCG‑CCGTTTTTCGCCCCAAGTATTTCACCGTAAATTGACTTTccccc                                                                       >  2:1019727/1‑150 (MQ=255)                                                                      ttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGGGCCG‑TTTT‑TCCCTACCGAAA‑CCGGGGGCGGTGGTGCTTTCTTTGCG‑CTGTTTTTTGCCCCAATTATTT‑TTCGTGAAATGACCTTTGCca                                                                       >  1:1067395/1‑150 (MQ=255)                                                                          ctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCTCTCCGGCCTT‑TTTT‑CCCCTACCCAAA‑CCGGGGGAGTTGTTGTTTTTTTTGCG‑CTGTTTGTCGCCCCAGTTTTTT‑TCCGTGAAAAGACCTTTGCCccac                                                                    >  2:762732/1‑149 (MQ=255)                                                                          ctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCTCCGGCCTT‑TTTT‑CCCCTCCCCAAA‑CCCGGGCCGTTGCTGTTTTTTTTGCG‑CCTTTTTTTCCCCCAGTTATTT‑TCCGTGAAAATACCTTCGCCACACtt                                                                  >  2:687866/1‑151 (MQ=255)                                                                          ctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCCCCGGCCTG‑TTTT‑CCCCCCCCGAAGCCCGGGGCCTTTG‑TGCTTTCTTTGCC‑CTGTTTTTTGCCCCCGGTTTTT‑TCCGTGGAAAGACCCTTGCCCCCCtt                                                                  >  2:380143/1‑151 (MQ=255)                                                                           tcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCCGCCCG‑TTTT‑CCCCTACCCAAA‑CCGGGGCCGTTGCTGCTTTCTTTTCG‑CTGTTTTTCGCCCCAGCTATTT‑ACCGTGAAATGACCTTTGCCCCACttt                                                                 >  2:795212/1‑151 (MQ=255)                                                                            ccGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACttt                                                                 <  1:1050620/150‑1 (MQ=255)                                                                            ccGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACtt                                                                  <  1:440541/149‑1 (MQ=255)                                                                                tATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTGTTTGTCCCGCCCG‑TTTATCCCCCACC‑AAA‑CCGGGGGCGGTGGTGTTTTCTTTTCGCCGGTTTTTCTCC‑CAGTTATTT‑TCCGTGAAAATAACTTTGCCCCACcttt                                                                >  2:698603/1‑146 (MQ=255)                                                                                tATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTTCTCCCGGCTG‑TTTA‑CCCCCCCCGAAA‑CCGGGGCCGTTTCTGTTTTCTATTCG‑CTGGTTTTTTCCACAGGTTTTT‑AACGGGAAATGACCTTTGCCACACttt                                                                 >  1:1084066/1‑146 (MQ=255)                                                                                 aTCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCAt                                                            <  1:826139/150‑1 (MQ=255)                                                                                     gggCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACt                                                       <  2:870243/151‑1 (MQ=255)                                                                                       gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCa                                                     <  2:380955/151‑1 (MQ=255)                                                                                          cTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCAtt                                                   <  2:687947/150‑1 (MQ=255)                                                                                            gTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTg                                                  <  2:272671/149‑1 (MQ=255)                                                                                            gTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCCG‑TTTA‑CCCCCTCCTAAG‑CCGGGGCCGCTGCTGCTTTTTTTGGG‑CTGTTTGTCCCCACAGCTACTT‑ACCGTGAAATGACCCTTTCCACACTTTTGCATGTAccc                                                      >  1:91574/1‑143 (MQ=255)                                                                                             tttCTTCCAGTCATCATTATTGGCGGGTTTTCCCCCGCCCCG‑TTTT‑ACCCTACCCCAG‑GCGGGGGCGGTGTTGCTTTCTTTGCG‑CTGTTTTTTGCCACCGTTATTT‑TTCCTGAAATGAACTTCGCCCCACCTTTGCATGTatccaat                                                   >  1:526063/1‑141 (MQ=255)                                                                                             tttCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGgc                                                <  1:651323/150‑1 (MQ=255)                                                                                             tttCTTCCAGTCATCATTATTGGCGGGGTTTCCTCCCGCCTG‑TTTT‑CCCCTACCCAAT‑CCGGGGGCGCTTCTGCTTTCTTTGGG‑CTGTTTTTTGCCACAGTTATTT‑TCCGTGGAATGACCTTCGCCCCACTTTTGCATGTACTCAtt                                                   >  2:1120723/1‑147 (MQ=255)                                                                                                  tcCAGTCATCATTATTGGTGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCgcgg                                            <  1:1141568/149‑1 (MQ=255)                                                                                                     aGTCATCATTATTGGCGGGGTTTGCTCCGGCCTG‑GTTT‑CCCCTACCCAAG‑ACGGGGCCGGTGCTGCTTTCTTTGCG‑CTGTTTTTTGCCCCAGTTATTT‑ACCGTGAAATGACCTTCGCCCCCCTTTTGCATGTACTCCTTTGCgcgg                                            >  2:1173752/1‑146 (MQ=255)                                                                                                        cctcatTATTGTCGGGTTTCGCTTCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGGGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAac                                     <  2:1187495/148‑1 (MQ=255)                                                                                                              tattGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTc                               <  1:1252404/150‑1 (MQ=255)                                                                                                               attGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAg                             <  1:48995/151‑1 (MQ=255)                                                                                                                 tGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGGGAAATGAACTTCGCCACACTTTGGCATGTACTCATTGGCGGGGCGAAAACCACCCCAg                             >  2:869094/1‑149 (MQ=255)                                                                                                                    cGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTATGTCGCCACAGTTATTT‑GCCGTGAAATGACCGTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAGACCACCTCAGTGGtt                        <  2:769608/151‑1 (MQ=255)                                                                                                                     gggTTTCGCTCCGGCCTG‑TTTT‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTa                       <  2:622210/151‑1 (MQ=255)                                                                                                                      ggTTTTGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTAt                      >  1:1009678/1‑151 (MQ=255)                                                                                                                       gTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATg                     >  2:1179516/1‑151 (MQ=255)                                                                                                                             cTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCt                <  2:802205/150‑1 (MQ=255)                                                                                                                                cGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCtggtg            >  2:971437/1‑151 (MQ=255)                                                                                                                                    cTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCt        >  2:47711/1‑151 (MQ=255)                                                                                                                                      g‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCg      >  1:1036090/1‑151 (MQ=255)                                                                                                                                           a‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGcac  >  2:1162562/1‑151 (MQ=255)                                                                                                                                           |                                                                                                                                                           GGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTG‑TTTA‑CCCCTACCGAAG‑CCGGGGCCGTTGCTGCTTTCTATGCG‑CTGTTTGTCGCCACAGTTATTT‑ACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCAC  >  NC_000913/3745372‑3745659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑