Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,611	G→A	Y146Y (TAC→TAT)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,611	0	G	A	95.8%	76.7 / NA	25	Y146Y (TAC→TAT)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/13);  minor base C (1/0);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTCCAGCTCATCAAAGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGC  >  NC_000913/69365‑69845                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            cGGGTCCAGCTCTTCACAGCACCTGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGCGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGCGTGGTCCCGGAAAGCAGCGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGCTTTTTGCCCAGAACCATTCGCTGGAATACATACCGCCAATCTAGc                                                                                                                                                                                                                                         >  2:24817/1‑251 (MQ=255)                                       cgccCCAGCTTTCGTGCCACAGAGTTTTATGCCCGGCGCTGCAACGTCCGCGCCGAATATCCTGGGGGCGGGTGGTACCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACgg                                                                                                                                                                                                   <  1:219271/251‑1 (MQ=255)                                                                   gcgCCCGGCGCGGCAACGCCCGCAACGAATACCCGGCGGGGGGGGGTTACGGGAAAGCAGAGCGGCCCCCCCGTCGCCCACCGCAATCCACGATGCGGCAGATGCCGCCAGGGCGGTGTCCTGGCGGGTCACAGGCGGTATTTTTGCCCAGACCCATTCGCGGGAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAg                                                                                                                                                                      <  2:257950/249‑1 (MQ=255)                                                                            cgtgaaaCGGCCGCGACGAATATCCTGCGGGCGGGTGGCACCGGAAAGCAGAGCTGCCACCCAGTCGAACAGCTCAACCCCCGAGGCGGCAGATTGCGCCACGGCGCTGTCCTGGGGAGTCACATGCAGGATTTTTGCCCAGACCCATCCGCGGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGTTCtt                                                                                                                                                             <  2:65206/246‑1 (MQ=255)                                                                                cAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTcc                                                                                                                                                          >  1:211181/1‑250 (MQ=255)                                                                                  aCGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCcac                                                                                                                                                        >  2:212759/1‑250 (MQ=255)                                                                                    gTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAAt                                                                                                                                                     >  2:65023/1‑251 (MQ=255)                                                                                         ggacgaATATCCTGCGGGCGGTGGGTACCGGAAAGCAGAGCGGGCACCAATTCGCACAGACAAATCAACGTGGCGGAAGATTGCCCCACGGCGCTGTCCTGGCGGGTAACATGAAGGATTTTTGCCCAGAACCATTCGCTGGAATAAACACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGaa                                                                                                                                                <  1:233659/250‑1 (MQ=255)                                                                                         cgacgaATAGCCGGGGGGGGGGTGGGACCGAAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGaa                                                                                                                                                <  2:95784/251‑1 (MQ=255)                                                                                                   ccTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTc                                                                                                                                      >  1:84797/1‑251 (MQ=255)                                                                                                                     cGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCCCGAGGCGGCAGATTGCGCCACGGCGCTGCCCTGGGGATTCACAGGCAGGATTTTTGCCCAGACCCATGCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTGCGGCAAACTcc                                                                                                                    <  1:202000/251‑1 (MQ=255)                                                                                                                          aGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATGCCGCCCCGGCGCTTTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGACCCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                               <  1:153864/251‑1 (MQ=255)                                                                                                                                                            ccgaggCGGCAGATTGCCCCACGGCGCTGTCCGGGCGATTCACAGGCAGGATTTTTGCCCAAACCCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGt                                                                             <  1:10732/246‑1 (MQ=255)                                                                                                                                                                tGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCgg                                                                         >  2:142599/1‑251 (MQ=255)                                                                                                                                                                                  ggCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGt                                                       >  1:38755/1‑251 (MQ=255)                                                                                                                                                                                       tGTCCTGGCGAGCCAAATGATGGTTTTTTGCCAAGAACCTTACCCTGAAAAAAACTCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAAc                                                   <  2:18064/250‑1 (MQ=255)                                                                                                                                                                                               cGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                          >  1:92513/1‑251 (MQ=255)                                                                                                                                                                                                       aTGCAGGATTTTTGCCCAGAACCATTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAg                                  <  1:65023/251‑1 (MQ=255)                                                                                                                                                                                                                 tttGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGAGGGAGTAGTCAACATTGCAAGGCGCGTGGCACAACCGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCCCAATACGAACATCGCGTTCGGGCTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACACCAGCTGCGCGCTg                        >  1:125842/1‑251 (MQ=255)                                                                                                                                                                                                                    gccaagaacatttCGCTGAAAGAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCTAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGtt                      <  1:179603/239‑1 (MQ=255)                                                                                                                                                                                                                               ttCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTAAGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGct          >  1:136312/1‑251 (MQ=255)                                                                                                                                                                                                                                    gggAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGGCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCa     <  1:10892/250‑1 (MQ=255)                                                                                                                                                                                                                                     gaaataaataCCGCCAATTTAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCGCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGGGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTTTTCGACGCTAAGCTCTGCa     <  1:70613/248‑1 (MQ=255)                                                                                                                                                                                                                                     gaaaaaaatcccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCa     <  2:84797/237‑1 (MQ=255)                                                                                                                                                                                                                                       aataaataACGCCAATATAGCGGGAGCGGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGACTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGTGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGc  >  2:102407/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            CGGGTCCAGCTCATCAAAGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGC  >  NC_000913/69365‑69845

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