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%	138.9 / 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 (21/14);  total (21/14)

CTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907354‑1907444                                              |                                              cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            >  1:1065758/1‑49 (MQ=255)  tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                            >  1:1617489/1‑48 (MQ=255)      ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                       <  1:2011439/49‑1 (MQ=255)      ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                        >  1:1166430/1‑48 (MQ=255)      ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                        <  1:1312184/48‑1 (MQ=255)         ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    >  1:759432/1‑49 (MQ=255)         ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGCTGAACCAc                                     <  1:1388486/48‑1 (MQ=37)          cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                    >  1:629380/1‑48 (MQ=255)            tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                 >  1:31192/1‑49 (MQ=255)              aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                               >  1:89455/1‑49 (MQ=255)              aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                               >  1:614737/1‑49 (MQ=255)              aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                <  1:608529/48‑1 (MQ=255)              aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                <  1:888912/48‑1 (MQ=255)                 cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                            >  1:1655192/1‑49 (MQ=255)                   gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          >  1:573189/1‑49 (MQ=255)                   gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          <  1:378826/49‑1 (MQ=255)                   gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                          <  1:997719/49‑1 (MQ=255)                    aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                         >  1:1056425/1‑49 (MQ=255)                        aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                     >  1:1134159/1‑49 (MQ=255)                           aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  <  1:1937042/49‑1 (MQ=255)                           aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                  <  1:737308/49‑1 (MQ=255)                             gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                <  1:2066647/49‑1 (MQ=255)                               cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              >  1:2028438/1‑49 (MQ=255)                               cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              >  1:1767200/1‑49 (MQ=255)                               cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc              >  1:1057482/1‑49 (MQ=255)                                 aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa            >  1:464113/1‑49 (MQ=255)                                   aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt          <  1:1613216/49‑1 (MQ=255)                                    ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt         >  1:1912728/1‑49 (MQ=255)                                     cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt         >  1:995329/1‑48 (MQ=255)                                     cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg        <  1:1590580/49‑1 (MQ=255)                                       ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa      >  1:1926743/1‑49 (MQ=255)                                        tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa       <  1:738126/47‑1 (MQ=255)                                        tAGACTTGTTGAACCACTTAACATGACCTTTAATCTTTGCCATTTGaa      >  1:201296/1‑48 (MQ=37)                                           aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  >  1:599416/1‑49 (MQ=255)                                           aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt  <  1:758387/49‑1 (MQ=255)                                              |                                              CTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT  >  NC_000913/1907354‑1907444

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