Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,611,582	A→G	16.0%	S267P (TCC→CCC)	rsmC ←	16S rRNA (guanine(1207)‑N(2))‑methyltransferase RsmC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,611,582	0	A	G	16.0%	170.3 / 2.9	66	S267P (TCC→CCC)	rsmC	16S rRNA (guanine(1207)‑N(2))‑methyltransferase RsmC
Reads supporting (aligned to +/- strand):  ref base A (5/49);  new base G (9/3);  total (14/52)
Fisher's exact test for biased strand distribution p-value = 9.26e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTCGAGAGCAG  >  NC_007779/4611346‑4611819                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    tAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                                    <  2:949835/238‑1 (MQ=255) tAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                                    <  1:268430/238‑1 (MQ=255)                          gTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCGCCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                                    <  1:841640/213‑1 (MQ=255)                                    cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGGTGATAAT‑AT‑CAAAACGA‑CCCTTTACCTCCGGAAAg                                                                                                                                                                                                             >  1:1410024/1‑239 (MQ=255)                                               ttGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGc                                                                                                                                                                                                  <  1:1007444/239‑1 (MQ=255)                                                tGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCa                                                                                                                                                                                                 <  1:595388/239‑1 (MQ=255)                                                      gcCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACt                                                                                                                                                                                           <  1:701881/239‑1 (MQ=255)                                                      gcCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACt                                                                                                                                                                                           <  1:1141037/239‑1 (MQ=255)                                                      gcCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGGTGATCAT‑AT‑CAAAAAGC‑CCTTTTACCTCCGAAAAGACCTTGGTGGCAAAAACt                                                                                                                                                                                           >  1:2076507/1‑239 (MQ=255)                                                       cccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTTGGGTGGAGCAT‑ATACAAACCAA‑CCTTT‑ACCTCGGGAAAGCGGCTGCTgggtcaggcac                                                                                                                                                                                          >  2:105882/1‑230 (MQ=255)                                                       cccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGAGCATAAT‑AAAAAGCA‑CCTTTCCCCCCGTAAAA‑ACGGTGCTgggagaaaatt                                                                                                                                                                                          >  2:701881/1‑230 (MQ=255)                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCAT‑AT‑CAAAAAGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCCAAGACTTc                                                                                                                                                                                         >  1:142659/1‑239 (MQ=255)                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGAGGATCAT‑AT‑CAAAACGA‑CCTCTCCCCTCGGAAAAAACGTTTCTGGCAAAGACtt                                                                                                                                                                                          >  1:1932371/1‑238 (MQ=255)                                                        ccTGGTTGCGCGATCACTTCGTGGAAGCCAAATGTCCCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCTCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                                         <  2:921972/239‑1 (MQ=255)                                                         cagtgtGCGCGATCGCTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTGTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAGAGACTTCa                                                                                                                                                                                        <  2:1067002/234‑1 (MQ=255)                                                           gTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAcc                                                                                                                                                                                      <  1:762120/239‑1 (MQ=255)                                                                     tCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCg                                                                                                                                                                           <  2:1092124/239‑1 (MQ=255)                                                                        cttcttGGAAGCAAAATGTCCCATCCAGCACGTCCGGATAAGGCAGGAAGGCCTCCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTGTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTg                                                                                                                                                                        <  2:467560/239‑1 (MQ=255)                                                                               gAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGAGGATCAT‑AT‑CAAAAAGA‑CCTTTCACCTCCGGAAAAACGGTGCTGGCAAAAACTTCACCCTCAAA‑ACCCTTGGGCGCaa                                                                                                                                                                 >  1:1421946/1‑239 (MQ=255)                                                                                 aGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGAACAT‑AT‑CAAAACGA‑CCCTTCACCCCGGGAAAAAAGGTGGTGGCAAAAAACTAACCCTCAAA‑ACCCTTTGCCGCCAgg                                                                                                                                                               >  1:209423/1‑238 (MQ=255)                                                                                   ccAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCttt                                                                                                                                                                                                                          <  1:1764325/179‑1 (MQ=255)                                                                                       aTGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTTCAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTgcg                                                                                                                                                         <  1:995656/239‑1 (MQ=255)                                                                                          tctcATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAcc                                                                                                                                                                                      <  1:1181743/208‑1 (MQ=255)                                                                                                         ccGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgcc                                                                                                                                       <  1:541494/239‑1 (MQ=255)                                                                                                                aGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                                <  1:502588/239‑1 (MQ=255)                                                                                                                aGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                                <  1:1671882/239‑1 (MQ=255)                                                                                                                 ggCGGGAAGGCGTTCTCTGCAATAGGGAGCTCGCCGCCCCTGTTAAGGTGACGCAGCGCCCCTCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AG‑CAAAACGA‑CCTTCCACCTCGGGAAAGACGTTGCTGGCAAATACTTCACCTTCAAC‑TCCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                                <  2:562471/238‑1 (MQ=255)                                                                                                                  gCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAgatgat                                                                                                                                                                                                                                               <  1:23253/130‑1 (MQ=255)                                                                                                                  gCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgcc                                                                                                                                       <  1:913553/230‑1 (MQ=255)                                                                                                                       aaGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑Ac                                                                                                                           <  1:208316/236‑1 (MQ=255)                                                                                                                       aaGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACa                                                                                                                          <  1:1227959/237‑1 (MQ=255)                                                                                                                           ctttcgcgtaaaTGCGCGGCTAGCCGCCGTTATTAAGAGGACGCTCCGCCCCGCGAATCGGCGTTTGCGCCGCATCCAGGCGCGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGAGGATCAT‑AT‑CAGAACTC‑CCTTGCGCCGCGGAAAAGACGTTGCTGGCAAAGACTTCCCCTACAAC‑ACCGTTGGCCGCAAGGGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGcg                                                                                                                     <  2:472866/228‑2 (MQ=255)                                                                                                                           cGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTGGGGGGTGATCAA‑AT‑CAAAACGACCCTTTCCCCCCCGCAAA‑AAGGTGTTTGCAAAAAATTCCCCTTCCACAACCCTTGTCCCCCA‑TGGTGGCCGGCTGCCTTCTCACCCCGGGGGCGAG‑ACATcccca                                                                                                                    >  2:253426/1‑235 (MQ=255)                                                                                                                                gCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCGCCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGC‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGt                                                                                                               <  1:1116045/239‑1 (MQ=255)                                                                                                                                    cAATTAGCAGGTGCCGGCCGCTATTAAGATGACGCACGGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTgaga                                                                                                           <  2:1270754/239‑1 (MQ=255)                                                                                                                                         cGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑Caaaa                                                                                                                                                                                                                                   <  1:1010254/117‑1 (MQ=255)                                                                                                                                         cGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCgg                                                                                                                                                      <  1:326990/192‑1 (MQ=255)                                                                                                                                                cgccgccacTGGTAAGGTGACGCGACGCTCCGCGAGGCGAGGTTGTGGCCTAATCCAGGCTGGTTTACATCCCATCCTGGAACGGCGGGTTGGGGTTGATAAT‑AG‑CAAAACGA‑CCTTGCGCCTCGGAAGAGAAGTTGCTGGTAATGACTTGTTCTGCAAC‑GCCGTTGGCCGCAAGTGTTGCGCGCCTTGCTTCTACCGCCGGGCCAGAG‑ACATCGCACAAGGTGAGACGAATTGTCgtc                                                                                               <  2:1263732/239‑3 (MQ=255)                                                                                                                                                      cgcTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATg                                                                                          <  1:1130350/238‑1 (MQ=255)                                                                                                                                                         tATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGt                                                                                                               <  1:351410/214‑1 (MQ=255)                                                                                                                                                               gATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCt                                                                                                                                             <  1:19827/179‑1 (MQ=255)                                                                                                                                                                     gcaccgcaccGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCCGGGCAGAG‑ACATCGCACAAGGGGAGACGAATTTTCGGCGAAAGGCGCGCAAAAGC‑a                                                                            >  1:1248919/1‑236 (MQ=255)                                                                                                                                                                       cccgccacgcgAATAAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGGGGCTGGCTTCTGCCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTg                                                                        <  2:318120/231‑1 (MQ=255)                                                                                                                                                                          gcaccgcGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACAACCCACAAGGGGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGaaa                                                                     >  1:1121442/1‑238 (MQ=255)                                                                                                                                                                                         gTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCcac                                                      <  1:105882/238‑1 (MQ=255)                                                                                                                                                                                         gTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCca                                                       <  1:2093065/237‑1 (MQ=255)                                                                                                                                                                                             gcgcCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCg                                                 <  1:1120149/239‑1 (MQ=255)                                                                                                                                                                                             gcgcCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCg                                                 <  1:566600/239‑1 (MQ=255)                                                                                                                                                                                              cgccgcATCAAGGTTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGTTCAT‑AT‑CAAAACGC‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTGGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCg                                                 <  2:142659/238‑1 (MQ=255)                                                                                                                                                                                                  gcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATc                                            <  1:1829100/239‑1 (MQ=255)                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAac‑a                                                                                                                                                                              <  1:307762/109‑1 (MQ=255)                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                                          <  1:1747157/239‑1 (MQ=255)                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                                          <  1:1244011/239‑1 (MQ=255)                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGCGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATcc                                           <  1:1645148/238‑1 (MQ=255)                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGATACTCAT‑AT‑CAAAACGG‑CCTTTCGCCTCGGAAAGGGCGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTGGCGCGGCTGGCTTCTGCCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCCCGCCCGCGCTACAGCCGACATCCa                                          <  2:1108040/239‑1 (MQ=255)                                                                                                                                                                                                     tCCAGGCTGGTTTGCATCCCTTGTGGGAAGGGGGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCGGGCAAAGACTTCACCTTCAAC‑CCCGTTGGCCGCAAGTGTTGCGGGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAg                                         <  2:1248919/239‑1 (MQ=255)                                                                                                                                                                                                       cAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCg                                                                                                                                                                           <  1:1510742/109‑1 (MQ=255)                                                                                                                                                                                                          gCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑GCCGTTGGCCGCAAGTGTTGCGCggctg                                                                                                                                                   <  1:810369/130‑1 (MQ=255)                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGAGCAT‑AT‑CAAAACGC‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCGTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccg                                                                                                                                         <  2:1076057/136‑1 (MQ=255)                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                                          <  1:1063130/229‑1 (MQ=255)                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTa                                  >  1:1504336/1‑237 (MQ=255)                                                                                                                                                                                                                ttGCATCCCATCGTGGAACGGCGGGGTGGGGAAGATAAT‑AT‑TAAAAAAA‑CCCTTCCACCCGGAAAAAAAGGTGCGGGCAAAAACCTCCCCTTCAAC‑ACCCTTGTCCGCCAGGGGTGCGCGGCAGGCGTTTTCCCCCGGGGGCGAGAAAATCCCCCAAGGGGGGACGAATTTTTGCGGAAGTGGGC‑CAAACGCAAACGGAAAGAAC‑CCCGCGCCCCACCCCACCTCCAGccacttacc                                >  2:92284/1‑229 (MQ=255)                                                                                                                                                                                                                  gCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTc                            <  1:937019/239‑1 (MQ=255)                                                                                                                                                                                                                       ccATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg                       <  1:1302175/239‑1 (MQ=255)                                                                                                                                                                                                                       ccATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg                       <  1:92284/239‑1 (MQ=255)                                                                                                                                                                                                                        cATCGTGGAACGGCGGGGTGGAGATGATCAT‑TT‑CAAAACGA‑CCTTTCACCCCGGCAAAGAGGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCCCGGCCGGCTCCTACCGCCGGGGCaga                                                                                                                               >  1:1076057/1‑136 (MQ=255)                                                                                                                                                                                                                             tGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGgcgcgc                                                                                    <  1:1667938/173‑1 (MQ=255)                                                                                                                                                                                                                               gAAGGGGGGGTTGGAGATGATCAT‑AT‑CAAACCGC‑CCTTTCCCCTCGGAAAGGCCGTTGCGGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGtt                <  2:1079356/238‑1 (MQ=255)                                                                                                                                                                                                                                      gggTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTCg         <  1:1110819/238‑1 (MQ=255)                                                                                                                                                                                                                                             agaTGAACAT‑AT‑CAACACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGTCTTCACCTTCAAC‑CCCGGTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTCGAGagcag  <  2:610194/238‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    TAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAT‑AT‑CAAAACGA‑CCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAAC‑ACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGC‑AACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTCGAGAGCAG  >  NC_007779/4611346‑4611819

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