Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,104	A→G	G315G (GGT→GGC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,104	0	A	G	92.5%	117.9 / ‑0.8	40	G315G (GGT→GGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (22/15);  total (24/16)
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

GATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGC  >  NC_000913/68968‑69237                                                                                                                                         |                                                                                                                                      gATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCt                                                                                                                                    >  2:296991/1‑140 (MQ=255)      tGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAAc                                                                                                                               >  2:539594/1‑140 (MQ=255)        ggcggcAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTg                                                                                                                             >  2:110237/1‑140 (MQ=255)                              cAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGttt                                                                                                       >  2:1025191/1‑140 (MQ=255)                                   gAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCgg                                                                                                  <  1:296991/140‑1 (MQ=255)                                              cAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATg                                                                                       <  2:46856/140‑1 (MQ=255)                                                 aCCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAAACAGAATGTCg                                                                                    <  2:582941/140‑1 (MQ=255)                                                       cGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTg                                                                              >  2:515683/1‑140 (MQ=255)                                                       cGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTg                                                                              >  1:431884/1‑140 (MQ=255)                                                        gTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTgg                                                                             <  1:539594/140‑1 (MQ=255)                                                           gATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAg                                                                          >  1:6051/1‑140 (MQ=255)                                                           gATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAg                                                                          >  2:271893/1‑140 (MQ=255)                                                                cACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCg                                                                     >  1:841181/1‑140 (MQ=255)                                                                 aCCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGa                                                                    <  1:110237/140‑1 (MQ=255)                                                                 aCCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGa                                                                    <  2:6051/140‑1 (MQ=255)                                                                 aCCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGa                                                                    <  1:271893/140‑1 (MQ=255)                                                                   cAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATa                                                                  >  1:828665/1‑140 (MQ=255)                                                                     aaCTCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAAc                                                                >  2:170312/1‑140 (MQ=255)                                                                      aCGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACt                                                               >  2:839441/1‑140 (MQ=255)                                                                      aCGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACt                                                               >  2:386847/1‑140 (MQ=255)                                                                             ttGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCa                                                        <  1:515683/140‑1 (MQ=255)                                                                                  cTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCg                                                   <  1:170312/140‑1 (MQ=255)                                                                                    gCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGtt                                                 >  2:811673/1‑140 (MQ=255)                                                                                     cTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTa                                                >  1:526744/1‑140 (MQ=255)                                                                                               gctAAATCCAGCCACCACGCTGCCATCAACCTGACCGCAACTACCTTTAACTTCCCGCTCGCCAACCCTCTGTTTGTTGGCAATCAGAATGTCGCAGGTCGAAGTACCGATAACTTTTACCAGCGCGTTAGGCTGtgcgc                                      >  2:876547/3‑140 (MQ=255)                                                                                               gATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCa                                                        <  1:469093/122‑1 (MQ=255)                                                                                               gATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCa                                                        >  2:469093/1‑122 (MQ=255)                                                                                                aTAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcc                                     >  1:938902/1‑140 (MQ=255)                                                                                                aTAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcc                                     >  1:541874/1‑140 (MQ=255)                                                                                                aTAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcc                                     <  1:839441/140‑1 (MQ=255)                                                                                                       cAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCa                              <  2:841181/140‑1 (MQ=255)                                                                                                         ggCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAAc                            <  2:938902/140‑1 (MQ=255)                                                                                                          gCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACt                           >  2:395666/1‑140 (MQ=255)                                                                                                              cacGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTgcgc                       <  2:526744/140‑1 (MQ=255)                                                                                                              cacGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTgcgc                       <  1:386847/140‑1 (MQ=255)                                                                                                                               gACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAAc      >  1:145945/1‑140 (MQ=255)                                                                                                                                 ccGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAAcgc    >  1:152917/1‑140 (MQ=255)                                                                                                                                  cGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAAcgcg   >  1:420002/1‑140 (MQ=255)                                                                                                                                      aaTGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACgcgc  <  2:569800/137‑1 (MQ=255)                                                                                                                                      aaTGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACgcgc  >  1:569800/1‑137 (MQ=255)                                                                                                                                         |                                                                                                                                      GATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGC  >  NC_000913/68968‑69237

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