Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,592	C→A	P294P (CCG→CCT)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,592	0	C	A	81.8%	33.1 / ‑1.3	22	P294P (CCG→CCT)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (7/11);  total (9/13)
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

GCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCC  >  NC_000913/2446351‑2446819                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    gCAGGGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCGTCACATCGGCATTTGGCGCCCGTTGCCGTAACAGGTGA‑TCCAGTAAAACGAGCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAgc                                                                                                                                                                                                                            <  2:32766/251‑1 (MQ=255)       ccAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTTAGCGccct                                                                                                                                                                                                                       >  2:81815/1‑248 (MQ=255)                       ggTTTTATTCATTTTTTCCCAGCGTGGAATATCAGTCTTCAAATCTGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTTCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTtgct                                                                                                                                                                                                      <  2:242963/251‑1 (MQ=255)                          tttATTCATTTTTTACCATCGTGGAAGATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCTTAACAGGTGA‑TCAAGTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACTGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTtgctgc                                                                                                                                                                                                    <  2:81699/250‑1 (MQ=255)                               tCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGATGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCgct                                                                                                                                                                                               <  1:158355/250‑1 (MQ=255)                               tCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCAGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATG‑TACGACCCGGCACGGTCATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCgctgg                                                                                                                                                                                             >  2:87618/1‑250 (MQ=255)                                          cAGCGTGGAATATCATTCTCCACATCGGCATTTTGCGCCCTTTGCAGTAACGGGTGA‑TCCATTAAAACTTCCGCCAGCACCGCTTCTGCGATCGGCACTTCGCGGATCTCGACACAGTGATAGTGACGGCCTTTGGTGATCATCGCAACTTCTT‑CGCCAAAGCTGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccgcc                                                                                                                                                                                   <  2:274776/251‑1 (MQ=255)                                                  aaTATCAGTCTTCCCATCGGCATTTTGCGCCCGTTGCCGTAAAAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGACAAAGCGGTTAATGGTACGACCCGACACGGTAATGCTGGAGGTAGGTTTCAGCGACATATGGGCTATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgc                                                                                                                                                                           >  1:52449/1‑251 (MQ=255)                                                       cAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGccgccc                                                                                                                                                                      >  1:30946/1‑251 (MQ=255)                                                                 aTCGGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGc                                                                                                                                                            >  1:74612/1‑251 (MQ=255)                                                                                  gTTGCCGTAACAG‑TGAACCCATGAAAACGAGCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTg                                                                                                                                           <  1:194332/251‑1 (MQ=255)                                                                                                                           cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCGGGCTGCCGCGCAGCGccacc                                                                                                    >  1:88561/1‑250 (MQ=255)                                                                                                                                     aTCGGCACTGCGCGGATCCCGCCACAGGGATCGTGACGGCCTTTGGTGATCATCTCATCTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGcc                                                                                         <  1:52292/251‑1 (MQ=255)                                                                                                                                       cgccccgtgccgaaTCCCGCCACAGGGACCGGGACGGCCTTTGGTGATCATCTCAACTTCTT‑CTCCAAAGGTGTTAATGGTCCGCCCCGGCACGGTAAGGCTGGAGGGAGTTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGt                                                                                       <  2:274529/238‑1 (MQ=255)                                                                                                                                            gcgcgtggaccccGACACAGGGATCGTGACGTCCTTTGGGGATCATCTCAATTTCTT‑CGCCAAAGCGGTTAAGGGTACGCCCCGGCGCGGGAATGCTGGGGGGAGTTTACAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCa                                                                                  <  2:117890/241‑1 (MQ=255)                                                                                                                                                gcggacccatacacAGGGATGGTGCCGGCCTTTGGGTATCATCTCAACTTCTC‑CGCAAATGCGGTAAATGGGACGCCCCGGCACGGTATTGCTGGGGGCCGGTTCCAGCCCCATATGGGCAATGATTTGCTGCCCGTTGCTGAAACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAAttt                                                                              <  2:30946/241‑1 (MQ=255)                                                                                                                                                          acccaGGGTTCGTGACGGCTGTGGGTGATCATCTCA‑GTTCTTGCGCCAACGCGGTTAATGGTTCGACCCGGCACGGTAATTCTGGGGGTAGGTATCAGCGCCATATGGGCATTGTTTTTCTGCCCGCTGCCTATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCGTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCttt                                                                    <  2:159201/248‑1 (MQ=255)                                                                                                                                                                     gTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACACGGCACGGTAATGCTGGCGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGGTGCTCTGGAAACCGTCTTTGGTGATTTCATCACGTTTCTGGCTGCCGCGCAGCGCAACCACGCCAACGCCGTCGCCAATTTCCACGACTTTCACcgcgatgg                                                         >  2:155228/1‑247 (MQ=255)                                                                                                                                                                             cTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGCGAATGCCGCCAGACTGGCTGCTCTGGAAACCGTCTCTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCCCCCCGTCAAAGCCGTCGCCAATTTCAACGGCTTTCACCCCGTTTATgataaac                                                 >  2:177665/1‑245 (MQ=255)                                                                                                                                                                                                  tctccGCCAAAGCGGTTAATGGTACGCCCCGGCACGGTAATGCTGGAGGTAGGTTCCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAg                              <  1:115779/246‑1 (MQ=255)                                                                                                                                                                                                       cccAAAGCGGTTAATG‑TACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATTGGCGATGTCAGCATCCAg                       <  1:87618/250‑1 (MQ=255)                                                                                                                                                                                                                              ccGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGcc  >  1:44976/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    GCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGA‑TCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTT‑CGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCC  >  NC_000913/2446351‑2446819

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