Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,499,158	T→C	V444A (GTT→GCT)	cysG →	siroheme synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,499,158	0	T	C	100.0%	107.5 / NA	32	V444A (GTT→GCT)	cysG	siroheme synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/19);  total (13/19)

GCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATT  >  NC_000913/3499029‑3499286                                                                                                                                  |                                                                                                                                 gCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGgctgt                                                                                                                               >  1:21951/1‑133 (MQ=255)            gTGGCAATTGGTGAAAACGGTACGACAGTCACACAGCGCGTTATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGcgcg                                                                                                                    <  2:193184/133‑1 (MQ=255)               gCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATa                                                                                                                 <  2:39142/133‑1 (MQ=255)                 aaTTGTCGTAAACGGGACTGCAATCACGCACCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATaaa                                                                                                               <  1:198270/133‑1 (MQ=255)                          aaaCGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTgg                                                                                                      <  1:99015/133‑1 (MQ=255)                           aaCGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGt                                                                                                     >  2:468998/1‑133 (MQ=255)                                tACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTATTGGCCTGCTCGATAAACTGAACTGGTTCTcc                                                                                                <  2:355926/133‑1 (MQ=255)                                      aGTCACGCAGCTCGTGATTGACGGTACGCTCACACAGCTGGGATAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAGCTGGTTCTCCAACCAt                                                                                          <  2:123872/133‑1 (MQ=255)                                                cgcgTGATGGACGGTCCGCTCATACAGCTGTGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAAcc                                                                                <  1:164098/133‑1 (MQ=255)                                                  cgTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAATTGAACTGGTTCTCCAAACATTGATTTAACCCg                                                                              >  2:164098/1‑133 (MQ=255)                                                  cgTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCg                                                                              >  1:443139/1‑133 (MQ=255)                                                                            tGGGAGAACTAGCGCTGCAAATGATCAGACCGACGCTAATTCTTATTGGTCGGGCTGTTGGCCTGCGCAATAAACTGAACTGGTTCTCCAACCATTAATTTACCCCGGGCCAGAGAATTCTGGCCTTCTTAAt                                                    <  2:163893/133‑1 (MQ=255)                                                                             gggAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATg                                                   <  2:87424/133‑1 (MQ=255)                                                                               gagaACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGAc                                                 <  2:332188/133‑1 (MQ=255)                                                                                agaACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGAct                                                <  1:93895/133‑1 (MQ=255)                                                                                            gcaAATGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCTAGAGAATTCTGGCCTTCTTAATGTCTCTTTTTatatat                                    >  2:320460/1‑133 (MQ=255)                                                                                               attGAACAGTTCATCGCTAATTATTATTCGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCa                                 <  2:244901/131‑1 (MQ=255)                                                                                               aaTGAACAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCa                                 <  2:443139/133‑1 (MQ=255)                                                                                               aaTGAACACTCCATCCCTAATTATTATTGGGCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCa                                 <  1:187135/133‑1 (MQ=255)                                                                                                 tGAACAGTACATCGCAAATTATTATTGGTCTGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGc                               <  1:247805/133‑1 (MQ=255)                                                                                                     cAGTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAAt                           <  2:21951/133‑1 (MQ=255)                                                                                                       gTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATaa                         >  1:419470/1‑133 (MQ=255)                                                                                                       gTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATa                          >  2:464859/1‑132 (MQ=255)                                                                                                       gTCCATCGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATa                          <  1:464859/132‑1 (MQ=255)                                                                                                             cGCTAATTATTATTGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGAGTCTTTTTGTATATTGAGCAAAGAAAACATa                   >  1:134965/1‑133 (MQ=255)                                                                                                                  attattattGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGcccc              >  1:111998/1‑133 (MQ=255)                                                                                                                  attattattGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGcccc              >  2:37417/1‑133 (MQ=255)                                                                                                                  attattattGCTCGGCCTGTTGGCATGCTCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGcccc              <  2:156251/133‑1 (MQ=255)                                                                                                                      ttattGGTCGTGCTGCTGGCCTGCTCGACAAATTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTtaa          <  2:134965/133‑1 (MQ=255)                                                                                                                      ttattGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTtaa          <  2:419470/133‑1 (MQ=255)                                                                                                                        attGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTtaata        >  2:303755/1‑133 (MQ=255)                                                                                                                         ttGGTCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTtaataa       >  2:96559/1‑133 (MQ=255)                                                                                                                             tCGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTGTTATATATGCAGTAAAGAAAACATATCCCCGTAATAAatat   >  2:142648/1‑133 (MQ=255)                                                                                                                              cGGGCTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAAtatt  >  1:107977/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                                 GCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGATTGACGGTACGCTCACACAGCTGGGAGAACTGGCGCAGCAAATGAACAGTCCATCGCTAATTATTATTGGTCGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATT  >  NC_000913/3499029‑3499286

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