Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	729,695	T→G	22.4%	V38G (GTG→GGG)	rhsC →	Rhs protein with putative toxin domain, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,695	0	T	G	22.4%	90.1 / 3.3	56	V38G (GTG→GGG)	rhsC	Rhs protein with putative toxin domain, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (19/22);  new base G (11/3);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 6.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGC  >  NC_000913/729475‑729933                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  tttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGAc                                                                                                                                                                                                                                 >  2:1152561/1‑239 (MQ=255)  ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACg                                                                                                                                                                                                                                >  2:703702/1‑239 (MQ=255)        aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGGCGGCCgg                                                                                                                                                                                                                           >  2:323673/1‑238 (MQ=255)                 aTCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGTCCGGCCCTCCGGTc                                                                                                                                                                                                                 >  2:199540/1‑239 (MQ=255)                 aTCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGAATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGGTCGGGGTGGCCCGGGGGGGGGGCGGCCGGCCATCCgggc                                                                                                                                                                                                                 >  1:1065437/1‑237 (MQ=14)                      aCAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCCGGCGGGGGACGCATCGGTGCCCCCCCCCGCGTGGCCCGTTCGGGGGGCCCCGGCGGGGGGGCGGCCCGCCCTcccggccctt                                                                                                                                                                                                             >  1:124016/1‑230 (MQ=14)                      aCAGATATCTCACTTAAGGCTTTCTTATAAATCTGCAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGTTGAGCGGAAAACCGGCGACGCGTCAGGGCGACATGCCGCAGTATCGCGG‑AGCATTGTCAAGGGTTCAGCGGGGGTACACCCCGGTGCCCCCACCGGCGTGGCCCGTTCGGTCGGCCCCCGCGGGGGGGCGCCCCGCCCTCCGGTCTAT‑cc                                                                                                                                                                                                           >  1:582675/1‑238 (MQ=14)                                    tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGGGGCCCGTTCGGGGTGCCCCGGCGGGGGGGCGTCCGGCCCTCCCGGCAAT‑CCGCCGCTCGGGGCaa                                                                                                                                                                                             >  2:429386/1‑239 (MQ=14)                                    tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCaa                                                                                                                                                                                             <  1:1152561/239‑1 (MQ=255)                                     aaGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGTTCAGCCGGGGGACGCCTCGGGGCCCCCCCCGGCGTGGCCTTTTCGGGGTGCCCCGCCGGGGTGACGTCCGGCCATCCGGGCAAT‑CCGCTGCCCGGTGCaaa                                                                                                                                                                                            >  2:322692/1‑239 (MQ=14)                                       ggCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTTAGCCGGGGAACGCATCGGGGCCCCCACCGGCGTGGACTGTTCGGGGTGGCCCGGCGGGGGGGAGGCCCGCCATCCCGTCCAT‑CCCCCTCTCgggggcaagg                                                                                                                                                                                          >  2:1164467/1‑231 (MQ=14)                                           ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGCGGCCTGTTCGGTGGGGCCCCGCGGGGGGGCGTCCGGCCCTCCCGCCCAA‑TCGCTGCCCGGTGCaaacggc                                                                                                                                                                                        >  2:517226/1‑233 (MQ=14)                                           ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGGGGGGCGGCCGGCCATCCCGGCAAT‑CCGCCGCCCCGTGCCAAGGGCCt                                                                                                                                                                                      >  2:275946/1‑239 (MQ=14)                                            tCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCtt                                                                                                                                                                                     <  2:1065437/239‑1 (MQ=255)                                                       tGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCCGGGTTCAGCCGGGGGCGCCACCGGGGCCCCCCCCGGCGGGGCCCTTTCGTGGTGCCCCTGGGGGGGGACGCCCCCCCCCCCGGGCCAT‑CCCCGCGCCGGGGCACAGGGCCTTccccggggaa                                                                                                                                                                           >  1:820634/1‑231 (MQ=2)                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCCTCCGGGCAAT‑CCTCTGCTCGGGGAACAGGGGCTTCCCGGGGaaaac                                                                                                                                                                         >  1:1047796/1‑237 (MQ=14)                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGTGCGTCCGGACCTCCCGGCAAT‑CCGCTGCTCGGGGCAAAGGGCCTTTCCCGTGAAA‑c                                                                                                                                                                         >  2:786178/1‑238 (MQ=14)                                                                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGGCGTCCGGCCCTCCGGGCCAT‑CCGCCGGTCGGGGCAAAGGGCCTTCCCGGGGAAA‑CCGGCATCcccc                                                                                                                                                              >  2:902813/1‑238 (MQ=14)                                                                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGAAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCCGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGGGGGGGGCCGTCCCGCCAACCCGCCAAT‑CCCCTCCTCGGGGCCAAGGGCCTTCCCGGTGAAA‑CCGAcctgcccc                                                                                                                                                              >  1:258906/1‑231 (MQ=14)                                                                       tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCCGGGGACGCATCCGTTGCCCCACCGGCGTGGCCTGTTCCGCTGGCCCCGGCGGGGGGGCGTCCGGCCCTACCGCCCAT‑CCGCGGCCCGGGGCACAGGCCCGTCCCGGTGGAA‑ACCACATCCCCCTgcc                                                                                                                                                          >  1:1071136/1‑238 (MQ=14)                                                                       tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGCAGCATTTTCCAGGGTTCAGCCGGGGTACGCACCGGGGCCCCCACCGGCGCGGCCTGTTCGGGGGGCCCCGGCGGGGGGGCGTCCGGCCCTCCCGGCCAT‑CCCCCGCCCCGGGCCAAGGTCCTTCCCGGGGGAC‑ACGCCAccgccccgcc                                                                                                                                                          >  1:651648/1‑238 (MQ=14)                                                                                    aaCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGt                                                                                                                                            <  1:74878/239‑1 (MQ=255)                                                                                      cccGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGGCCGGCCATCCGGTCAAT‑CCGCCGCTCGGTGCAAAGGGCCTTCCCGGTGGAA‑CCGACATCGCCCCGCCCCGGCCCCTGCCCtt                                                                                                                                           >  2:1130031/1‑238 (MQ=14)                                                                                          aTAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCCGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGTCCGGCCACCCCGGCAAT‑CCCGCGCTCGGGGCCAAGGGCCTTCCCGGGCAAACCCCACTTCCCCCCccccggccccgcccgctccat                                                                                                                                        >  1:414813/1‑226 (MQ=14)                                                                                           tAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATtct                                                                                                                                     <  2:461157/239‑1 (MQ=255)                                                                                                cAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGGTCGGTGGGCCCCCGCGGGGGGGCGTCCGGCCCACCCGGCCAT‑CCGCCGCCCGGGGCCAAGGGCCCTGCCCGGGAAA‑CCCACATCCCCCCTCCCCGCCCCCTGCCCTTCCCTCTCTcc                                                                                                                                 >  1:774653/1‑238 (MQ=14)                                                                                                        aGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGGCGTCCGGCCATCCCGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTTTTCTCCCCCACCTAc                                                                                                                        >  2:905534/1‑239 (MQ=14)                                                                                                                                     cAGGGCGACATGGCGCAGGATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCgcgc                                                                                           <  2:820634/238‑1 (MQ=255)                                                                                                                                           gACATGACGCAGTATGGCGGGAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTgg                                                                                     <  1:703702/238‑1 (MQ=255)                                                                                                                                             cATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCCGTTCGGGGGGCCCCGGCGGGGTGGCGTCCGGGCATCCCGGCAATCCCCCGGCCGGGGGCAAAGGCCTTCCCCGGGGAAC‑CCAACA‑CCCCCCGCCCCGCCCCCTGCCCTTTCTTTTCTCCCCCCCCCCACGCAGTTTACCGGCCAAAAACGCCC‑CTCCCGTggggg                                                                                  >  1:1206995/1‑238 (MQ=255)                                                                                                                                             cATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGGGGCCCGTTCGGTGTGCCCCGGCGGGGGGGCGGCCGGCCCTCCCGGCAAT‑CCCCCGCTCGGGGCAAAGGGCCCTCCCCGGGGAA‑CCCACATCCCCCCGCCCGGCCCCCTGCCGCTCCTTCTCTCCCCCCCCCACAGCAG‑TTACCGGGCCAAAACGCCCGCGCCCggggggg                                                                                  >  1:900048/1‑238 (MQ=255)                                                                                                                                                  cGCAGTACGGCGGTGGCATTGTCCAGGGTTCAGCCGGGGTCCGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGCGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAC‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCt                                                                              <  1:517226/238‑1 (MQ=255)                                                                                                                                                                        ccAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTg                                                       <  2:124016/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  2:651648/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  2:878114/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  2:774653/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  2:900048/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  1:1176948/239‑1 (MQ=37)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCCG‑TTACCGGCCCAAAACCCCCGCGCCGGGGGGGGGCCTCGGCCCCCGCTGGAAAATGCCCGCgggtctc                                               >  1:54128/1‑233 (MQ=11)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGGCCTTCCCGGTGCAA‑ACGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCCG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGCAAATGCCTGCCGataac                                               >  1:305124/1‑236 (MQ=255)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGGCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGCTCATTCTCTCCCGCCCCCACAGCGG‑TTCCCGGACCAAAACGCCCCCCCCGGCTGGGTGCCTCCGGCCCTGCCGGaccatccctcctgtgatc                                               >  1:870381/1‑221 (MQ=11)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCCGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACCACTCCCTGCCCGGCCCGCTCCCGTTCATTCTCTCCCGCCCCTCCAACCG‑TTCACCGACAAAAACCCCCCCCCCGGTCGGGAGCCTCCGCCCCGGCGGGAAAATGCCTGCGGATATc                                               >  1:446632/1‑238 (MQ=11)                                                                                                                                                                                        gTACGCATCGGTGCCCCCACCGGCGTGCCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCAGCCATCCGGACAAA‑CCGCTGCTCGGCGCAAAGGGACTTCCCGGCGAAA‑CCGACAACGCCCTGCCCGGCCCGCCACCCTTCCTTCGCCCCCGCCCCTACATCAG‑TTCCCGCGCAAAAACGCCCCCGCCGGCGGGGCGACTGGGCCCCGGGTGGGCAATGCATGCGGactgccccttgt                                        >  1:961750/1‑227 (MQ=255)                                                                                                                                                                                            gCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTg                                   <  2:1047796/239‑1 (MQ=35)                                                                                                                                                                                            gCAGCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTg                                   <  2:661519/239‑1 (MQ=21)                                                                                                                                                                                             acagggTGCCCCCACGGGCGTGGCCTGTTCGGGGTGCCCAGGCTGGGTGACGTCCGGCCAGCCGGTCAAC‑CCGCTGGTCGGTGCAAAGGTCCTTCCCGGTGAGA‑CCGACATCGCCCTGCCAAGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGGGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTg                                   <  1:832383/234‑1 (MQ=11)                                                                                                                                                                                              tgcGGGGCCCCCACCGGCGGGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcg                                 <  1:516107/237‑1 (MQ=21)                                                                                                                                                                                                   tGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGGGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGCGCTCCGCTTCCAGCTGCGCGGTa                            >  2:616177/1‑239 (MQ=11)                                                                                                                                                                                                    gCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATaa                           <  2:788305/239‑1 (MQ=33)                                                                                                                                                                                                     ccccccccGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATaa                           <  2:305124/238‑1 (MQ=18)                                                                                                                                                                                                     cccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATaa                           <  1:275946/238‑1 (MQ=33)                                                                                                                                                                                                      ccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAaca                         <  2:414813/239‑1 (MQ=33)                                                                                                                                                                                                        ccACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacaca                       <  1:199540/239‑1 (MQ=18)                                                                                                                                                                                                        caaCCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacaca                       <  2:446632/237‑1 (MQ=18)                                                                                                                                                                                                                gTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt               <  2:1071136/239‑1 (MQ=33)                                                                                                                                                                                                                 tGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt               <  2:926453/238‑1 (MQ=33)                                                                                                                                                                                                                             tgtgCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcggc  >  2:1059852/1‑239 (MQ=33)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  TTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAA‑CCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGC  >  NC_000913/729475‑729933

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