Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,006,859	A→C	S200R (AGC→CGC)	ygeW →	putative carbamoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,006,859	0	A	C	100.0%	155.0 / NA	42	S200R (AGC→CGC)	ygeW	putative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/21);  total (21/21)

GGCGCAGGCAACGCCTATATGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTT  >  NC_000913/3006643‑3007062                                                                                                                                                                                                                         |                                                                                                                                                                                                            ggCGCAGGCAACGCCTATATGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATg                                                                                                                                                                                                       >  1:347948/1‑223 (MQ=255)   cgcAGGCAACGCCTATATGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTAGTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAAcc                                                                                                                                                                                                <  2:642148/228‑1 (MQ=255)        gCAACGCCGAAAGGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                                                                                                                                           <  1:476801/228‑1 (MQ=255)           aCGCCTATATGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGTCCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTg                                                                                                                                                                                        <  2:135472/228‑1 (MQ=255)                  tatGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTAGGGCAAACCGc                                                                                                                                                                                              >  1:227545/1‑215 (MQ=255)                  attGCGTGAAGTTGGCCCTGCATTTGAGGACGGTTACAAGCTGGGTTTACTGCCCCAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGct                                                                                                                                                                                             <  2:227198/214‑1 (MQ=255)                                                cGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGca                                                                                                                                                                            >  1:292399/1‑203 (MQ=255)                                                cGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCg                                                                                                                                                    >  2:284714/1‑227 (MQ=255)                                                 ggTTACAAGCAGGGTNTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGcat                                                                                                                                                                           <  2:291929/203‑1 (MQ=255)                                                             ggTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctctct                                                                                                                                                                                         >  1:126808/1‑177 (MQ=255)                                                              gtgtACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTg                                                                                                                                                                                        <  2:126622/177‑1 (MQ=255)                                                                       cacaGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATc                                                                                                                            >  1:205961/1‑228 (MQ=255)                                                                       cacaGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATc                                                                                                                            >  1:598451/1‑228 (MQ=255)                                                                             gTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAg                                                                                                                      <  2:597500/228‑1 (MQ=255)                                                                              tCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACca                                                                                                                                                                                   >  1:10196/1‑166 (MQ=255)                                                                              tCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTAt                                                                                                                                                >  1:65985/1‑201 (MQ=255)                                                                               ccGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACcac                                                                                                                                                                                  <  2:10184/166‑1 (MQ=255)                                                                               ccGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATg                                                                                                                                               <  2:65890/201‑1 (MQ=255)                                                                                        gTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctc                                                                                                                                                                                            >  2:441128/1‑147 (MQ=255)                                                                                         tGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                                                                                                                                           <  1:441847/147‑1 (MQ=255)                                                                                                   aTGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAATACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTa                                                                                                                                                 >  2:666182/1‑179 (MQ=255)                                                                                                    tGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTAt                                                                                                                                                <  1:667267/179‑1 (MQ=255)                                                                                                        atatTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTg                                                                                           >  1:6296/1‑228 (MQ=255)                                                                                                                                   cTGACCTGGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGcc                                                                 <  2:205646/227‑1 (MQ=255)                                                                                                                                      aCCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAAcc                                                                                                                                                                                                >  1:654496/1‑97 (MQ=255)                                                                                                                                       ccTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATAACCTGGGCCTACTCTCCACGCTATGGCAAACCg                                                                                                                                                                                               <  2:653437/97‑1 (MQ=255)                                                                                                                                                           acTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGT‑ACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCAc                                       <  1:285166/228‑1 (MQ=255)                                                                                                                                                             tttGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTAc                                                                                                                 >  1:586306/1‑153 (MQ=255)                                                                                                                                                              ttGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACg                                                                                                                <  2:585366/153‑1 (MQ=255)                                                                                                                                                               tGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGc                                                                                                                                                                                              >  1:155639/1‑74 (MQ=255)                                                                                                                                                                ggTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGct                                                                                                                                                                                             <  2:155414/74‑1 (MQ=255)                                                                                                                                                                       tGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAg                                                                   >  1:423062/1‑189 (MQ=255)                                                                                                                                                                        ggAAACCCTGAAAGGTACACAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGc                                                                  <  2:422373/189‑1 (MQ=255)                                                                                                                                                                                         aaaaaaTCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTAcc                                                                                                                                                                                    >  1:140619/1‑58 (MQ=255)                                                                                                                                                                                          aaaaaTCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACca                                                                                                                                                                                   <  2:140415/58‑1 (MQ=255)                                                                                                                                                                                                 gCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGtt  >  2:446914/1‑228 (MQ=255)                                                                                                                                                                                                                 aCTCTCCACGCTATGGCAAACCGCTCTCTGTACCACaa                                                                                                                                                                                >  2:347970/1‑38 (MQ=255)                                                                                                                                                                                                                 aCTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGc                                                                  >  2:402625/1‑148 (MQ=255)                                                                                                                                                                                                                  ctctCCACGCTATGGCAAACCGCTCTCTGTACCACAAg                                                                                                                                                                               <  1:348556/38‑1 (MQ=255)                                                                                                                                                                                                                  ctctCCACGCTATGGCAAACCGCTCTCTGTACCACAAg                                                                                                                                                                               >  1:431360/1‑38 (MQ=255)                                                                                                                                                                                                                  ctctCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGcc                                                                 <  1:403283/148‑1 (MQ=255)                                                                                                                                                                                                                   tctcCACGCTATGGCAAACCGCTCTCTGTACCACAAgg                                                                                                                                                                              <  2:430658/38‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                            GGCGCAGGCAACGCCTATATGCGTGAAGTTGGCGCTGCACTTGATGACGGTTACAAGCAGGGTGTACTGCCACAGCGTCCGGCTTTAGTGAACCTGCAATGCGATATTGACCACCCGACTCAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGACCTGATCCCGGATGTGGTTGAAGTGGCGAAAAACAATGCTAAAGCCTCCGGTGGTAGCTTCCGTCAGGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTT  >  NC_000913/3006643‑3007062

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