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

TGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907357‑1907445                                           |                                               tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             <  1:1721622/46‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                          <  1:284664/49‑1 (MQ=255)   ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:2471822/49‑1 (MQ=255)      ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      <  1:832955/48‑1 (MQ=255)      ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    <  1:558467/50‑1 (MQ=255)      ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     <  1:1037752/49‑1 (MQ=255)       cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                   <  1:777116/50‑1 (MQ=255)        aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  >  1:1409478/1‑50 (MQ=255)          cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                >  1:893470/1‑50 (MQ=255)          cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                >  1:1555550/1‑50 (MQ=255)           aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                <  1:407468/49‑1 (MQ=255)           aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                <  1:1906561/49‑1 (MQ=255)           aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:1887771/1‑50 (MQ=255)              cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            <  1:1300581/50‑1 (MQ=255)                gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:1440811/50‑1 (MQ=255)                gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:2350245/50‑1 (MQ=255)                gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:444458/1‑50 (MQ=255)                gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:1951691/1‑50 (MQ=255)                 aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         <  1:1897522/50‑1 (MQ=255)                 aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:1114263/1‑50 (MQ=255)                  gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        <  1:2702148/50‑1 (MQ=255)                  gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:1384713/1‑49 (MQ=255)                   tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:2341574/1‑49 (MQ=255)                   tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:103853/1‑49 (MQ=255)                   tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                       <  1:338751/50‑1 (MQ=255)                     aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:910367/49‑1 (MQ=255)                     aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     <  1:440678/50‑1 (MQ=255)                        aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:70598/50‑1 (MQ=255)                         aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:2059176/50‑1 (MQ=255)                          gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                >  1:2540463/1‑50 (MQ=255)                             aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:1684161/1‑50 (MQ=255)                                aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:2071664/50‑1 (MQ=255)                                aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:72265/50‑1 (MQ=255)                                  cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        <  1:1508550/50‑1 (MQ=255)                                     tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      >  1:1400402/1‑49 (MQ=255)                                     tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      >  1:928945/1‑49 (MQ=255)                                     tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa       <  1:1700475/48‑1 (MQ=255)                                     tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa       <  1:1410598/48‑1 (MQ=255)                                        aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  <  1:1743214/50‑1 (MQ=255)                                           |                                               TGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907357‑1907445

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