Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,923,535	T→C	N58D (AAT→GAT)	yebY ←	DUF2511 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,923,535	0	T	C	100.0%	118.8 / NA	34	N58D (AAT→GAT)	yebY	DUF2511 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/17);  total (17/17)

CATTTTTTCGCCGGGATTGTTAGGATCATCAATCTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATTAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGCGCCGCTGCAGAACTGGT  >  NC_000913/1923407‑1923670                                                                                                                                   |                                                                                                                                        cATTTTTTCGCCGGGATTGTTAGGATCATCAATCTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATAc                                                                                                                               <  1:959253/139‑1 (MQ=255)        tCGCCGGGATTGTTAGGATCATCAATCCGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCAt                                                                                                                        >  1:679660/1‑139 (MQ=255)         cGCCGGGATTGTTAGGATCATCAATCTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCAtt                                                                                                                       >  2:620087/1‑139 (MQ=255)                       ggATCATCAATCCGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAg                                                                                                         <  2:679660/139‑1 (MQ=255)                            atcaATCTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACa                                                                                                                                  <  1:487322/109‑1 (MQ=255)                            atcaATCTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACa                                                                                                                                  >  2:487322/1‑109 (MQ=255)                                  cTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACa                                                                                                                              <  1:1002117/107‑1 (MQ=255)                                  cTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACa                                                                                                                              >  2:1002117/1‑107 (MQ=255)                                      aTCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCgg                                                                                                             >  1:467380/1‑120 (MQ=255)                                      aTCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCgg                                                                                                             <  2:467380/120‑1 (MQ=255)                                                    ttGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGcc                                                                            >  1:165699/1‑139 (MQ=255)                                                        gCATTGGTCTTCCCACTGGCAACTTCCTTTTGTGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCa                                                                          >  2:772514/1‑139 (MQ=255)                                                        gCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCAtt                                                                        <  2:693039/139‑1 (MQ=255)                                                            tGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTc                                                                    >  1:1033266/1‑139 (MQ=255)                                                             ggTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCt                                                                   <  2:1033266/139‑1 (MQ=255)                                                                 ttCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACa                                                                                                                              <  2:994859/76‑1 (MQ=255)                                                                 ttCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACa                                                                                                                              >  1:994859/1‑76 (MQ=255)                                                                 ttCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTAc                                                               >  1:1131157/1‑139 (MQ=255)                                                                         ggCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCa                                                       >  1:791420/1‑139 (MQ=255)                                                                         ggCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCa                                                       >  1:584987/1‑139 (MQ=255)                                                                           cAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCaaa                                                     <  2:165699/139‑1 (MQ=255)                                                                               ttCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGa                                                 <  2:1131157/139‑1 (MQ=255)                                                                                          tgCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTa                                        <  1:200213/139‑1 (MQ=255)                                                                                                    tttAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCg                              >  1:600654/1‑139 (MQ=255)                                                                                                     ttAAAGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCgg                             <  2:791420/139‑1 (MQ=255)                                                                                                         aGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGca                         >  1:785626/1‑139 (MQ=255)                                                                                                         aGGATACTGCACGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGca                         <  1:772514/139‑1 (MQ=255)                                                                                                                   aCGAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGcgccg               >  1:790335/1‑139 (MQ=255)                                                                                                                     gAGGGTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGCGCCGCt             >  2:91151/1‑139 (MQ=255)                                                                                                                        ggTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGt                                                                     <  2:856487/80‑1 (MQ=255)                                                                                                                        ggTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGt                                                                     >  1:856487/1‑80 (MQ=255)                                                                                                                        ggTACTCGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGCGCCGCTGCa          <  2:584987/139‑1 (MQ=255)                                                                                                                              cGGATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGCGCCGCTGCAGAACTg    <  2:600654/139‑1 (MQ=255)                                                                                                                                gATCAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGCGCCGCTGCAGAACTGGt  <  2:785626/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        CATTTTTTCGCCGGGATTGTTAGGATCATCAATCTGAATCACCGAAATGGGTTGGGCATTGGTCTTCCCACTGGCAACTTCCTTT‑‑TGTGCGATATCGTTTAAAGGATACTGCACGAGGGTACTCGGATTAATAACATACAAAGCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGACTGACGGTAATCACCTGCGGTGCAGCCAGCGCCGCTGCAGAACTGGT  >  NC_000913/1923407‑1923670

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