Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,398,384	C→T	intergenic (‑6/+3)	yhdE ← / ← mreD	nucleoside triphosphate pyrophosphatase YhdE/cell shape determining protein MreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,398,384	0	C	T	100.0%	99.0 / NA	29	intergenic (‑6/+3)	yhdE/mreD	nucleoside triphosphate pyrophosphatase YhdE/cell shape determining protein MreD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/19);  total (10/19)

CAGCACTTCTCCGTTCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGCGAGGTTGCGGAAAAGCTGGTATTTCAGCGCCACCAGGTAAGCAATGAT  >  NC_000913/3398139‑3398629                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cAGCACTTCTCCGTTCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTa                                                                                                                                                                                                                                                  <  2:160242/251‑1 (MQ=255)    cACTTCTCCGTTCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTg                                                                                                                                                                                                                                               <  1:91503/251‑1 (MQ=255)              ttCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAAc                                                                                                                                                                                                                                     >  1:110097/1‑251 (MQ=255)                            aGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACtt                                                                                                                                                                                                                        <  1:44349/250‑1 (MQ=255)                               aTCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACGGCTGACGGACTTTgcg                                                                                                                                                                                                                    <  1:212365/251‑1 (MQ=255)                                    cACCCAGCACCGGGAGAGCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCATAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCa                                                                                                                                                                                                               <  2:46182/251‑1 (MQ=255)                                         aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCa                                                                                                                                                                                                            >  2:16966/1‑249 (MQ=255)                                                    aTCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACgg                                                                                                                                                                                               <  1:85665/251‑1 (MQ=255)                                                        ttCGCAGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATAATGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGcc                                                                                                                                                                                             >  1:89690/1‑249 (MQ=255)                                                            ccGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCa                                                                                                                                                                                       <  1:141169/251‑1 (MQ=255)                                                                           cTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGactacac                                                                                                                                                                        >  1:188713/1‑251 (MQ=255)                                                                               ccGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTACGCGCTCTCCTGCGGCTGAAGCTGCTCCTCAATGCACGTAACAATACGTTCAAAGGTCACGCCACGTTGCGCAAGTAACTCCTGACGACGCGGAAAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCACGAAAATCCACGGCCAGAGCACCCCATTGACTACactact                                                                                                                                                                    >  1:168316/1‑251 (MQ=255)                                                                                          cgcgcgCCAGACGCACAACATACGGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAGCAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTGCACTGCTCCAGAACACtt                                                                                                                                                         <  2:159630/251‑1 (MQ=255)                                                                                            cgcgcCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTcc                                                                                                                                                       <  2:48519/251‑1 (MQ=255)                                                                                                         caacaTACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAAGACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAgaga                                                                                                                                          <  2:132163/251‑1 (MQ=255)                                                                                                            caTACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGt                                                                                                                                       <  1:44337/251‑1 (MQ=255)                                                                                                             aTACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGtt                                                                                                                                      >  1:166163/1‑251 (MQ=255)                                                                                                                       gcgcTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTaa                                                                                                                             >  1:210445/1‑250 (MQ=255)                                                                                                                              ccTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACTGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAAGCACTAAAAActct                                                                                                                      <  2:87284/250‑1 (MQ=255)                                                                                                                                          gctgctCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACaa                                                                                                         <  2:161573/251‑1 (MQ=255)                                                                                                                                                   cAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGGAACTCCTGACTACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTAGTGCACTGCAAACTGCTGACGGCCTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTGCACTACTCCATAACGCTTCCGGTCTGAAAGAGACCTTAATCACGAAAAACTCTGCCCAGAAAGCAAAAATATcca                                                                                                <  2:145856/251‑1 (MQ=255)                                                                                                                                                       gCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATccacca                                                                                             <  1:137525/250‑1 (MQ=255)                                                                                                                                                           gTAACAATAAGTTCAAAGGTCACGCCAAGTTGCGCCAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGc                                                                                        >  1:83172/1‑251 (MQ=255)                                                                                                                                                                          aaGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACCTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgacg                                                                         <  1:12019/251‑1 (MQ=255)                                                                                                                                                                                                gTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATg                                                   >  1:93131/1‑251 (MQ=255)                                                                                                                                                                                                         gacgacGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAATCTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGCGAGGtt                                            <  2:41342/249‑1 (MQ=255)                                                                                                                                                                                                             acgCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGCGAGGTTGCGGa                                       >  1:61008/1‑250 (MQ=255)                                                                                                                                                                                                                 ggAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGCGAGGTTGCGGAAAAGc                                  <  1:241819/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aGAAATCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGCGAGGTTGCGGAAAAGCTGGTATTTCAGCGCCACCAGGTAAGCAatgat  <  2:110097/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CAGCACTTCTCCGTTCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAATGCGAGGTTGCGGAAAAGCTGGTATTTCAGCGCCACCAGGTAAGCAATGAT  >  NC_000913/3398139‑3398629

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