Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,368,168	T→C	R30R (CGT→CGC)	pspA →	phage shock protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,368,168	0	T	C	100.0%	60.8 / NA	19	R30R (CGT→CGC)	pspA	phage shock protein A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/12);  total (7/12)

GCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGA  >  NC_000913/1368039‑1368296                                                                                                                                  |                                                                                                                                 gcagGACAGTCCTGATCGTAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCtg                                                                                                                               <  1:242553/133‑1 (MQ=255)  cagGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCtga                                                                                                                              <  1:87593/133‑1 (MQ=255)   agGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAGCATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCtgat                                                                                                                             <  1:328232/133‑1 (MQ=255)    ggACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCtgatg                                                                                                                            <  2:95565/133‑1 (MQ=255)       cAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATc                                                                                                                         <  1:51283/133‑1 (MQ=255)           ccTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAgg                                                                                                                     <  1:21461/133‑1 (MQ=255)                    aGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAg                                                                                                            >  1:66130/1‑133 (MQ=255)                           aaGAGGACCACATTATGGTTATGTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACt                                                                                                     <  1:249310/133‑1 (MQ=255)                                   aacaTTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTgg                                                                                                   >  1:89165/1‑127 (MQ=255)                                   aacaTTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTgg                                                                                                   <  2:89165/127‑1 (MQ=255)                                     caTTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTa                                                                                           <  2:66130/133‑1 (MQ=255)                                         aTGGGTATTTTTCCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGtt                                                                                       <  1:484252/133‑1 (MQ=255)                                               aTTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGAGCCAGGAGAAGGAAGATACACTGGTTGAAGTACGTTCTACtt                                                                                 >  2:140322/1‑133 (MQ=255)                                                                              ccAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCt                                                  >  1:441264/1‑133 (MQ=255)                                                                                    tCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCg                                            <  2:414290/133‑1 (MQ=255)                                                                                                       aGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGAGTCGACGCATTGAATAAGCGTCg                         >  1:217318/1‑133 (MQ=255)                                                                                                       aGCGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCg                         <  2:124318/133‑1 (MQ=255)                                                                                                         cGGAAGATCCACAGAAACTGGTTCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGCTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGgc                       >  1:513449/1‑133 (MQ=255)                                                                                                             aGATCCACAGAAACTAGTTCGCCTGATGTTCCAGGAGATGGTAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGt                   <  2:161889/133‑1 (MQ=255)                                                                                                                              ttCGCCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGa  >  1:109296/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                                 GCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTGATGATCCAGGAGATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGA  >  NC_000913/1368039‑1368296

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