Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,405	A→G	T311T (ACT→ACC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,405	0	A	G	96.7%	93.7 / ‑3.9	30	coding (933/1503 nt)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (12/17);  minor base . (0/1);  total (12/18)
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.91e-01

CGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAC  >  NC_000913/67158‑67650                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       cGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGtt                                                                                                                                                                                                                                                     <  2:45018/251‑1 (MQ=255)              accGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCtt                                                                                                                                                                                                                                       <  1:72544/250‑1 (MQ=255)                                     gtcGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCctgc                                                                                                                                                                                                                >  2:112773/3‑251 (MQ=255)                                            aGGATCGTCCTTACCTCCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTGGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCa                                                                                                                                                                                                         <  2:286185/251‑1 (MQ=255)                                               aTCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGAc                                                                                                                                                                                                      <  1:42078/251‑1 (MQ=255)                                                       tACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACg                                                                                                                                                                                              >  1:318244/1‑251 (MQ=255)                                                        accaccAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGAAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACgg                                                                                                                                                                                             >  1:138349/1‑251 (MQ=255)                                                             caATACCGAGCTGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCCGTACGGcccg                                                                                                                                                                                         >  2:195320/1‑248 (MQ=255)                                                             caATACCGAGATGCTGAACGTCTAGGATCGTTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCTCCTTCATGATGCGAAGCAGGGCTGC‑GTTTTGCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGa                                                                                                                                                                                        <  1:141467/251‑1 (MQ=255)                                                                 aCCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGAcctg                                                                                                                                                                                    <  2:121731/251‑3 (MQ=255)                                                                       aTGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCt                                                                                                                                                                              <  2:316291/251‑1 (MQ=255)                                                                                                       cggcaATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAggtggtg                                                                                                                                              <  2:99167/251‑1 (MQ=255)                                                                                                            aTCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGaa                                                                                                                                         <  1:220313/251‑1 (MQ=255)                                                                                                                                cATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCTTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGtt                                                                                                                     <  1:116385/251‑1 (MQ=255)                                                                                                                                        aGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTCCGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGaa                                                                                                              >  2:26193/1‑250 (MQ=255)                                                                                                                                                                   cGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCTCCTTCATGATGCGAAGCAGGGCGGCGGTGTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTGCAAAGGTGGTGGGGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATAcgc                                                                                     <  2:320392/248‑1 (MQ=255)                                                                                                                                                                   cGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGcgcc                                                                                  >  2:153783/1‑251 (MQ=255)                                                                                                                                                                    gAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCCCCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTTGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGcgccg                                                                                 <  1:199984/251‑1 (MQ=255)                                                                                                                                                                    gAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCGTCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGcgccg                                                                                 >  1:88556/1‑251 (MQ=255)                                                                                                                                                                                   gTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCt                                                                  <  1:303278/251‑1 (MQ=255)                                                                                                                                                                                             aaGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGtttttc                                                        >  2:231156/1‑250 (MQ=255)                                                                                                                                                                                             aaGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGCCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGtttttc                                                        >  2:231140/1‑250 (MQ=255)                                                                                                                                                                                                 aggTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCg                                                    >  1:15908/1‑251 (MQ=255)                                                                                                                                                                                                            tGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTtgtgt                                         >  2:12473/1‑251 (MQ=255)                                                                                                                                                                                                                 aCCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAg                                    <  2:48745/251‑1 (MQ=255)                                                                                                                                                                                                                                   cATGATGCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTagcag                   <  1:56018/250‑1 (MQ=255)                                                                                                                                                                                                                                       atgCGAAGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTc              <  1:39374/251‑1 (MQ=255)                                                                                                                                                                                                                                             aGCAGTGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCTCCTTGTGCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCGTTCGTGCTc        <  2:272495/251‑1 (MQ=255)                                                                                                                                                                                                                                             aGCAGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCATACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCTAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCTCCGTGGAGTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc        <  1:317909/251‑1 (MQ=255)                                                                                                                                                                                                                                                aGGGCGGCGGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTCCGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCCCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGCTTTGTGTTGCAGGCCTCATGGTGTAGCAGTCTTCGTACTCATc     >  2:66398/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gcggcgGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAc  <  2:15908/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTT‑CCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAC  >  NC_000913/67158‑67650

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