Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	551,332	G→A	*298* (TAG→TAA)	ybcF →	putative carbonate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	551,332	0	G	A	100.0%	122.8 / NA	39	*298* (TAG→TAA)	ybcF	putative carbonate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

ATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATG  >  NC_000913/551199‑551464                                                                                                                                      |                                                                                                                                     aTGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                 >  1:219896/1‑139 (MQ=255)     ttCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCa                                                                                                                             <  2:1192577/139‑1 (MQ=255)         aTGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAg                                                                                                                         >  2:978376/1‑139 (MQ=255)         aTGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAg                                                                                                                         <  1:978376/139‑1 (MQ=255)                        aCGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGAt                                                                                                          >  2:258607/1‑139 (MQ=255)                                        aTGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGt                                                                                          <  2:355943/139‑1 (MQ=255)                                             aGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAt                                                                                     >  1:243712/1‑139 (MQ=255)                                              gAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    <  1:258607/139‑1 (MQ=255)                                                aGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGc                                                                                                   <  1:599708/122‑1 (MQ=255)                                                aGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGc                                                                                                   >  2:599708/1‑122 (MQ=255)                                                     tGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGcc                                                                                                              <  1:366799/106‑1 (MQ=255)                                                     tGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGcc                                                                                                              >  2:366799/1‑106 (MQ=255)                                                               gcgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGc                                                                   >  2:70231/1‑139 (MQ=255)                                                                cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                 <  2:132351/124‑1 (MQ=255)                                                                cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                 >  1:132351/1‑124 (MQ=255)                                                                cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                  >  2:463555/1‑139 (MQ=255)                                                                 gTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                 <  2:602234/75‑1 (MQ=255)                                                                 gTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTa                                                                                                                                 >  1:602234/1‑75 (MQ=255)                                                                         gggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGc                                                         <  2:219896/139‑1 (MQ=255)                                                                           gCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAgg                                                                                                                               <  1:636988/67‑1 (MQ=255)                                                                           gCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAgg                                                                                                                               >  2:636988/1‑67 (MQ=255)                                                                                aTCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGt                                                  <  1:1259478/139‑1 (MQ=255)                                                                                    cgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAc                                                                         <  1:169933/112‑1 (MQ=255)                                                                                    cgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAc                                                                         >  2:169933/1‑112 (MQ=255)                                                                                        ttGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGATTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGCTGCACAATg                                          >  2:763712/1‑139 (MQ=255)                                                                                         tGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGcc                                                                                                              >  1:1095254/1‑70 (MQ=255)                                                                                         tGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGcc                                                                                                              <  2:1095254/70‑1 (MQ=255)                                                                                         tGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    >  1:769881/1‑96 (MQ=255)                                                                                         tGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    <  2:769881/96‑1 (MQ=255)                                                                                           aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    >  1:696223/1‑94 (MQ=255)                                                                                           aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                    <  2:696223/94‑1 (MQ=255)                                                                                                         ggcgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGa                         <  1:463555/139‑1 (MQ=255)                                                                                                                    gACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATc              >  2:960363/1‑139 (MQ=255)                                                                                                                             ttCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACa     <  2:243712/139‑1 (MQ=255)                                                                                                                               cGCTGTAATCGTAGGCATTAGACATTTGTGCGTGATGCGACGCTTGACGCGTCTTATCAGGCCGACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt               <  2:842143/127‑1 (MQ=255)                                                                                                                               cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   >  2:1134052/1‑139 (MQ=255)                                                                                                                               cGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt   <  1:960363/139‑1 (MQ=255)                                                                                                                               cGCTGTAAACGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt               >  1:842143/1‑127 (MQ=255)                                                                                                                                 cTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATg  <  1:199034/138‑1 (MQ=255)                                                                                                                                 cTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATg  >  2:199034/1‑138 (MQ=255)                                                                                                                                      |                                                                                                                                     ATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATG  >  NC_000913/551199‑551464

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