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%	158.4 / NA	50	*399* (TAG→TAA)	hemY	putative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/26);  total (24/26)

ATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGT  >  NC_000913/3986552‑3986812                                                                                                                                       |                                                                                                                               aTGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAattact                                                                                                                            >  1:297348/1‑139 (MQ=255)         aaCGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                                    <  1:1356993/139‑1 (MQ=255)         aaCGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                                    <  1:221825/139‑1 (MQ=255)                  tGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           <  2:460074/139‑1 (MQ=255)                  tGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCGGCTCCGTAATAAGAGCAGGCGGTTTTTTATGTATCAGGAAGGCCCAGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGc                                                                                                           <  2:297348/139‑1 (MQ=255)                                   gtCATACAGACCTGTTTTAGCGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAAcc                                                                                          >  2:1189647/1‑139 (MQ=255)                                   gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGtt                                                                                                     <  1:337894/128‑1 (MQ=255)                                   gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGtt                                                                                                     >  2:337894/1‑128 (MQ=255)                                   gtCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAAcc                                                                                          >  1:462518/1‑139 (MQ=255)                                                 ttttAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCa                                                                               <  1:1497670/136‑1 (MQ=255)                                                 ttttAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCa                                                                               >  2:1497670/1‑136 (MQ=255)                                                                 gTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCgg                                                            <  1:1189647/139‑1 (MQ=255)                                                                  taataaGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                           <  2:1134422/139‑1 (MQ=255)                                                                       agagCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTtgt                                                       <  1:1362222/138‑1 (MQ=255)                                                                       agagCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTtgt                                                       >  2:1362222/1‑138 (MQ=255)                                                                          gCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                     <  2:236698/121‑1 (MQ=255)                                                                          gCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                     >  1:236698/1‑121 (MQ=255)                                                                             ggCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                     >  2:117147/1‑118 (MQ=255)                                                                             ggCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                     <  1:117147/118‑1 (MQ=255)                                                                                            aTCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCa                                                                                                          >  2:810315/1‑66 (MQ=255)                                                                                            aTCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCa                                                                                                          <  1:810315/66‑1 (MQ=255)                                                                                                    ggCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATg                         >  2:163314/1‑139 (MQ=255)                                                                                                     gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCa                                                                                              >  1:1485724/1‑69 (MQ=255)                                                                                                     gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCa                                                                                              <  2:1485724/69‑1 (MQ=255)                                                                                                     gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  2:1446729/1‑139 (MQ=255)                                                                                                     gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  2:1321473/1‑139 (MQ=255)                                                                                                     gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  1:981635/1‑139 (MQ=255)                                                                                                      ccccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCg                       <  2:981635/139‑1 (MQ=255)                                                                                                       cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGt                      >  1:189790/1‑139 (MQ=255)                                                                                                       cccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGt                      >  1:460401/1‑139 (MQ=255)                                                                                                           gAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGtt                                                                                                     <  1:938133/56‑1 (MQ=255)                                                                                                           gAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGtt                                                                                                     >  2:938133/1‑56 (MQ=255)                                                                                                                gCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGt                                                                                                                   >  1:718550/1‑37 (MQ=255)                                                                                                                gCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGt                                                                                                                   <  2:718550/37‑1 (MQ=255)                                                                                                                gCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        <  2:639346/128‑1 (MQ=255)                                                                                                                gCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  1:639346/1‑128 (MQ=255)                                                                                                                 cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                 <  1:562407/70‑1 (MQ=255)                                                                                                                 cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                 >  2:562407/1‑70 (MQ=255)                                                                                                                 cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        >  1:1453854/1‑127 (MQ=255)                                                                                                                 cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                        <  2:1453854/127‑1 (MQ=255)                                                                                                                   tGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCAt                                                                              >  1:473063/1‑71 (MQ=255)                                                                                                                   tGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCAt                                                                              <  2:473063/71‑1 (MQ=255)                                                                                                                        ccGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGa     >  2:1619097/1‑139 (MQ=255)                                                                                                                        ccGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGa     <  1:1619097/139‑1 (MQ=255)                                                                                                                        ccGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGa     <  1:1321473/139‑1 (MQ=255)                                                                                                                           ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGc                                                     <  1:220305/88‑1 (MQ=255)                                                                                                                           ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGc                                                     >  2:220305/1‑88 (MQ=255)                                                                                                                           ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACgt  <  1:1446729/139‑1 (MQ=255)                                                                                                                               agaAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCgg                                                            <  1:1318951/77‑1 (MQ=255)                                                                                                                               agaAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCgg                                                            >  2:1318951/1‑77 (MQ=255)                                                                                                                                    gAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACgt  <  2:1309046/130‑1 (MQ=255)                                                                                                                                       |                                                                                                                               ATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGT  >  NC_000913/3986552‑3986812

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: ‑