Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	729,695	T→G	23.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	23.4%	86.4 / 3.7	54	V38G (GTG→GGG)	rhsC	Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base G (11/3);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 3.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGC  >  NC_000913/729475‑729933                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 tttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGAc                                                                                                                                                                                                                                >  1:1152561/1‑239 (MQ=255)  ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACg                                                                                                                                                                                                                               >  1:703702/1‑239 (MQ=255)        aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGGCGGCCgg                                                                                                                                                                                                                          >  1:323673/1‑238 (MQ=255)                 aTCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGTCCGGCCCTCCGGTc                                                                                                                                                                                                                >  1:199540/1‑239 (MQ=255)                 aTCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGAATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGGTCGGGGTGGCCCGGGGGGGGGGCGGCCGGCCATCCgggc                                                                                                                                                                                                                >  2:1065437/1‑237 (MQ=14)                      aCAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCCGGCGGGGGACGCATCGGTGCCCCCCCCCGCGTGGCCCGTTCGGGGGGCCCCGGCGGGGGGGCGGCCCGCCCTcccggccctt                                                                                                                                                                                                            >  2:124016/1‑230 (MQ=14)                      aCAGATATCTCACTTAAGGCTTTCTTATAAATCTGCAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGTTGAGCGGAAAACCGGCGACGCGTCAGGGCGACATGCCGCAGTATCGCGG‑AGCATTGTCAAGGGTTCAGCGGGGGTACACCCCGGTGCCCCCACCGGCGTGGCCCGTTCGGTCGGCCCCCGCGGGGGGGCGCCCCGCCCTCCGGTCTAT‑cc                                                                                                                                                                                                          >  2:582675/1‑238 (MQ=14)                                    tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGGGGCCCGTTCGGGGTGCCCCGGCGGGGGGGCGTCCGGCCCTCCCGGCAAT‑CCGCCGCTCGGGGCaa                                                                                                                                                                                            >  1:429386/1‑239 (MQ=14)                                    tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCaa                                                                                                                                                                                            <  2:1152561/239‑1 (MQ=255)                                     aaGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGTTCAGCCGGGGGACGCCTCGGGGCCCCCCCCGGCGTGGCCTTTTCGGGGTGCCCCGCCGGGGTGACGTCCGGCCATCCGGGCAAT‑CCGCTGCCCGGTGCaaa                                                                                                                                                                                           >  1:322692/1‑239 (MQ=14)                                       ggCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTTAGCCGGGGAACGCATCGGGGCCCCCACCGGCGTGGACTGTTCGGGGTGGCCCGGCGGGGGGGAGGCCCGCCATCCCGTCCAT‑CCCCCTCTCgggggcaagg                                                                                                                                                                                         >  1:1164467/1‑231 (MQ=14)                                           ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGCGGCCTGTTCGGTGGGGCCCCGCGGGGGGGCGTCCGGCCCTCCCGCCCAA‑TCGCTGCCCGGTGCaaacggc                                                                                                                                                                                       >  1:517226/1‑233 (MQ=14)                                           ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGGGGGGCGGCCGGCCATCCCGGCAAT‑CCGCCGCCCCGTGCCAAGGGCCt                                                                                                                                                                                     >  1:275946/1‑239 (MQ=14)                                            tCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCtt                                                                                                                                                                                    <  1:1065437/239‑1 (MQ=255)                                                       tGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCCGGGTTCAGCCGGGGGCGCCACCGGGGCCCCCCCCGGCGGGGCCCTTTCGTGGTGCCCCTGGGGGGGGACGCCCCCCCCCCCGGGCCAT‑CCCCGCGCCGGGGCACAGGGCCTTccccggggaa                                                                                                                                                                          >  2:820634/1‑231 (MQ=2)                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCCTCCGGGCAAT‑CCTCTGCTCGGGGAACAGGGGCTTCCCGGGGaaaac                                                                                                                                                                        >  2:1047796/1‑237 (MQ=14)                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGTGCGTCCGGACCTCCCGGCAAT‑CCGCTGCTCGGGGCAAAGGGCCTTTCCCGTGAAAc                                                                                                                                                                         >  1:786178/1‑238 (MQ=14)                                                                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGGCGTCCGGCCCTCCGGGCCAT‑CCGCCGGTCGGGGCAAAGGGCCTTCCCGGGGAAACCGGCATCcccc                                                                                                                                                              >  1:902813/1‑238 (MQ=14)                                                                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGAAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCCGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGGGGGGGGCCGTCCCGCCAACCCGCCAAT‑CCCCTCCTCGGGGCCAAGGGCCTTCCCGGTGAAACCGAcctgcccc                                                                                                                                                              >  2:258906/1‑231 (MQ=14)                                                                       tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCCGGGGACGCATCCGTTGCCCCACCGGCGTGGCCTGTTCCGCTGGCCCCGGCGGGGGGGCGTCCGGCCCTACCGCCCAT‑CCGCGGCCCGGGGCACAGGCCCGTCCCGGTGGAAACCACATCCCCCTgcc                                                                                                                                                          >  2:1071136/1‑238 (MQ=14)                                                                       tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGCAGCATTTTCCAGGGTTCAGCCGGGGTACGCACCGGGGCCCCCACCGGCGCGGCCTGTTCGGGGGGCCCCGGCGGGGGGGCGTCCGGCCCTCCCGGCCAT‑CCCCCGCCCCGGGCCAAGGTCCTTCCCGGGGGACACGCCAccgccccgcc                                                                                                                                                          >  2:651648/1‑238 (MQ=14)                                                                                    aaCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGt                                                                                                                                            <  2:74878/239‑1 (MQ=255)                                                                                      cccGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGGCGGCCGGCCATCCGGTCAAT‑CCGCCGCTCGGTGCAAAGGGCCTTCCCGGTGGAACCGACATCGCCCCGCCCCGGCCCCTGCCCtt                                                                                                                                           >  1:1130031/1‑238 (MQ=14)                                                                                           tAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATtct                                                                                                                                     <  1:461157/239‑1 (MQ=255)                                                                                                cAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGGTCGGTGGGCCCCCGCGGGGGGGCGTCCGGCCCACCCGGCCAT‑CCGCCGCCCGGGGCCAAGGGCCCTGCCCGGGAAACCCACATCCCCCCTCCCCGCCCCCTGCCCTTCCCTCTCTcc                                                                                                                                 >  2:774653/1‑238 (MQ=14)                                                                                                        aGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGGCGTCCGGCCATCCCGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTTTTCTCCCCCACCTAc                                                                                                                        >  1:905534/1‑239 (MQ=14)                                                                                                                                     cAGGGCGACATGGCGCAGGATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCgcgc                                                                                           <  1:820634/238‑1 (MQ=255)                                                                                                                                           gACATGACGCAGTATGGCGGGAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTgg                                                                                     <  2:703702/238‑1 (MQ=255)                                                                                                                                             cATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCCGTTCGGGGGGCCCCGGCGGGGTGGCGTCCGGGCATCCCGGCAATCCCCCGGCCGGGGGCAAAGGCCTTCCCCGGGGAACCCAACA‑CCCCCCGCCCCGCCCCCTGCCCTTTCTTTTCTCCCCCCCCCCACGCAGTTTACCGGCCAAAAACGCCC‑CTCCCGTggggg                                                                                  >  2:1206995/1‑238 (MQ=255)                                                                                                                                             cATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGGGGCCCGTTCGGTGTGCCCCGGCGGGGGGGCGGCCGGCCCTCCCGGCAAT‑CCCCCGCTCGGGGCAAAGGGCCCTCCCCGGGGAACCCACATCCCCCCGCCCGGCCCCCTGCCGCTCCTTCTCTCCCCCCCCCACAGCAG‑TTACCGGGCCAAAACGCCCGCGCCCggggggg                                                                                  >  2:900048/1‑238 (MQ=255)                                                                                                                                                  cGCAGTACGGCGGTGGCATTGTCCAGGGTTCAGCCGGGGTCCGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGCGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAACCCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCt                                                                              <  2:517226/238‑1 (MQ=255)                                                                                                                                                                        ccAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTg                                                       <  1:124016/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  1:900048/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  1:878114/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  1:774653/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  1:651648/239‑1 (MQ=37)                                                                                                                                                                               tCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                <  2:1176948/239‑1 (MQ=37)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCCG‑TTACCGGCCCAAAACCCCCGCGCCGGGGGGGGGCCTCGGCCCCCGCTGGAAAATGCCCGCgggtctc                                               >  2:54128/1‑233 (MQ=11)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGGCCTTCCCGGTGCAAACGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCCG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGCAAATGCCTGCCGataac                                               >  2:305124/1‑236 (MQ=255)                                                                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCCGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACCACTCCCTGCCCGGCCCGCTCCCGTTCATTCTCTCCCGCCCCTCCAACCG‑TTCACCGACAAAAACCCCCCCCCCGGTCGGGAGCCTCCGCCCCGGCGGGAAAATGCCTGCGGATATc                                               >  2:446632/1‑238 (MQ=11)                                                                                                                                                                                        gTACGCATCGGTGCCCCCACCGGCGTGCCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCAGCCATCCGGACAAA‑CCGCTGCTCGGCGCAAAGGGACTTCCCGGCGAAACCGACAACGCCCTGCCCGGCCCGCCACCCTTCCTTCGCCCCCGCCCCTACATCAG‑TTCCCGCGCAAAAACGCCCCCGCCGGCGGGGCGACTGGGCCCCGGGTGGGCAATGCATGCGGactgccccttgt                                        >  2:961750/1‑227 (MQ=255)                                                                                                                                                                                            gCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTg                                   <  1:1047796/239‑1 (MQ=35)                                                                                                                                                                                            gCAGCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTg                                   <  1:661519/239‑1 (MQ=21)                                                                                                                                                                                             acagggTGCCCCCACGGGCGTGGCCTGTTCGGGGTGCCCAGGCTGGGTGACGTCCGGCCAGCCGGTCAAC‑CCGCTGGTCGGTGCAAAGGTCCTTCCCGGTGAGACCGACATCGCCCTGCCAAGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGGGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTg                                   <  2:832383/234‑1 (MQ=11)                                                                                                                                                                                              tgcGGGGCCCCCACCGGCGGGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcg                                 <  2:516107/237‑1 (MQ=21)                                                                                                                                                                                                   tGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGGGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGCGCTCCGCTTCCAGCTGCGCGGTa                            >  1:616177/1‑239 (MQ=11)                                                                                                                                                                                                    gCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATaa                           <  1:788305/239‑1 (MQ=33)                                                                                                                                                                                                     ccccccccGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATaa                           <  1:305124/238‑1 (MQ=18)                                                                                                                                                                                                     cccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATaa                           <  2:275946/238‑1 (MQ=33)                                                                                                                                                                                                      ccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAaca                         <  1:414813/239‑1 (MQ=33)                                                                                                                                                                                                        ccACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacaca                       <  2:199540/239‑1 (MQ=18)                                                                                                                                                                                                        caaCCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacaca                       <  1:446632/237‑1 (MQ=18)                                                                                                                                                                                                                gTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt               <  1:1071136/239‑1 (MQ=33)                                                                                                                                                                                                                 tGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt               <  1:926453/238‑1 (MQ=33)                                                                                                                                                                                                                             tgtgCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcggc  >  1:1059852/1‑239 (MQ=33)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 TTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAG‑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 < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑