Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,264	G→A	I358I (ATC→ATT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,264	0	G	A	92.9%	80.3 / ‑1.2	28	I358I (ATC→ATT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (12/14);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGAC  >  NC_000913/67020‑67487                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                cGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTTAACTGCGAGGATCGGTTTCTCGTCTACGGCAATCGAc                                                                                                                                                                                                                           <  2:133775/251‑1 (MQ=255)    aaGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCAACCGGCAGTTTCGGAAGGGAGTGCGGTGTTTTCCCCGTGTCGATACAGTTAACAAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTCCCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACggg                                                                                                                                                                                                                        >  1:175104/1‑251 (MQ=255)      gCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCa                                                                                                                                                                                                                      <  2:131893/251‑1 (MQ=255)                           gcTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGcc                                                                                                                                                                                                 <  1:191711/251‑1 (MQ=255)                                  acaGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCAcc                                                                                                                                                                                          <  2:247413/251‑1 (MQ=255)                                                      aGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGa                                                                                                                                                                       >  1:190460/1‑250 (MQ=255)                                                         ttCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGGAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCATGTCATTACCTTTCTCGAAGt                                                                                                                                                                    <  2:29188/250‑1 (MQ=255)                                                                 ggAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCCTTACCTTTCTCGAAGTGATAGgtgt                                                                                                                                                           >  1:238944/1‑251 (MQ=255)                                                                               tcgccGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATaaa                                                                                                                                              <  1:108188/247‑1 (MQ=255)                                                                                cACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGGCCTCCATAAagg                                                                                                                                            <  1:196243/251‑1 (MQ=255)                                                                                cACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGGATTGAAGATCAGTCGGGCAGGATCATCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAagg                                                                                                                                            >  2:37119/1‑251 (MQ=255)                                                                                        cGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTgccg                                                                                                                                    <  2:190460/251‑1 (MQ=255)                                                                                          atacAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTgccgcc                                                                                                                                  <  1:68093/248‑1 (MQ=255)                                                                                                                   gATCGCCGAGATCAAGCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTc                                                                                                         <  1:176331/251‑1 (MQ=255)                                                                                                                             atcaatcaGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGa                                                                                                >  1:191999/1‑250 (MQ=255)                                                                                                                               caatcagGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGCAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATCTGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGGGAAGc                                                                                             >  2:218903/1‑251 (MQ=255)                                                                                                                                                   gctggaacGGTCTGGGTATTGAAGATACGTCCGGCTGGCTCGTCCTTACCAACAATACCGTGCTGCTCAAAGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCCCCCGGTCATTACCTTTCTCGAAGTGATAGGCGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGATGCCGGGCGGCGGTTTTCCAGTc                                                                         >  2:30397/8‑251 (MQ=255)                                                                                                                                                        accGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCtt                                                                    <  1:206233/250‑1 (MQ=255)                                                                                                                                                         ccGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTc                                                                   >  2:218459/1‑251 (MQ=255)                                                                                                                                                                 ggTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTTCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCaa                                                           <  1:160457/251‑1 (MQ=255)                                                                                                                                                                 ggTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCa                                                            >  1:190014/1‑250 (MQ=255)                                                                                                                                                                 ggTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTAACTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGGGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCAAGTCGCATTCGCCCGCa                                                            >  1:190023/1‑250 (MQ=255)                                                                                                                                                                      tGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGcc                                                       <  1:229514/250‑1 (MQ=255)                                                                                                                                                                                     cAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCACAGCCGTAACCCTGCTGCATc                                       >  2:35738/1‑251 (MQ=255)                                                                                                                                                                                     cAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATc                                       >  1:18233/1‑251 (MQ=255)                                                                                                                                                                                                   ccaccaATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGc                         <  1:244163/251‑1 (MQ=255)                                                                                                                                                                                                                         cGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTACAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGa   <  2:183522/251‑1 (MQ=255)                                                                                                                                                                                                                          gTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGAc  >  2:119163/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                CGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGAC  >  NC_000913/67020‑67487

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