Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,950,099	C→T	intergenic (‑130/‑223)	yifB ← / → ilvL	putative magnesium chelatase YifB/ilvXGMEDA operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,950,099	0	C	T	100.0%	99.1 / NA	30	intergenic (‑130/‑223)	yifB/ilvL	putative magnesium chelatase YifB/ilvXGMEDA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/16);  total (14/16)

AAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTA  >  NC_000913/3949860‑3950278                                                                                                                                                                                                                                                |                                                                                                                                                                                    aaGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATg                                                                                                                                                                          <  2:397116/251‑1 (MQ=255)  aGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGc                                                                                                                                                                         >  2:339079/1‑251 (MQ=255)               tAAGCCTGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGCCAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCCAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAACAATGCAGCGGACACAgg                                                                                                                                                             >  2:265497/1‑250 (MQ=255)                   ccGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGa                                                                                                                                                         >  1:267703/1‑250 (MQ=255)                                tACTGATATGTACCTCAACAGTGATCGGGGACGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGAAATAAGCCCTACTTGAGCCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAggg                                                                                                                                           >  1:700748/1‑251 (MQ=255)                                    gATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGaac                                                                                                                                        <  2:371949/250‑1 (MQ=255)                                         gTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGTAACAACAt                                                                                                                                   <  1:465090/250‑1 (MQ=255)                                         gTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAAAAACAtt                                                                                                                                  >  2:491237/1‑251 (MQ=255)                                            ccTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCTGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGTCCGCTAATTCGTGAATTGTAGTTGCTGATTCCTGTTTATTTGGGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATTTTCGCAAAAATGCAGCGGACAAAGGAGGAACTACTAGGAAGGTTACAACATTCAt                                                                                                                               <  2:734020/251‑1 (MQ=255)                                                aaCAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTg                                                                                                                           >  1:685416/1‑251 (MQ=255)                                                  cAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGa                                                                                                                          >  1:278261/1‑250 (MQ=255)                                                    gTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGGATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAAt                                                                                                                       <  1:689875/251‑1 (MQ=255)                                                    gTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAAt                                                                                                                       <  1:47655/251‑1 (MQ=255)                                                       aTCGGGGTCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTTAGTTACCATTATGTGCATAAGATAGCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGTAATGATGTTCGCAAAAATGCATCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGa                                                                                                                    <  1:114566/251‑1 (MQ=255)                                                          gggggCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                                                                 <  1:12640/251‑1 (MQ=255)                                                            gggCGCATTTGCTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATAGCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGGTTATTGGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTGCTCAAAAATGGAGCGGACAAAGGAGGAACTACGAGGAAGGGAACAACATTCAGACTGAAATAGAAtttt                                                                                                               <  2:637626/251‑1 (MQ=255)                                                                     tACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTGTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACGGAAATTTAATTTTTGtcactc                                                                                                       <  1:129117/250‑1 (MQ=255)                                                                        ccccAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTTAATTTTTTTCACTCACta                                                                                                   <  1:554778/250‑1 (MQ=255)                                                                                     cgGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATAGCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                                      <  1:723782/251‑1 (MQ=255)                                                                                        gTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaa                                                                                    >  2:228522/1‑250 (MQ=255)                                                                                                aaTTGACAGTGACATAAGCCTTCCTTGAGTCACCATTATGTGCATAAGATAGCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTGCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                                                           <  1:93332/251‑1 (MQ=255)                                                                                                      cAGTGACATAAGCCCGCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATGTATAtt                                                                     <  2:463784/251‑1 (MQ=255)                                                                                                               aaGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTACTGTTTATTTGTGCAAGTGCAGTTGAGTTATTCTGGCGGTGGAATGATGTTCCCACAAATGCAGCGGACAAAGGATGAACTACGAGGGAGGGAAAAACATTCATAATACAATTGAATTTTTTTTACTCACTATTTTCTTTTTAAAAAAAAACAATTTATATTGAAattatt                                                            >  1:714136/1‑251 (MQ=255)                                                                                                                   cctcctTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc                                                        <  1:549081/251‑1 (MQ=255)                                                                                                                      cctTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACCACAATTTATATTGAAATTATTAAACGc                                                      >  1:328253/1‑250 (MQ=255)                                                                                                                          gAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGGTGTTCGCAAAACTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTCAAACACAACAATTTATATTGAAATTATTAAACGcatca                                                  >  1:713725/1‑250 (MQ=255)                                                                                                                           aGTCACCATTATGTGCATAAGATATCGCAGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTAATTCCTGTTTATTTGTGCATGCGAAGTTGAGTTGTTCTGGCGGTTGAATGATGTTCGCAAAAATGCAGCGGACAACGAATGAACTACGAGGAAGGGGACAACATTCCTACTGAAATTGAATTTTTTTCACTCACTCTTTTATTTTTTAAAAACACAAATTTTTATTGAAATTATTATACGcatcat                                                 >  2:273178/1‑250 (MQ=255)                                                                                                                                 cATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGACGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                                          >  1:686672/1‑251 (MQ=255)                                                                                                                                        tgCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGCGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaa                                    >  1:154501/1‑250 (MQ=255)                                                                                                                                                                         cGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAATCCTATGGTa  <  1:33195/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                    AAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTA  >  NC_000913/3949860‑3950278

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