Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,382	G→T	A57E (GCG→GAG)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,867,382	0	G	T	100.0%	131.0 / NA	36	A57E (GCG→GAG)	rpoS	RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/21);  total (15/21)

CGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCT  >  NC_000913/2867305‑2867459                                                                              |                                                                              cgcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGt                                                                            >  2:517143/1‑81 (MQ=255)   cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTcc                                                                          <  2:517100/81‑1 (MQ=255)    gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCa                                                                         <  1:517101/81‑1 (MQ=255)     aCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCaa                                                                        <  1:517103/81‑1 (MQ=255)     aCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCaa                                                                        <  2:517102/81‑1 (MQ=255)               tAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCtgt                                                              <  1:517104/81‑1 (MQ=255)                 aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCtgtgt                                                            <  1:517105/81‑1 (MQ=255)                        tcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTccc                                                     >  1:517155/1‑81 (MQ=255)                         cttcGGCCGTTATCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTccc                                                     <  1:517106/80‑1 (MQ=255)                         cttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCt                                                    >  2:517160/1‑81 (MQ=255)                            cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCt                                                    <  1:517107/78‑1 (MQ=255)                             ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                >  1:517166/1‑81 (MQ=255)                               ccGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATa                                              <  1:517108/81‑1 (MQ=255)                                cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                             <  2:517109/81‑1 (MQ=255)                                 gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTTCAACACACGCTGTGTGGCTCCCTGCGATaa                                             >  2:517167/1‑80 (MQ=255)                                 gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                             >  2:517148/1‑80 (MQ=255)                                 gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                             >  1:517150/1‑80 (MQ=255)                                     aCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                        <  2:517110/81‑1 (MQ=255)                                       aGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTcc                                      >  1:517151/1‑81 (MQ=255)                                         tggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCctc                                    <  1:517111/81‑1 (MQ=255)                                            tgAATAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                 <  2:517112/81‑1 (MQ=255)                                              aaTAACCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCgg                               <  2:517113/81‑1 (MQ=255)                                                 aaCCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGcc                             >  1:517164/1‑80 (MQ=255)                                                  aCCAATCTCACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCa                            <  1:517114/80‑1 (MQ=255)                                                          cACCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                   >  2:517157/1‑81 (MQ=255)                                                           aCCAAGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                     <  2:517115/78‑1 (MQ=255)                                                           aCCAAGGTAAAGATGAGTCTCATCCAGCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                     <  2:517116/78‑1 (MQ=255)                                                               aGGTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTg              >  2:517161/1‑81 (MQ=255)                                                                 gTAAAGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTggg            >  1:517170/1‑81 (MQ=255)                                                                   aaaGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGtt          >  1:517162/1‑81 (MQ=255)                                                                    aaGCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTc         >  1:517174/1‑81 (MQ=255)                                                                      gCTGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCt       <  2:517117/81‑1 (MQ=255)                                                                        tGAGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGt     >  1:517156/1‑81 (MQ=255)                                                                          aGTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc   <  1:517118/81‑1 (MQ=255)                                                                           gTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTct  <  2:517119/81‑1 (MQ=255)                                                                           gTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTct  <  2:517120/81‑1 (MQ=255)                                                                              |                                                                              CGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCT  >  NC_000913/2867305‑2867459

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