Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,986,686	C→T	*399* (TAG→TAA)	hemY ←	putative protoheme IX synthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,986,686	0	C	T	100.0%	139.9 / NA	45	*399* (TAG→TAA)	hemY	putative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/17);  total (28/17)

ACCATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGT  >  NC_000913/3986549‑3986812                                                                                                                                          |                                                                                                                               acCATGGACACAACGTTGAGTGAAGCACCC‑CTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAatta                                                                                                                              >  2:140893/1‑139 (MQ=255)   cATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCAGTAATAAGAGCAGGCGCTTTTTTATGTATCAGGAAGGCACCGGAGGTGCTTGCCTCCGGGTGAGAAGGAattac                                                                                                                             >  1:1077414/1‑139 (MQ=255)     tGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGGGAAGGAATTACtg                                                                                                                           >  2:165962/1‑139 (MQ=255)     tGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtg                                                                                                                           >  1:1048511/1‑139 (MQ=255)              cGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtgtg                                                                                                                         >  1:812917/1‑132 (MQ=255)              cGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtgtg                                                                                                                         <  2:812917/132‑1 (MQ=255)                  gAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTc                                                                                                              >  2:269477/1‑139 (MQ=255)                     tGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGTGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           <  2:282060/139‑1 (MQ=255)                                 ttgttgTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           >  1:175161/1‑127 (MQ=255)                                 ttgttgTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           <  2:175161/127‑1 (MQ=255)                                      gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtgtg                                                                                                                         >  2:314322/1‑108 (MQ=255)                                      gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtgtg                                                                                                                         <  1:314322/108‑1 (MQ=255)                                      gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTg                                                                                                            >  2:488984/1‑121 (MQ=255)                                      gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTg                                                                                                            <  1:488984/121‑1 (MQ=255)                                      gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAAcc                                                                                          >  1:674054/1‑139 (MQ=255)                                       tCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGc                                                                                                                       <  2:1148197/109‑1 (MQ=255)                                       tCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTATTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGc                                                                                                                       >  1:1148197/1‑109 (MQ=255)                                               aCCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                 >  2:971335/1‑139 (MQ=255)                                               aCCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                 >  1:650044/1‑139 (MQ=255)                                                 cTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCa                                                                               >  2:1053716/1‑139 (MQ=255)                                                              tGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                 >  1:515720/1‑124 (MQ=255)                                                              tGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                 <  2:515720/124‑1 (MQ=255)                                                                     taataaGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                           >  2:1104100/1‑139 (MQ=255)                                                                             gCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGGTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAg                                                   >  1:927297/1‑139 (MQ=255)                                                                              cAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                  <  1:140893/139‑1 (MQ=255)                                                                              cAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCGCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                  <  1:565700/139‑1 (MQ=255)                                                                               aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAatta                                                                                                                              >  1:249889/1‑62 (MQ=255)                                                                               aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAatta                                                                                                                              <  2:249889/62‑1 (MQ=255)                                                                               aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                      >  2:927900/1‑102 (MQ=255)                                                                               aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                      <  1:927900/102‑1 (MQ=255)                                                                                ggCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCt                                                >  1:295913/1‑139 (MQ=255)                                                                                              tATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCgg                                  >  1:585077/1‑139 (MQ=255)                                                                                                    gAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGcc                            <  2:295913/139‑1 (MQ=255)                                                                                                      aGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCAt                          >  1:1140679/1‑139 (MQ=255)                                                                                                        gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           <  2:1150123/56‑1 (MQ=255)                                                                                                        gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           >  1:1150123/1‑56 (MQ=255)                                                                                                        gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  1:268328/1‑139 (MQ=255)                                                                                                        gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  1:370899/1‑139 (MQ=255)                                                                                                        gCCCCGGAGGTGCTTGCCTCCGGGTAAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  1:905359/1‑139 (MQ=255)                                                                                                         ccccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCg                       <  2:370899/139‑1 (MQ=255)                                                                                                         ccccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCg                       <  2:1140679/139‑1 (MQ=255)                                                                                                          cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGt                      >  2:321131/1‑139 (MQ=255)                                                                                                                          tCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGcc                            <  1:503309/117‑1 (MQ=255)                                                                                                                          tCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGcc                            >  2:503309/1‑117 (MQ=255)                                                                                                                          tCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCgg      >  1:641656/1‑139 (MQ=255)                                                                                                                           ccGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCCTAGCTGCAGCATCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGa     >  1:132423/1‑139 (MQ=255)                                                                                                                            cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                      <  2:229868/57‑1 (MQ=255)                                                                                                                            cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                      >  1:229868/1‑57 (MQ=255)                                                                                                                              ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACgt  <  2:268328/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                               ACCATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGT  >  NC_000913/3986549‑3986812

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