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%	108.5 / NA	28	E13K (GAG→AAG)	cspC	stress protein, member of the CspA family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/16);  total (12/16)

TCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAA  >  NC_000913/1907353‑1907443                                               |                                             tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                            >  1:1760713/1‑49 (MQ=255)     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                           <  1:976642/46‑1 (MQ=255)     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                         <  1:34215/48‑1 (MQ=255)       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                      <  1:1861671/49‑1 (MQ=255)          ccATCAGCCGGAGTAATGAAGCCACAACCTTTAGACTTGTTGAACCACt                                   >  1:171337/1‑49 (MQ=38)           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                  <  1:1634480/49‑1 (MQ=255)           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                  <  1:1246381/49‑1 (MQ=255)                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                          >  1:907350/1‑49 (MQ=255)                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                          <  1:671426/49‑1 (MQ=255)                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                          <  1:2052762/49‑1 (MQ=255)                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                         >  1:1157644/1‑49 (MQ=255)                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                        <  1:1762241/49‑1 (MQ=255)                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                       <  1:973931/49‑1 (MQ=255)                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                      <  1:1690898/49‑1 (MQ=255)                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                       >  1:1423318/1‑48 (MQ=255)                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                      >  1:1720366/1‑48 (MQ=255)                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                 <  1:1336295/49‑1 (MQ=255)                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt               >  1:1295497/1‑49 (MQ=255)                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg              >  1:2105620/1‑49 (MQ=255)                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg              <  1:234067/49‑1 (MQ=255)                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg              <  1:117121/49‑1 (MQ=255)                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt         >  1:737673/1‑49 (MQ=255)                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt        <  1:235162/49‑1 (MQ=255)                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt        >  1:1284695/1‑49 (MQ=255)                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt        >  1:811187/1‑49 (MQ=255)                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa     >  1:49093/1‑48 (MQ=255)                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa  <  1:1884361/49‑1 (MQ=255)                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa   <  1:1584260/48‑1 (MQ=255)                                               |                                             TCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAA  >  NC_000913/1907353‑1907443

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