Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	451,643	G→A	intergenic (‑33/+427)	cyoA ← / ← ampG	cytochrome o ubiquinol oxidase subunit II/muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	451,643	0	G	A	100.0%	153.8 / NA	50	intergenic (‑33/+427)	cyoA/ampG	cytochrome o ubiquinol oxidase subunit II/muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/29);  total (21/29)

ATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAG  >  NC_000913/451504‑451778                                                                                                                                            |                                                                                                                                        aTCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCAttt                                                                                                                                      <  1:476986/143‑1 (MQ=255)    tGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTaaa                                                                                                                                   <  1:403998/143‑1 (MQ=255)    cgTCCGTTGGGGTCTAGCAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACTACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTaaa                                                                                                                                   <  2:657558/142‑1 (MQ=255)         tttGGGATCTAACAGCGCAGAATTACAGCCACTGAGAAATACAGTGGCTGCAAAAAATGACAACCATCCCAAACTTTTCTTGTATATCCTGAGTCTCATTTAACGACCTCAATTCCACGGCACCTGGTGGCATTTAAAgtgtg                                                                                                                              >  1:369262/1‑143 (MQ=255)            tggATCTAACATCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAAATTTTTTTGTAGTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTgg                                                                                                                           <  1:310147/142‑1 (MQ=255)              gATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCCATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGc                                                                                                                         >  2:715784/1‑143 (MQ=255)                tCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCAt                                                                                                                       >  1:528873/1‑143 (MQ=255)                      aGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACg                                                                                                                 <  1:293721/143‑1 (MQ=255)                      aGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACg                                                                                                                 <  1:293723/143‑1 (MQ=255)                      aGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCAGTTTACg                                                                                                                 >  2:252320/1‑143 (MQ=255)                       ccgcAGAATTACAGCCACTGAGCAATACAGTGCCCGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACgg                                                                                                                <  1:258807/142‑1 (MQ=255)                           aGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAg                                                                                                            <  1:584563/143‑1 (MQ=255)                              agtACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGtt                                                                                                         <  2:783085/141‑1 (MQ=255)                               ttACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTa                                                                                                        <  1:273959/143‑1 (MQ=255)                                      cACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACt                                                                                                 <  1:682149/143‑1 (MQ=255)                                       aCTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGa                                                                                                              >  2:156130/1‑129 (MQ=255)                                       aCTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGa                                                                                                              <  1:156130/129‑1 (MQ=255)                                       aCTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTg                                                                                                <  1:794565/143‑1 (MQ=255)                                         tGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTa                                                                                              >  2:753909/1‑143 (MQ=255)                                                     tGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTaaa                                                                                  <  1:242746/143‑1 (MQ=255)                                                     tGCCTGCAAATAACGACAACCATCCCAAACTTTTATTTTATTTCCTTATTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTaaa                                                                                  <  1:257806/143‑1 (MQ=255)                                                         tGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTcc                                                                              >  1:727118/1‑143 (MQ=255)                                                          gCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCt                                                                             <  2:574257/143‑1 (MQ=255)                                                           cagagaaggagacccATCCCAAACTGTGAGTGGAGTGCCGGTGTCGCATTTAGCGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATGTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTg                                                                            <  1:588817/130‑1 (MQ=255)                                                                    aCAACCCTCCGAAACTTTTATTGGTTTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGa                                                                   <  2:82655/143‑1 (MQ=255)                                                                       aCCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGAATCTGGTGGCATTTAAAGTGTGTGGGCATTATACCGGAAGGTTACATTACTGTCAACATGATTAAATTCCTGTTCCTTGGAtct                                                                >  1:615137/1‑141 (MQ=255)                                                                            cccAAACTTTTGTTGTATTTCCTGACTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTg                                                           <  1:698785/143‑1 (MQ=255)                                                                              cAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTTTGAGCATTTTACGGGAAGGATACATTACTGTAAACATGATTAAATTACTATTACGTGAATTTGGCTggc                                                         >  2:21031/1‑141 (MQ=255)                                                                                  ctttgcttGTGTTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGGGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGt                                                     <  2:185519/137‑1 (MQ=255)                                                                                   ttttATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGATACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGcctgttc                                                    >  1:339672/1‑138 (MQ=255)                                                                                    tttATTGTGTTTCCTGAGTCTCATTTAACGAGCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTca                                                   <  2:703832/143‑1 (MQ=255)                                                                                          gTATTTCCTGGGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATg                                             <  2:355484/143‑1 (MQ=255)                                                                                                  tGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTCTAAACATGCTTACATTTGAGCTGCGTGGATTTGGCTGCCAGGTCACATATGTTGCaaaa                                     >  2:802810/1‑143 (MQ=255)                                                                                                   tggTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAAc                                    <  1:552001/141‑1 (MQ=255)                                                                                                       ctcATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACatat                                <  1:767576/143‑1 (MQ=255)                                                                                                        tcATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATc                               >  1:72420/1‑143 (MQ=255)                                                                                                           tttAACGACCTCAATTCCACGGGATCTGGTGGCATTTATAGAGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCGCATACGTTGCAAAACATATCAgg                            >  2:453351/1‑143 (MQ=255)                                                                                                                 gACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGAtt                      <  1:795494/143‑1 (MQ=255)                                                                                                                     tCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTTGCTGCCAGGTCACATATGTGGCAAAACATATCAGGTTGATTGTaa                  >  1:353258/1‑143 (MQ=255)                                                                                                                       aaTTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCGTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATCTCAGGTTGCTTTTaaga                >  2:497894/1‑143 (MQ=255)                                                                                                                         ttCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCCCCCATGCACCAAAACATATCAGGTTGATTTTaagaag              >  1:780129/1‑143 (MQ=255)                                                                                                                          tCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGGTTTTAAGAAGc             >  1:749055/1‑143 (MQ=255)                                                                                                                             ccGGGATTTTGTGGCATTTAAAGTCTGTGTGCATTTTACGTGAAGGTTACATTACTGTAAGCATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGt          <  2:471285/142‑1 (MQ=255)                                                                                                                             aCGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTACTGTTACATGGATTTGGCTGCCAGGTCACATATATTGCAAAACATATCAGGTTGATTTTAAGCAGCTGt          >  1:453090/1‑143 (MQ=255)                                                                                                                              cGGGATCGGGTTGCATTTAAAGTGTGTGGGCACTTTCCGTGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGATGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGtt         <  2:365195/143‑1 (MQ=255)                                                                                                                               cgcaTCTGGTGCCAGTCAAAGCGTGTGGGCATTTTACCGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTg        <  1:499655/140‑1 (MQ=255)                                                                                                                                 gATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCa      >  1:355991/1‑143 (MQ=255)                                                                                                                                  aTCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGACGCTGTTGCaa     >  2:115456/1‑143 (MQ=255)                                                                                                                                    cTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCaaaa   <  1:385583/143‑1 (MQ=255)                                                                                                                                     tggtggCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAg  >  2:480070/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                        ATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAG  >  NC_000913/451504‑451778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑