Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,849	G→A	G400G (GGC→GGT)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,849	0	G	A	96.9%	104.4 / NA	32	G400G (GGC→GGT)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/16);  total (15/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.92e-01

CGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGG  >  NC_000913/68612‑69073                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 cGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCCGCGGAGCGTCGGTAGAGAGTATAAGATAGGCAATTACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGATAAACCAgg                                                                                                                                                                                                                     >  2:179603/1‑250 (MQ=255)     cGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAg                                                                                                                                                                                                                 <  2:116465/251‑1 (MQ=255)      gATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGACGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGc                                                                                                                                                                                                                >  1:63708/1‑251 (MQ=255)          cAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGcacc                                                                                                                                                                                                             >  1:246115/1‑250 (MQ=255)                gTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCa                                                                                                                                                                                                      <  2:119758/251‑1 (MQ=255)                     agcagGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTg                                                                                                                                                                                                  >  1:101773/1‑250 (MQ=255)                                                            cgccCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCAATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCGTCGGTCAGCGcc                                                                                                                                                          <  2:19348/251‑1 (MQ=255)                                                                     ccATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCAGGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAgt                                                                                                                                                 <  2:104083/251‑1 (MQ=255)                                                                            gTTATTGACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTct                                                                                                                                          >  2:118998/1‑251 (MQ=255)                                                                                   aCGGCGATCCCCTGATCGGTAAAGCACTCCAGGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCtggctggc                                                                                                                                   <  2:251584/251‑1 (MQ=255)                                                                                        gATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCCGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTgg                                                                                                                              >  1:161859/1‑250 (MQ=255)                                                                                            cccTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGAtttt                                                                                                                          >  1:64535/1‑250 (MQ=255)                                                                                                   cGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGcttt                                                                                                                    <  1:204381/250‑1 (MQ=255)                                                                                                   cGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTc                                                                                                                   <  1:54707/251‑1 (MQ=255)                                                                                                       aaaGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCCGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTTCGCTTTCAGtt                                                                                                               >  2:157619/1‑251 (MQ=255)                                                                                                                                     ggCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCg                                                                                 >  1:200062/1‑251 (MQ=255)                                                                                                                                      gCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCgg                                                                                <  2:101773/251‑1 (MQ=255)                                                                                                                                       cGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGc                                                                               >  2:81336/1‑251 (MQ=255)                                                                                                                                              agcaATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGctga                                                                        <  2:37911/251‑4 (MQ=255)                                                                                                                                                             gAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAAcc                                                         <  2:41746/251‑1 (MQ=255)                                                                                                                                                                 agcagcGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGc                                                     <  2:234313/251‑1 (MQ=255)                                                                                                                                                                      cGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGa                                                <  1:60662/251‑1 (MQ=255)                                                                                                                                                                        gAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGat                                               >  1:197638/1‑250 (MQ=255)                                                                                                                                                                                  aGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACg                                    <  1:134874/251‑1 (MQ=255)                                                                                                                                                                                     gAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCg                                 >  1:65206/1‑251 (MQ=255)                                                                                                                                                                                       ggTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCGGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGAt                               <  1:203892/251‑1 (MQ=255)                                                                                                                                                                                            aGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGcc                          <  1:151060/251‑1 (MQ=255)                                                                                                                                                                                             gATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCt                         <  2:197638/251‑1 (MQ=255)                                                                                                                                                                                                    aaTCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTcc                   <  1:264768/250‑1 (MQ=255)                                                                                                                                                                                                     aTCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTcc                   >  1:9127/1‑249 (MQ=255)                                                                                                                                                                                                                  ggCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATcc     >  1:70768/1‑250 (MQ=255)                                                                                                                                                                                                                    cGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGGGACGGATTTTTGGCCCATGCCTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATCTCACAAAACGCCGATTGGCCTGCTTCCCGAACGATAAAGCCAgg  >  2:202000/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 CGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGG  >  NC_000913/68612‑69073

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