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%	162.7 / NA	52	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/31);  total (21/31)

CTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGG  >  NC_000913/451506‑451779                                                                                                                                          |                                                                                                                                         cTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTaa                                                                                                                                     >  1:88977/1‑143 (MQ=255)    tCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTACCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCACTTCCACGGGATCTGGTGGCATTTAAGgt                                                                                                                                  >  2:31468/1‑143 (MQ=255)      cTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAgtgt                                                                                                                                >  1:647107/1‑143 (MQ=255)           ggATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTggg                                                                                                                           >  1:52076/1‑143 (MQ=255)            gATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGGCAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGc                                                                                                                          <  2:598949/143‑1 (MQ=255)            gATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCCAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGc                                                                                                                          >  1:478501/1‑143 (MQ=255)              tCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATGTGGTGGCATTTAAAGTGTGTGGTCAt                                                                                                                        >  2:290960/1‑143 (MQ=255)                 aaCAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTATTATCGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTCAAAGTGTGTGGGCAtttt                                                                                                                     >  1:202335/1‑143 (MQ=255)                     gcgcAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACgg                                                                                                                 <  2:457544/143‑1 (MQ=255)                        cAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGaa                                                                                                              <  1:649671/143‑1 (MQ=255)                         aGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAg                                                                                                             >  2:157207/1‑143 (MQ=255)                         aGAATTACAGACACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAg                                                                                                             <  2:597042/143‑1 (MQ=255)                              tACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTAc                                                                                                        >  2:613484/1‑143 (MQ=255)                                  gCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTa                                                                                                    <  2:647107/143‑1 (MQ=255)                                   ccACTTAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTAc                                                                                                   <  1:307247/143‑1 (MQ=255)                                       tGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTa                                                                                               >  2:344881/1‑143 (MQ=255)                                         aGCAAGACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTaaa                                                                                             <  2:118071/143‑1 (MQ=255)                                              tACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGa                                                                                        <  2:430403/143‑1 (MQ=255)                                               aCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGAt                                                                                       >  2:503149/1‑143 (MQ=255)                                                        gcgaatcatGACTACCATCCCAAACTTTTGTTGTATTTCCTGCGTCTGATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCt                                                                              <  1:6970/140‑1 (MQ=255)                                                        gCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCt                                                                              <  2:327591/143‑1 (MQ=255)                                                        gCAAAGAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCt                                                                              <  2:190842/143‑1 (MQ=255)                                                               atGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACtt                                                                       <  2:715366/143‑1 (MQ=255)                                                                 gACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCGCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTGCGGGGAGGTCGCAGTACTGTAAACATGATTAAATTCCTGTTACTTgg                                                                     >  1:737639/1‑143 (MQ=255)                                                                  aCAACCGTCCCAATCTTTTATTGTATTTCCTGAGTCTCATATAACGACCTCAAGTCCACGGGATCGGGTGGCATTTAAAGTGTGTGGTCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGa                                                                    <  1:285849/143‑1 (MQ=255)                                                                  aCAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCGTGTTACTTGGa                                                                    >  2:83009/1‑143 (MQ=255)                                                                   cAACCATCCCACCCTCTTATTGTATTTCCTGAGTCCCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGAt                                                                   <  2:86820/143‑1 (MQ=255)                                                                    aaCCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGTTGGCATTTAAAGTGTGTGGGAATTTTACGGGAAGGTTACATCACTGTAAACATGATTAAATTCCTGTTACTTgagtt                                                                  >  2:547680/1‑139 (MQ=255)                                                                      ccATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTg                                                                <  1:503149/143‑1 (MQ=255)                                                                         tCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACGTTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCt                                                             >  1:545150/1‑143 (MQ=255)                                                                           ccAGACCTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGc                                                           <  2:208891/143‑1 (MQ=255)                                                                              aaCTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAg                                                        <  2:733876/143‑1 (MQ=255)                                                                                 ttttATTGTATTTCCTGAGTCTCATGTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTc                                                     <  1:15564/143‑1 (MQ=255)                                                                                        gTATTTCCTGAGTCTCATTTAACGACATCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGAATCTTAAGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATg                                              >  1:228267/1‑143 (MQ=255)                                                                                          agtgacggAGTCTAATTTAACGACCTCAATGTCACGGGATCGGGTGGCATTTAAAGTGTGTTGGCAGTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGtt                                            <  2:30197/136‑1 (MQ=255)                                                                                          aTTGCCTGAGTCTCATTTAACGACTGCACTTCATGTCGTGTTTTTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGtt                                            <  1:737566/143‑1 (MQ=255)                                                                                                tGAGTCTCATTTAGCGACCTCAATTCCACGGGATCTGGTGGCATGTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCaaaa                                      <  2:582385/143‑1 (MQ=255)                                                                                                tGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCaaaa                                      <  2:266180/143‑1 (MQ=255)                                                                                                tGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCaaaa                                      <  1:241292/143‑1 (MQ=255)                                                                                                   gTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACat                                   <  2:279421/143‑1 (MQ=255)                                                                                                         tttAACGACCTCAATTTCACGGGATCTGGTGGCATTTAAAGTGTTTGCGTATTTTTCGGGAAGGTTACATTACTGTAAACATGATTAAATTGCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAgg                             <  1:175305/143‑1 (MQ=255)                                                                                                           tAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTg                                           <  2:128918/127‑1 (MQ=255)                                                                                                           tAACGACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTg                                           >  1:128918/1‑127 (MQ=255)                                                                                                               gACCTCAATTCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTGGAtt                       >  2:420942/1‑143 (MQ=255)                                                                                                                 ccTCAATTCCACGGGAGCTGGCGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACAGGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGAtttt                     <  2:142642/143‑1 (MQ=255)                                                                                                                   ttaaTTCCACGGGCTGTGGTGGCATGTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATGCCTGCTACTTGGCTTTGGCTGCGAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTaa                   <  2:144774/141‑1 (MQ=255)                                                                                                                        tCCACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGc              <  2:239329/143‑1 (MQ=255)                                                                                                                          cACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTACATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTg            >  1:50799/1‑143 (MQ=255)                                                                                                                          cACGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTg            >  2:262384/1‑143 (MQ=255)                                                                                                                            cGGGATCTGGTGGCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGtt          <  2:484284/143‑1 (MQ=255)                                                                                                                                 tgtggtggTATTTAAAGTGTGTGGCCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATGTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCaaa     <  1:360586/141‑1 (MQ=255)                                                                                                                                    ggtggCATTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAgg  >  2:442849/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                         CTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGG  >  NC_000913/451506‑451779

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: ‑