Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,983,288	C→T	*279* (TAG→TAA)	kduI ←	hexuronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,983,288	0	C	T	100.0%	145.5 / NA	46	*279* (TAG→TAA)	kduI	hexuronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/24);  total (22/24)

CAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAACTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGC  >  NC_000913/2983155‑2983419                                                                                                                                      |                                                                                                                                    caGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCATCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAgcg                                                                                                                                <  2:633983/139‑1 (MQ=255)              ccAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACg                                                                                                                   <  1:235394/139‑1 (MQ=255)               cAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACgg                                                                                                                  <  1:876285/139‑1 (MQ=255)                aGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTGACGGc                                                                                                                 >  2:449874/1‑139 (MQ=255)                    ccAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTGACGGCcaca                                                                                                             <  1:449874/139‑1 (MQ=255)                      aaCGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATg                                                                                                           <  1:24307/139‑1 (MQ=255)                       aCGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATgg                                                                                                          >  1:878777/1‑139 (MQ=255)                           cATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCa                                                                                                      <  2:878777/139‑1 (MQ=255)                                         gtccagtATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACg                                                                                                                   <  1:573083/112‑1 (MQ=255)                                         gtccagtATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACg                                                                                                                   >  2:573083/1‑112 (MQ=255)                                         gtcaagtATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACACTACTATATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGAcc                                                                                        >  1:550139/1‑139 (MQ=255)                                                         aGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGAc                                                                                         >  2:348810/1‑122 (MQ=255)                                                         aGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGAc                                                                                         <  1:348810/122‑1 (MQ=255)                                                           tGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATg                                                                      >  2:883509/1‑139 (MQ=255)                                                             acgacCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcac                                                                                                              >  2:560360/1‑97 (MQ=255)                                                             acgacCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCCCAAATCTTTAACGGCcac                                                                                                              <  1:560360/97‑1 (MQ=255)                                                                                cGATAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcaca                                                                                                             <  1:649741/79‑1 (MQ=255)                                                                                cGAGGGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcaca                                                                                                             >  2:649741/1‑79 (MQ=255)                                                                                 gagagaAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCtt                                                >  2:471085/1‑139 (MQ=255)                                                                                             aCTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCttt                                                                                                                       <  1:772804/56‑1 (MQ=255)                                                                                             aCTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCttt                                                                                                                       >  2:772804/1‑56 (MQ=255)                                                                                                         aaCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTaa                                                                                                                     <  1:39621/46‑1 (MQ=255)                                                                                                         aaCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTaa                                                                                                                     >  2:39621/1‑46 (MQ=255)                                                                                                         aaCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTgg                                                                                     <  1:328274/78‑1 (MQ=255)                                                                                                         aaCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTgg                                                                                     >  2:328274/1‑78 (MQ=255)                                                                                                         aaCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGaca                                     <  2:467352/126‑1 (MQ=255)                                                                                                         aaCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGaca                                     >  1:467352/1‑126 (MQ=255)                                                                                                          aCAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGa                       <  1:471085/139‑1 (MQ=255)                                                                                                              tACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGa                                                               >  2:998001/1‑95 (MQ=255)                                                                                                              tACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGa                                                               <  1:998001/95‑1 (MQ=255)                                                                                                                 cTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTc                <  1:883509/139‑1 (MQ=255)                                                                                                                   tATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCatat                                                                                                   >  2:828311/1‑54 (MQ=255)                                                                                                                   tATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCatat                                                                                                   <  1:828311/54‑1 (MQ=255)                                                                                                                   tATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGa                                                               <  1:287695/90‑1 (MQ=255)                                                                                                                   tATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGa                                                               >  2:287695/1‑90 (MQ=255)                                                                                                                     ttCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCATCGACCATGccc                                                                   >  1:997101/1‑84 (MQ=255)                                                                                                                     ttCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCATCGACCATGccc                                                                   <  2:997101/84‑1 (MQ=255)                                                                                                                     ttCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGt                                                        >  2:584342/1‑95 (MQ=255)                                                                                                                     ttCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGt                                                        <  1:584342/95‑1 (MQ=255)                                                                                                                          ttATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTcc                                         >  2:487551/1‑105 (MQ=255)                                                                                                                          ttATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTcc                                         <  1:487551/105‑1 (MQ=255)                                                                                                                          ttATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGaca                                     >  1:45959/1‑109 (MQ=255)                                                                                                                          ttATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGaca                                     <  2:45959/109‑1 (MQ=255)                                                                                                                           tATGCCCACAATTAGCGTAAATCTTTAACGGCCACATg                                                                                                           >  1:578363/1‑38 (MQ=255)                                                                                                                           tATGCCCACAATTAGCGCAAATCTTTAACGGCCACATg                                                                                                           <  2:578363/38‑1 (MQ=255)                                                                                                                               cccACAATTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGc  >  2:767989/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                    CAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAACTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGC  >  NC_000913/2983155‑2983419

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