Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,471	A→T	P289P (CCT→CCA)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,471	0	A	T	100.0%	91.5 / NA	28	P289P (CCT→CCA)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (10/18);  total (10/18)

ATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGT  >  NC_000913/67228‑67716                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      aTGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCt                                                                                                                                                                                                                                                <  2:316291/251‑1 (MQ=255)                                 cggcaATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAggtggtg                                                                                                                                                                                                                <  2:99167/251‑1 (MQ=255)                                      aTCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGaa                                                                                                                                                                                                           <  1:220313/251‑1 (MQ=255)                                                          cATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCTTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGtt                                                                                                                                                                                       <  1:116385/251‑1 (MQ=255)                                                                  aGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTCCGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGaa                                                                                                                                                                                >  2:26193/1‑250 (MQ=255)                                                                                             cGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCTCCTTCATGATGCGAAGCAGGGCGGCGGTGTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTGCAAAGGTGGTGGGGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATAcgc                                                                                                                                                       <  2:320392/248‑1 (MQ=255)                                                                                             cGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGcgcc                                                                                                                                                    >  2:153783/1‑251 (MQ=255)                                                                                              gAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCCCCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTTGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGcgccg                                                                                                                                                   <  1:199984/251‑1 (MQ=255)                                                                                              gAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCGTCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGcgccg                                                                                                                                                   >  1:88556/1‑251 (MQ=255)                                                                                                             gTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCt                                                                                                                                    <  1:303278/251‑1 (MQ=255)                                                                                                                       aaGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGtttttc                                                                                                                          >  2:231156/1‑250 (MQ=255)                                                                                                                       aaGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGCCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGtttttc                                                                                                                          >  2:231140/1‑250 (MQ=255)                                                                                                                           aggTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCg                                                                                                                      >  1:15908/1‑251 (MQ=255)                                                                                                                                      tGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTtgtgt                                                                                                           >  2:12473/1‑251 (MQ=255)                                                                                                                                           aCCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAg                                                                                                      <  2:48745/251‑1 (MQ=255)                                                                                                                                                             cATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTagcag                                                                                     <  1:56018/250‑1 (MQ=255)                                                                                                                                                                 atgCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTc                                                                                <  1:39374/251‑1 (MQ=255)                                                                                                                                                                       aGCAGTGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCTCCTTGTGCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCGTTCGTGCTc                                                                          <  2:272495/251‑1 (MQ=255)                                                                                                                                                                       aGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCATACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCTAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCTCCGTGGAGTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc                                                                          <  1:317909/251‑1 (MQ=255)                                                                                                                                                                          aGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTCCGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCCCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGCTTTGTGTTGCAGGCCTCATGGTGTAGCAGTCTTCGTACTCATc                                                                       >  2:66398/1‑251 (MQ=255)                                                                                                                                                                             gcggcgGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAc                                                                    <  2:15908/251‑1 (MQ=255)                                                                                                                                                                                                        gCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGAt                                         >  2:65110/1‑251 (MQ=255)                                                                                                                                                                                                         cAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATg                                        <  2:31717/251‑1 (MQ=255)                                                                                                                                                                                                           aaGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGGTTCAGACCGTGCAAAGCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCTCCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGa                                      <  2:320629/251‑1 (MQ=255)                                                                                                                                                                                                                      cggcggcATCAGCCGCTGTACGGCAAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCCCGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCt                           <  1:171186/246‑1 (MQ=255)                                                                                                                                                                                                                                       gTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTGTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAAcc          <  2:233169/251‑1 (MQ=255)                                                                                                                                                                                                                                              cAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGCACGCGTGGAAGCCACCTTGTTCCAGGAACCGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCACCCGCCCAgg   >  1:203874/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGt  <  2:65910/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      ATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGT  >  NC_000913/67228‑67716

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