Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,596	(T)6→7	coding (742/1503 nt)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,590	1	.	T	100.0%	129.2 / NA	37	R250R (CGA→AGA)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (17/20);  total (17/20)

TGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑TTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGAC  >  NC_000913/67457‑67722                                                                                                                                       |                                                                                                                                     tGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑tttttt                                                                                                                                <  1:730609/140‑1 (MQ=255)   tACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑tttttttg                                                                                                                              <  2:1596/140‑3 (MQ=255)   tACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑tttttttg                                                                                                                              >  1:78817/1‑138 (MQ=255)     cGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑tttttttgcc                                                                                                                            <  1:166874/140‑5 (MQ=255)      ggCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑tttttttgccg                                                                                                                           <  1:1452/140‑6 (MQ=255)        ccAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTg                                                                                                                          >  1:548131/1‑140 (MQ=255)           gACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGAt                                                                                                                       >  1:6444/1‑140 (MQ=255)            aCCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGAtt                                                                                                                      >  1:762135/1‑140 (MQ=255)                ggAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTtgtg                                                                                                                  <  2:762417/140‑1 (MQ=255)                  aaGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTtgtgtg                                                                                                                <  2:772322/140‑1 (MQ=255)                      tGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAg                                                                                                            <  2:776423/140‑1 (MQ=255)                       gTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAgg                                                                                                           <  2:692419/140‑1 (MQ=255)                            aGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCa                                                                                                      <  2:548131/140‑1 (MQ=255)                             gACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCAt                                                                                                     >  1:286024/1‑140 (MQ=255)                              aCCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATg                                                                                                    <  2:266178/140‑1 (MQ=255)                                  tGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGc                                                                                                              <  2:141431/126‑1 (MQ=255)                                  tGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGc                                                                                                              >  1:141431/1‑126 (MQ=255)                                  tGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGgtgt                                                                                                >  1:94238/1‑140 (MQ=255)                                  tGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGgtgt                                                                                                >  2:94603/1‑140 (MQ=255)                                  tGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGgtgt                                                                                                >  2:908528/1‑140 (MQ=255)                                  tGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGgtgt                                                                                                >  2:901244/1‑140 (MQ=255)                                     aaaTCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTagc                                                                                             <  2:909193/140‑1 (MQ=255)                                       aTCTTCAAAGGTGGTGGTGAACGCGTGGAAGCAACCTTGTTCCAGTAAACGCTTCAGCCCCGGCTCAATACGCTCCGCTTCCAGCACGTACTGTCGTTTTTTTGCCGTGGTATTGTGTGGCAGACGTGAGGGAGTagcag                                                                                           >  2:519680/1‑140 (MQ=255)                                              aaGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGt                                                                                    <  2:12480/140‑1 (MQ=255)                                                    gtggtgAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCAt                                                                              >  2:840202/1‑140 (MQ=255)                                                         gAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGAcc                                                                         <  1:728050/140‑1 (MQ=255)                                                          aaCGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCa                                                                        <  1:87227/140‑1 (MQ=255)                                                           aCGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAg                                                                       <  1:94603/140‑1 (MQ=255)                                                                    aGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAAc                                                              <  2:286024/140‑1 (MQ=255)                                                                         ccTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGc                                                         >  2:734144/1‑140 (MQ=255)                                                                                   ggAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGAt                                               >  2:709887/1‑140 (MQ=255)                                                                                           ttCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTca                                       <  1:840202/140‑1 (MQ=255)                                                                                           ttCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTca                                       <  2:78817/140‑1 (MQ=255)                                                                                               tCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAACTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTcaccac                                   <  1:901244/140‑1 (MQ=255)                                                                                                   cAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGc                               <  1:73997/140‑1 (MQ=255)                                                                                                     gCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTca                                       <  2:148346/130‑1 (MQ=255)                                                                                                     gCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTca                                       >  1:148346/1‑130 (MQ=255)                                                                                                     gCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCAc                             >  1:565182/1‑140 (MQ=255)                                                                                                                gccgcTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCaa                  <  1:908528/140‑1 (MQ=255)                                                                                                                 ccgcTTCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAAc                 <  2:94238/140‑1 (MQ=255)                                                                                                                     ttCCAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGcc             >  1:610403/1‑140 (MQ=255)                                                                                                                       ccAGCACGTTCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCa           >  2:191118/1‑140 (MQ=255)                                                                                                                                tCTGTCGTTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGAc  >  2:19637/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                     TGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCG‑TTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGAC  >  NC_000913/67457‑67722

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