Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,047,936	T→C	17.7%	T295A (ACC→GCC)	yccM ←	putative 4Fe‑4S membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,047,936	0	T	C	17.7%	41.2 / 2.3	22	T295A (ACC→GCC)	yccM	putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base C (2/2);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCAT  >  CP009273/1047684‑1048160                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 caaaatATTGGGCATTTGAGACCACTACCCCGTGCTGGTTGATGTTCCGCTTTTTGCTAATTTAGTGCCCAATCATTGGCCCTTTTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCATAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACAt                                                                                                                                                                                                            <  2:5361/275‑1 (MQ=255)                 tgtgAACATCACCCCGTGCGAGGTGATGTTCAGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTAAGGATACCGGACATGACTCTACGACAGTCa                                                                                                                                                                                            >  1:56973/1‑275 (MQ=255)                                               ccGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGt                                                                                                                                                                                                                                  <  1:161357/207‑1 (MQ=255)                                                                                                                ccGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGAAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATa                                                                                                                                                                                          >  2:404815/1‑182 (MQ=255)                                                                                                                    cACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGc                                                                                                                                                                                      <  1:404820/182‑1 (MQ=255)                                                                                                                            tttGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCg                                                                                                                                                                                    <  1:136683/176‑1 (MQ=255)                                                                                                                                     cTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTTCGAATTTGCCACAGTCGAGAAAACTTTCGGCATTGAGCCGAATc                                                                        >  1:19821/1‑275 (MQ=255)                                                                                                                                                  tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATATTGAGGCTACCGGACATGACTCTACGCAAGTCATACa                                                                                                                                                                                        >  2:22762/1‑150 (MQ=255)                                                                                                                                                      tAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATATTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCg                                                                                                                                                                                    <  1:22762/150‑1 (MQ=255)                                                                                                                                                              cACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTc                                                                                                                                                                                             >  2:361261/1‑133 (MQ=255)                                                                                                                                                                  gCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATAc                                                                                                                                                                                         <  1:361266/133‑1 (MQ=255)                                                                                                                                                                                             gTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACAt                                                                                                                                                                              >  1:159207/1‑117 (MQ=255)                                                                                                                                                                                              tCAGTACCCAGCCCGACAACGCAGAGGTTTTTTTATTGGCGCCATTCTTTTCAAGTGAAAATGTCAATGTTTCGGCTACCGGACATTACTCTACGCAAGTCATACAGCCGGTAAATTCCTCTGTTCTGACCTGAATTAATTTATCGACCGGGATCCTGGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATaa               <  2:19821/275‑1 (MQ=255)                                                                                                                                                                                                 ttAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCa                                                                                                                                                                          <  2:159204/117‑1 (MQ=255)                                                                                                                                                                                                         ccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATACCACTGTTCGTACCTGAATTAATTTATc                                                                                                                                                 >  2:304111/1‑134 (MQ=255)                                                                                                                                                                                                         ccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGCCAATGTTGAg                                                                                                                                                                                                                        >  2:175641/1‑63 (MQ=255)                                                                                                                                                                                                             aacaagcaanttcTTTTTTATGTGACTCAGGTGTTTGCAGTGAAAAGGTCAATTTTGAGGCTACCGGACATGACTCTACGCTAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACGAGTATCCGTGATGGGCAATTGTTTGCGCATTTGCGACAGTGGATACAACTTTAGGCAGTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCAGCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCAt  <  2:56973/261‑1 (MQ=255)                                                                                                                                                                                                             aCAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATACCACTGTTCGTACCTGAATTAATTTATCGAcc                                                                                                                                             <  1:304116/134‑1 (MQ=255)                                                                                                                                                                                                             aCAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGCCAATGTTGAGGCTa                                                                                                                                                                                                                    <  1:175645/63‑1 (MQ=255)                                                                                                                                                                                                              cAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGCCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGaattaat                                                                                                                                                      >  2:476409/1‑124 (MQ=255)                                                                                                                                                                                                                  gCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGCCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCCCTGTTCGTCCCTGAATTAATTTAt                                                                                                                                                  <  1:476417/124‑1 (MQ=255)                                                                                                                                                                                                                           tttttATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATg                                                                                                                                                                                                                             >  1:637440/1‑40 (MQ=255)                                                                                                                                                                                                                               tATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGa                                                                                                                                                                                                                         <  2:637429/40‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 CAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCAT  >  CP009273/1047684‑1048160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑