Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,759	G→A	G193G (GGC→GGT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,759	0	G	A	93.1%	90.9 / ‑5.4	29	G193G (GGC→GGT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/4);  minor base C (0/2);  total (23/6)
Fisher's exact test for biased strand distribution p-value = 3.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCC  >  NC_000913/67514‑67984                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  aaCGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATc                                                                                                                                                                                                                              <  2:146176/251‑1 (MQ=255)        ggAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGAcc                                                                                                                                                                                                                       >  1:289940/1‑251 (MQ=255)           aGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGcc                                                                                                                                                                                                                    <  1:171681/251‑1 (MQ=255)                                 cTTCATCCCAAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCAGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTAATCATCGACCAGAGCGTTAACCTCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCCACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCa                                                                                                                                                                                               >  1:160380/1‑250 (MQ=255)                                                         cgcTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCACCCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCATAACGGCAGACTTTCAGATGACggg                                                                                                                                                                      >  1:117764/1‑251 (MQ=255)                                                                aGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAACAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTg                                                                                                                                                               >  1:290793/1‑251 (MQ=255)                                                                             gTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGcc                                                                                                                                                  >  2:92845/1‑251 (MQ=255)                                                                               ttttcaccGTGGACTTGTGTTCCAGGCGTCATGGTGTCGCAGCTTTCGTCCTCATCGACCAGCGCGTTAAAATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGTAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACTGCCACTTCACGCATGTTATCGCCAAAACGGACGACTTTCTGATGACTGTCTTCCTGTTTAGAAACCGCCTg                                                                                                                                                >  2:132376/7‑251 (MQ=255)                                                                                   cgCCGTGGTTTTGTTTTGCAGGCGTCATGGTGTAGCAGCTTGCGTACTCATCGACCAGCGCTTTAACATCCCCGTCGCTGATGGAGTCCACCCCCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGc                                                                                                                                            <  1:112773/250‑1 (MQ=255)                                                                                    gCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGACAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGc                                                                                                                                            >  1:285511/1‑250 (MQ=255)                                                                                      cGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATc                                                                                                                                         >  1:42455/1‑251 (MQ=255)                                                                                           aTTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAgg                                                                                                                                     >  1:111170/1‑250 (MQ=255)                                                                                                         cGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCAt                                                                                                                      >  1:12915/1‑251 (MQ=255)                                                                                                          gTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATg                                                                                                                     >  1:273348/1‑251 (MQ=255)                                                                                                               ggTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCtt                                                                                                                >  1:56705/1‑251 (MQ=255)                                                                                                                 tgtAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCTACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGt                                                                                                              >  2:218656/1‑251 (MQ=255)                                                                                                                              gTACTCATCGCCCGGCGCGTTAACACCGCCTCCGCTGAGGGAGTTCCCCCCCGGCACCAGATCGCCACCCCCCCAGGTTTTGAGGGAAAACCCAAACTTGTTCTGTGCGGCAACTTTACCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGt                                                                                                 <  1:231156/251‑1 (MQ=255)                                                                                                                                                      aTCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCt                                                                         >  1:319688/1‑251 (MQ=255)                                                                                                                                                           cGTCGCTGATGGAGTTAACCACCTGCACCAGATCGCCAACCGCCCAGGAATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGc                                                                    >  1:255883/1‑251 (MQ=255)                                                                                                                                                            gTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGc                                                                    >  1:243629/1‑250 (MQ=255)                                                                                                                                                               gCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAc                                                                >  1:245537/1‑251 (MQ=255)                                                                                                                                                                 tGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGa                                                                                                                                               >  2:131474/1‑170 (MQ=255)                                                                                                                                                                 tGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcgc                                                              >  1:292659/1‑251 (MQ=255)                                                                                                                                                                 tGATGGAGTCCCCCACCTGCACCAGACCGCCACCCCCCCAGGTATTGACGGAAAAACCAACCTGGATCTGGGCGCCAATTTTACCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGa                                                                                                                                               <  1:131474/170‑1 (MQ=255)                                                                                                                                                                  gATGGAGTTCACCACCTGCACCAGCTCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGCCTTTCAGATGACGGGTATCCTGTTTAGAACCCGCCTGACGACTCCAGAAGCCGATACGCTCATGGGCATGTTTATCCTGCCAGTGACCGGTAACGACGGCCTGGTGCTGAATCATAcgcc                                                             >  2:118011/1‑250 (MQ=255)                                                                                                                                                                                                   cGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCAtgtg                            >  1:113268/1‑251 (MQ=255)                                                                                                                                                                                                                         ccGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCCTGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGCCGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGCACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATaa       >  2:29449/1‑250 (MQ=255)                                                                                                                                                                                                                            accTTGATCTGTGCGACAACTTAACCACCACGGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTc   <  1:26193/249‑1 (MQ=255)                                                                                                                                                                                                                              cTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTcc  >  2:54253/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  AACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCC  >  NC_000913/67514‑67984

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