Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,328,842	T→C	intergenic (‑127/‑121)	greA ← / → dacB	transcription elongation factor GreA/peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,328,842	0	T	C	100.0%	149.4 / NA	43	intergenic (‑127/‑121)	greA/dacB	transcription elongation factor GreA/peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/23);  total (20/23)

CTTGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTGCGCCGTAGTATGACGGCTCGATTCCAGGTTG  >  NC_000913/3328727‑3328949                                                                                                                    |                                                                                                            ctTGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACAtt                                                                                            >  2:418746/1‑133 (MQ=255) ctTGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGCGTAAAAGGCCCAGGCCTGCTAGAATTCACAtt                                                                                            >  2:133384/1‑133 (MQ=255)  ttGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTTCAAGTCCCAGGCCATCTACAATTCACAttt                                                                                           >  1:202077/1‑133 (MQ=255)  ttGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACAttt                                                                                           >  1:1304/1‑133 (MQ=255)   tGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACAtttt                                                                                          >  2:306729/1‑133 (MQ=255)   tGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACAtttt                                                                                          >  2:6452/1‑133 (MQ=255)        atTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAg                                                                                     <  1:418746/133‑1 (MQ=255)          tCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGTGTACAAGTGCCACGCCGGCTACAATGCACATTTTGATAGTc                                                                                   >  2:158988/1‑133 (MQ=255)            cTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCAt                                                                                 >  1:431312/1‑133 (MQ=255)                 aGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTAc                                                                            >  2:45114/1‑133 (MQ=255)                  gggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTCCAAGTCACATTTTGATAGTCATTTTAcc                                                                           >  2:179563/1‑133 (MQ=255)                  gggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTAcc                                                                           <  2:376057/133‑1 (MQ=255)                   ggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTAccc                                                                          >  2:315317/1‑133 (MQ=255)                   ggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTAccc                                                                          >  2:260114/1‑133 (MQ=255)                   ggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTAccc                                                                          >  1:475947/1‑133 (MQ=255)                   ggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCACGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTTCCAGGCCAGCTAAAATTCACATTTTTATAGTCATTTTAccc                                                                          >  2:194753/1‑133 (MQ=255)                   ggCAAGTCCTCACCCGCCCCGGCTGTTCGGTCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTAccc                                                                          <  1:131585/133‑1 (MQ=255)                              acccacccCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCg                                                               >  1:389422/1‑133 (MQ=255)                              acccacccCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCg                                                               >  1:367588/1‑133 (MQ=255)                              acccacccCGGCTGTGCGCCCTGCACGGCAAAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCg                                                               <  1:194753/133‑1 (MQ=255)                                  acccCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAgg                                                           >  2:438380/1‑133 (MQ=255)                                   ccccGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAggg                                                          >  2:397760/1‑133 (MQ=255)                                   ccccGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAggg                                                          <  2:325300/133‑1 (MQ=255)                                   ccccGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAggg                                                          <  1:293866/133‑1 (MQ=255)                                   ccccGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGAACAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAggg                                                          <  2:183534/133‑1 (MQ=255)                                     ccGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTc                                                        <  2:389422/133‑1 (MQ=255)                                     ccGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTc                                                        <  2:367588/133‑1 (MQ=255)                                      cGGCTGTTCGGCCTGCCCGGCATATCTCGTTAACGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCa                                                       <  2:202077/133‑1 (MQ=255)                                      cGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCa                                                       <  1:438380/133‑1 (MQ=255)                                        gCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATc                                                     <  1:62745/133‑1 (MQ=255)                                          tGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGt                                                   <  1:258874/133‑1 (MQ=255)                                          tGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGt                                                   <  1:6452/133‑1 (MQ=255)                                          tGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGt                                                   <  1:397760/133‑1 (MQ=255)                                                      ccGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAgg                                       >  1:94624/1‑133 (MQ=255)                                                         gCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCg                                    <  1:306729/133‑1 (MQ=255)                                                            tAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTg                                 <  2:431312/133‑1 (MQ=255)                                                            tAACTCCTCACCGGGCCCCGAATCAAAAAACTACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTg                                 <  1:133384/133‑1 (MQ=255)                                                              aCTCCTTAACGGGACTAGTAGCTAAAAAATCCCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTgcg                               <  1:158988/133‑1 (MQ=255)                                                               cTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTgcgc                              <  1:260114/133‑1 (MQ=255)                                                               cTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTgcgc                              >  1:359305/1‑133 (MQ=255)                                                                 ccTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTGCGCCg                            <  1:67193/133‑1 (MQ=255)                                                                    tACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTGCGCCgtag                         >  2:296063/1‑133 (MQ=255)                                                                                    aaaaaaaTACCGACCGGGTTACCAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTGCGCCGTAGTATGACGGCTCGATTc         <  1:179563/133‑1 (MQ=255)                                                                                           tACCGACCCGGGTACAAGTCCCAGGCCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTGCGCCGTAGTATGACGGCTCGATTCCAGgttg  <  1:315317/133‑1 (MQ=255)                                                                                                                    |                                                                                                            CTTGAATATTCCTGATAGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTTACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCATTTTACCCTGAAGTTCCCGAAGGGTCATCGTTTACTTTATAGGGCGTTGCGCCGTAGTATGACGGCTCGATTCCAGGTTG  >  NC_000913/3328727‑3328949

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