Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,636,119:1	+A	11.4%	coding (85/1014 nt)	rodZ ←	MreB assembly cytoskeletal protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,636,119	1	.	A	11.4%	113.3 / 11.6	35	coding (85/1014 nt)	rodZ	MreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base . (16/15);  new base A (2/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACAGTATAAGACCGTCTGTGTGC  >  CP009273/2635846‑2636319                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                         gggcaccggcgNNNGTCATCAGGCNNCCGTCGCGTTTTTTTCGGCCTTTCCCGAGGGGAAAACTCTGCATCGGCGCACCTTTTGCAGCCCGCAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCATTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCGTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACt                                                                                                                                                                                                           <  2:507430/264‑1 (MQ=255)                   tcaGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAAGGTTGAAGCAAGATCGGCGGNGGCCTTATCTTCTTCAATGTCGGGGACCGTGGAAAACTTTAGGGCAAGTCGCTCGGCAACGGCCTGCGGAAA‑AAGGCATAgttctgtctct                                                                                                                                                                                       >  1:507440/1‑267 (MQ=255)                                                taagaaacagAAAAACTCTGCATCGGCGCTACGTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCCGACGCGCATAAGAGCGGATATATCCTCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTaa                                                                                                                                                           <  1:324268/265‑1 (MQ=255)                                                                            tagcatgtgaatcccTAAGTGGATCCTGTTTTTCCAGCCCTGGCATCAGTTCGTCTTCTGGAATATGTACCATACGCGCAGAAGAGCGGATATATCCGCGTAGGAATGTTGAAGCAAGATCGGCGGGTTCCTTATCTTCTTCAATGTAGCGTACCGTGGAAACCTGCAGGCAAAGTCGCTTGGCAACGGCCTGCTGACTAAAGTCCTAGTTGTTCGCGAGCATTACGCAGGGGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCtt                                                                                                                                <  2:217240/262‑1 (MQ=255)                                                                            gtcagattGCAGACCCGAGGGGCGCCTGTTTTTCCCGCCCTGGCAGCAGTTTTTCTTCTGGAATAGGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAAGGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTAGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCtt                                                                                                                                <  2:117430/268‑1 (MQ=255)                                                                                             aGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGCTATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTAGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAAGTCCTAGTTGTTCGCGAGCATTACGCCGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTTAGTGGCTTCAGTattaattctgtctc                                                                                                              >  1:217244/1‑269 (MQ=255)                                                                                             aGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTAGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAAGTCCTAGTTGTTCGCGAGCATTACGAAGGCGAGCGCCGGTAGTAAGTGCTTCCTTTTGGGtgtgcgtg                                                                                                                                    >  1:117431/1‑253 (MQ=255)                                                                                                                     ggCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑a                                                                                                                                                                                                         >  2:162420/1‑159 (MQ=255)                                                                                                                         gcagTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTc                                                                                                                                                                                                     <  1:162423/159‑1 (MQ=255)                                                                                                                                                         gcATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑a                                                                                                                                                                                                         <  2:490017/123‑1 (MQ=255)                                                                                                                                                                  cGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTTAGTATTCCTTCGCTGCCGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATTCACACCCACAACCCCAAACAACCgttaa                                        >  1:35636/1‑275 (MQ=255)                                                                                                                                                                   ggATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAgg                                                                                                                                                                           >  2:516127/1‑143 (MQ=255)                                                                                                                                                                     atatatCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCa                                                                                                                                                                             >  1:386851/1‑139 (MQ=255)                                                                                                                                                                       atatCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAg                                                                                                                                                                       <  1:516137/143‑1 (MQ=255)                                                                                                                                                                         atCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCg                                                                                                                                                                         <  2:386846/139‑1 (MQ=255)                                                                                                                                                                               gcAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGTCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCCTTACGCAgg                                                                                                                                                                           >  2:493713/1‑131 (MQ=255)                                                                                                                                                                                   gAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGTCTGCTGACT‑AAGTCCGAGTTGTTCGCGAGCATTACGCCGGCGAg                                                                                                                                                                       <  1:493721/131‑1 (MQ=255)                                                                                                                                                                                        ttGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGa                                                                             >  1:552497/1‑216 (MQ=255)                                                                                                                                                                                            aGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACca                                                                         <  2:552486/216‑1 (MQ=255)                                                                                                                                                                                                         ggTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTACTCGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGTTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACTATATAAAACCCTCtgtgttc  >  1:27314/1‑272 (MQ=255)                                                                                                                                                                                                                  tcttctTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGttt                                                                                                     >  1:440730/1‑166 (MQ=255)                                                                                                                                                                                                                      cttcAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAAt                                                                                                 <  2:440724/166‑1 (MQ=255)                                                                                                                                                                                                                            tGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCTGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGc                                                                                                               >  1:590189/1‑146 (MQ=255)                                                                                                                                                                                                                            tGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTaa                                                                                                                                                           >  2:293273/1‑102 (MQ=255)                                                                                                                                                                                                                                gcgTACCGTGGAAACCTTCAGGCAAAGTCGCTCTGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGt                                                                                                           <  2:590178/146‑1 (MQ=255)                                                                                                                                                                                                                                gcgTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGc                                                                                                                                                       <  1:293278/102‑1 (MQ=255)                                                                                                                                                                                                                                 cgTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTaa                                                                                                                                                           >  1:431956/1‑97 (MQ=255)                                                                                                                                                                                                                                     ccGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGc                                                                                                                                                       <  2:431950/97‑1 (MQ=255)                                                                                                                                                                                                                                        tGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTcgcg                                                                                                                                                                                        >  2:302851/1‑61 (MQ=255)                                                                                                                                                                                                                                            aaCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCa                                                                                                                                                                                    <  1:302856/61‑1 (MQ=255)                                                                                                                                                                                                                                             aCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTattcattc                                                                                                                    >  2:299585/1‑124 (MQ=255)                                                                                                                                                                                                                                                 tCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTg                                                                                                                <  1:299590/124‑1 (MQ=255)                                                                                                                                                                                                                                                  cAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGc                                                                                                               >  2:285053/1‑124 (MQ=255)                                                                                                                                                                                                                                                      cAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGt                                                                                                           <  1:285058/124‑1 (MQ=255)                                                                                                                                                                                                                                                                             ctgACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgcgtg                                                                                                                                    >  2:632998/1‑76 (MQ=255)                                                                                                                                                                                                                                                                             ctgACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTatt                                                                                                                         >  2:470746/1‑87 (MQ=255)                                                                                                                                                                                                                                                                                 cT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCtt                                                                                                                                <  1:633009/76‑1 (MQ=255)                                                                                                                                                                                                                                                                                 cT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTattcatt                                                                                                                     <  1:470754/87‑1 (MQ=255)                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                         CCAGCCAAGTGAAGGTCATCAGCCAGCCGTCGCGTTTTTTGCGGCGTTTACCGAGGGAAAAACTCTGCATCGGCGCAACTTTTGCAGCCCGAAGTGGAGCCTGCTTTTCCAGCCCTGGCAGCAGTTCTTCTTCTGGAATATGTACCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACT‑AAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACAGTATAAGACCGTCTGTGTGC  >  CP009273/2635846‑2636319

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: ‑