Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,907,399	C→T	E13K (GAG→AAG)	cspC ←	stress protein, member of the CspA family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,907,399	0	C	T	100.0%	94.7 / NA	27	E13K (GAG→AAG)	cspC	stress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/16);  total (11/16)

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAA  >  NC_000913/1907352‑1907443                                                |                                             aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                            <  1:2038737/50‑1 (MQ=255)      tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                           <  1:274382/46‑1 (MQ=255)        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                      <  1:2119232/49‑1 (MQ=255)        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                     <  1:1428641/50‑1 (MQ=255)        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                     <  1:1545847/50‑1 (MQ=255)           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                    <  1:1188273/48‑1 (MQ=255)            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                  >  1:139403/1‑49 (MQ=255)                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                              <  1:1190089/49‑1 (MQ=255)                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                              <  1:2079041/49‑1 (MQ=255)                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                             <  1:1693433/50‑1 (MQ=255)                  ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                           >  1:1826813/1‑50 (MQ=255)                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                          >  1:362595/1‑50 (MQ=255)                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                          >  1:1462344/1‑50 (MQ=255)                    ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                         >  1:917590/1‑50 (MQ=255)                       gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                      <  1:1881805/50‑1 (MQ=255)                          aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                   >  1:1389839/1‑50 (MQ=255)                           tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                  >  1:2245970/1‑50 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt               <  1:272247/50‑1 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt               >  1:1348382/1‑50 (MQ=255)                               gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg              <  1:1463677/50‑1 (MQ=255)                                  aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa           >  1:1940206/1‑50 (MQ=255)                                    aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt         <  1:1370470/50‑1 (MQ=255)                                     aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt        <  1:1905823/50‑1 (MQ=255)                                      ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt        >  1:900520/1‑49 (MQ=255)                                       cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa      <  1:1460787/50‑1 (MQ=255)                                           aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa  >  1:2261889/1‑50 (MQ=255)                                           aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa   <  1:95602/49‑1 (MQ=255)                                                |                                             ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAA  >  NC_000913/1907352‑1907443

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