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

CTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907354‑1907444                                              |                                              cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            >  1:1735286/1‑49 (MQ=255)     gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                        >  1:2181500/1‑49 (MQ=255)       tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                      >  1:716075/1‑49 (MQ=255)       tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                      >  1:382356/1‑49 (MQ=255)        gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                     >  1:1231109/1‑49 (MQ=255)        gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                     >  1:2203175/1‑49 (MQ=255)         ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    <  1:704566/49‑1 (MQ=255)          cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                   <  1:1129540/49‑1 (MQ=255)           aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                   >  1:1483991/1‑48 (MQ=255)           aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                  <  1:1747110/49‑1 (MQ=255)               gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                              <  1:2565926/49‑1 (MQ=255)                ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                              >  1:1238977/1‑48 (MQ=255)                ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                             <  1:2649456/49‑1 (MQ=255)                ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                             <  1:2783992/49‑1 (MQ=255)                 cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:2268453/1‑49 (MQ=255)                 cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:703281/1‑49 (MQ=255)                 cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:2118374/1‑49 (MQ=255)                  ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                           <  1:963496/49‑1 (MQ=255)                   gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          >  1:2715486/1‑49 (MQ=255)                    aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                         >  1:861368/1‑49 (MQ=255)                     gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                        >  1:1696468/1‑49 (MQ=255)                      tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                       <  1:1519115/49‑1 (MQ=255)                        aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                     >  1:939427/1‑49 (MQ=255)                        aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                     <  1:1859200/49‑1 (MQ=255)                           aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  <  1:903629/49‑1 (MQ=255)                             gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                >  1:1020125/1‑49 (MQ=255)                                aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              >  1:729456/1‑48 (MQ=255)                                   aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt          >  1:1791774/1‑49 (MQ=255)                                    ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt         <  1:2354523/49‑1 (MQ=255)                                        tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa      >  1:824575/1‑48 (MQ=255)                                         aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa     >  1:761392/1‑48 (MQ=255)                                          gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa   >  1:2858518/1‑49 (MQ=255)                                          gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa    <  1:396661/48‑1 (MQ=255)                                          gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa    <  1:2204436/48‑1 (MQ=255)                                           aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  >  1:991059/1‑49 (MQ=255)                                              |                                              CTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907354‑1907444

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