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

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445                                                |                                               aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              >  1:1854385/1‑50 (MQ=255)  tCTTTGCTGCCATCAGCCGGAGTAATGAAGACAAAACCTTTAGACTTGtt                                             >  1:325109/1‑50 (MQ=25)        ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                       >  1:1543738/1‑50 (MQ=255)           ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     <  1:396964/49‑1 (MQ=255)            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    >  1:187846/1‑49 (MQ=255)            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    >  1:1487513/1‑49 (MQ=255)            cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    >  1:1667819/1‑49 (MQ=255)              tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                 >  1:1200582/1‑50 (MQ=255)               cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                >  1:95613/1‑50 (MQ=255)                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  1:1330885/50‑1 (MQ=255)                aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:1498046/1‑50 (MQ=255)                 gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  1:904836/50‑1 (MQ=255)                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            >  1:1609472/1‑50 (MQ=255)                   cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            >  1:764655/1‑50 (MQ=255)                    ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  1:396500/1‑50 (MQ=255)                    ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           <  1:2237744/50‑1 (MQ=255)                     gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:433090/50‑1 (MQ=255)                      aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         <  1:1224927/50‑1 (MQ=255)                          aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:290621/49‑1 (MQ=255)                             aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:1467080/50‑1 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:495983/50‑1 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:2129426/50‑1 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:430540/50‑1 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 >  1:1235649/1‑50 (MQ=255)                              aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 >  1:1010656/1‑50 (MQ=255)                                 cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc              <  1:250295/50‑1 (MQ=255)                                 cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:1075707/1‑49 (MQ=255)                                  aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:191983/1‑50 (MQ=255)                                    aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           <  1:2221771/50‑1 (MQ=255)                                    aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           <  1:2200228/50‑1 (MQ=255)                                    aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           >  1:489805/1‑50 (MQ=255)                                       cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        <  1:1691010/50‑1 (MQ=255)                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      <  1:780310/49‑1 (MQ=255)                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa       <  1:1744580/48‑1 (MQ=255)                                          tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa       <  1:1722127/48‑1 (MQ=255)                                           aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa     <  1:1374848/49‑1 (MQ=255)                                             aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  >  1:400675/1‑50 (MQ=255)                                                |                                               ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445

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