Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	453,447	A→G	G33G (GGT→GGC)	ampG ←	muropeptide:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	453,447	0	A	G	100.0%	140.6 / NA	40	G33G (GGT→GGC)	ampG	muropeptide:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/20);  total (20/20)

TGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAG  >  NC_000913/453329‑453580                                                                                                                       |                                                                                                                                      tgtAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGt                                                                                                                  >  2:1069998/1‑140 (MQ=255) tgtAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGt                                                                                                                  >  1:739489/1‑140 (MQ=255)        gTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCgg                                                                                                           >  2:1064316/1‑140 (MQ=255)           cATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAg                                                                                                        >  2:856266/1‑140 (MQ=255)                   gTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAAccc                                                                                                >  2:801073/1‑140 (MQ=255)                   gTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAAccc                                                                                                >  1:675446/1‑140 (MQ=255)                       ccAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGa                                                                                            <  1:587047/140‑1 (MQ=255)                         agagGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATc                                                                                          <  2:1284/140‑1 (MQ=255)                           agGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAg                                                                                        <  2:675446/140‑1 (MQ=255)                           agGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCACCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAg                                                                                        <  2:133542/140‑1 (MQ=255)                                tttAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAgg                                                                                             <  2:983568/130‑1 (MQ=255)                                tttAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAgg                                                                                             >  1:983568/1‑130 (MQ=255)                                                   ccTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTg                                                                >  2:1010921/1‑140 (MQ=255)                                                   ccTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTg                                                                >  2:521470/1‑140 (MQ=255)                                                       ccAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCgg                                                                      >  2:971943/1‑130 (MQ=255)                                                       ccAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCgg                                                                      <  1:971943/130‑1 (MQ=255)                                                        cAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGaaaaa                                                           <  1:1010921/140‑1 (MQ=255)                                                          gagagaAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATa                                                         <  2:266417/140‑1 (MQ=255)                                                            gagaAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAAcgc                                                                        >  1:613988/1‑123 (MQ=255)                                                            gagaAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAAcgc                                                                        <  2:613988/123‑1 (MQ=255)                                                            gagaAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACg                                                       <  1:521470/140‑1 (MQ=255)                                                             agaagaAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGaaaaa                                                           <  1:620475/135‑1 (MQ=255)                                                             agaagaAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGaaaaa                                                           >  2:620475/1‑135 (MQ=255)                                                                    aCCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAat                                                  >  1:419660/1‑137 (MQ=255)                                                                    aCCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAat                                                  <  2:419660/137‑1 (MQ=255)                                                                    aCCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATAtt                                               >  2:883618/1‑140 (MQ=255)                                                                      cAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGa                                             <  1:856266/140‑1 (MQ=255)                                                                      cAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGa                                             <  2:438291/140‑1 (MQ=255)                                                                                  gaTCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGc                                 <  2:658408/140‑1 (MQ=255)                                                                                  gaTCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGc                                 <  1:801073/140‑1 (MQ=255)                                                                                               gACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCg                    >  2:566469/1‑140 (MQ=255)                                                                                               gACCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCg                    >  1:238732/1‑140 (MQ=255)                                                                                                aCCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCgg                   >  2:833998/1‑140 (MQ=255)                                                                                                aCCGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCgg                   >  1:448210/1‑140 (MQ=255)                                                                                                  cGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCt                 <  2:471739/140‑1 (MQ=255)                                                                                                  cGTCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCt                 <  1:1069998/140‑1 (MQ=255)                                                                                                    tCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGc                                 <  2:297923/122‑1 (MQ=255)                                                                                                    tCATCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGc                                 >  1:297923/1‑122 (MQ=255)                                                                                                      aTCCAGGCCTGTAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCa             <  1:883618/140‑1 (MQ=255)                                                                                                                 tAAGGTGCCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCagag  >  2:601391/1‑140 (MQ=255)                                                                                                                       |                                                                                                                                      TGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAG  >  NC_000913/453329‑453580

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