Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,400	G→A	L217L (CTG→TTG)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,400	0	G	A	97.0%	98.9 / NA	33	L217L (CTG→TTG)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/10);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAAGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCA  >  NC_000913/69154‑69626                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   cGCAGGTGGAAGTACAGATAACTCTTACGGGTGTGTTAGGCTGTGTTCCTGCGCCAACTGCGCCCATAGGGGAGTCAATCGCGCCGCCGGAAATCACCACAGTTTCAGACCGGCCGAGACGCTACGCCCATTCCTGGCATAAGGTGCCCACCGGAATATCGGCGGTCCAGGGGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCATCTCAGCAAAGAAAC‑TGGATGGCGGCAAgcct                                                                                                                                                                                                                                <  2:204988/251‑2 (MQ=255)       tGGAAGCACCGAAAAATTTTAAAAGCGCGATAGGCTGGGCGCCTGCGCAAACTGCGAACATATGGAAGTCAAACGCGCAGCCGGAAATCACCACGCTGTCAGGAAGGACGAGACGAGGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGAGGCAGTGAAAAGAGGGGAAGGGAAATGGCGATTGAGGATAGGGGCCAGCTAATCAAAGAAAC‑AGGCGGGCGGCAAGCCGcccc                                                                                                                                                                                                                            >  2:220687/1‑249 (MQ=255)                cGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGc                                                                                                                                                                                                                 >  1:78771/1‑251 (MQ=255)                 gATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGcc                                                                                                                                                                                                                >  1:286239/1‑251 (MQ=255)                                    ttAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGgc                                                                                                                                                                                             >  2:264497/1‑251 (MQ=255)                                     tAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGgcg                                                                                                                                                                                            >  1:24294/1‑251 (MQ=255)                                           gtgCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGATATTTATGCCCGGCGCTGCaa                                                                                                                                                                                      <  1:271667/251‑1 (MQ=255)                                                   gcgccAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCcgcg                                                                                                                                                                              >  1:194320/1‑251 (MQ=255)                                                                 cATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTgcg                                                                                                                                                                <  1:35156/251‑1 (MQ=255)                                                                     tGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAAGATCCTgcgggcg                                                                                                                                                            >  2:117766/1‑251 (MQ=255)                                                                                   gccgccGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGaaa                                                                                                                                              <  1:243390/251‑1 (MQ=255)                                                                                    ccgccgGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGaa                                                                                                                                               >  1:157950/1‑249 (MQ=255)                                                                                                               aggcCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGGGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATCACCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTc                                                                                                                  >  2:56350/1‑251 (MQ=255)                                                                                                                    gagaCGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAAGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATcc                                                                                                             >  2:156250/1‑251 (MQ=255)                                                                                                                         gCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGAt                                                                                                        >  1:195946/1‑251 (MQ=255)                                                                                                                                 cATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGc                                                                                                   >  2:56335/1‑248 (MQ=255)                                                                                                                                 cATTACGGGCATAAGGTGCCCACCGGAATATCGGCAGTGCAGGAGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGGTTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCATGCGGGCGGGGGGTACCGGAAAGCAGAGCTGGCACCCAGGCGCACAgatacatccctggtgcggccg                                                                                                 >  2:11957/1‑230 (MQ=255)                                                                                                                                     ccGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGCTGCGGCAGATTgc                                                                                            >  2:113877/1‑251 (MQ=255)                                                                                                                                        ggCATAAGGTGCCCACCGGAATATCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGGTTGCGcc                                                                                         >  2:267197/1‑251 (MQ=255)                                                                                                                                                              atCGGCAGTCCAGGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTc                                                                    >  1:265946/1‑250 (MQ=255)                                                                                                                                                                       ccaggtgtCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑GGGCTGGCGGCGAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTTGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACTGCGCTGGCCTGTCGAGTCACATGCAgg                                                           <  2:39172/246‑1 (MQ=255)                                                                                                                                                                                 gTGAACAGCGGGGAAGGCAAATGGCGAGTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCGGGCGAGTCACATGCAGGATTTTTGCCCa                                                <  1:102263/251‑1 (MQ=255)                                                                                                                                                                                   gAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGccc                                                 >  1:278539/1‑248 (MQ=255)                                                                                                                                                                                   gAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAg                                               <  1:224235/250‑1 (MQ=255)                                                                                                                                                                                         ggggtaaGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCAtt                                        <  1:198875/246‑1 (MQ=255)                                                                                                                                                                                          ggggAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTc                                       >  1:112163/1‑251 (MQ=255)                                                                                                                                                                                            ggAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TACCTGGCGGCAAGCCGCCCAAGCTTTCGTGCCAAAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGCATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGAGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGc                                     >  1:188903/1‑251 (MQ=255)                                                                                                                                                                                                ggCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGa                                 <  1:222163/251‑1 (MQ=255)                                                                                                                                                                                                 gCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGCTTTTTGCCCAGAACCATTCGCTGGa                                 >  2:150910/1‑250 (MQ=255)                                                                                                                                                                                                       ggCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGaataaata                          <  1:131143/251‑1 (MQ=255)                                                                                                                                                                                                                ggaccGGGTCCAGCTCA‑AAAAGAAACCTGGCGGGGGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCGGGCGCAGCAACGCCCGAGAAGAATATCCTGCGGGCGGGTGGTACCGAAAAAAAGTGCGGCCCCCCAGTCGCACCGCTAAACCCTCGATGCGGCAGATTGCGCCACGGCGCTGTCCGGGCGATTCACATGCAGTATTTTGGCCCAGCCCCTTTCGCGGGAATAAATGCCGCCAAta                 <  2:218858/247‑2 (MQ=255)                                                                                                                                                                                                                                atcaAAGAAAC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCa  >  2:2122/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaaC‑TGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTg                                                    >  1:22736/1‑193 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   CGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCAGGCAGGCCGAGACGCTGCGCCCATTCCGGGCATAAGGTGCCCACCGGAATATCGGCAGTCCAAGTGTCAGTGAACAGCGGGGAAGGCAAATGGCGATTGAGGATCGGGTCCAGCTCATCAAAGAAAC‑TGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCA  >  NC_000913/69154‑69626

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