Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 395,831 G→A R153C (CGT→TGT)  ampH ← peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913395,8310GA95.5% 52.8 / ‑3.7 22R153C (CGT→TGT) ampHpeptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (10/11);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

GCGGCACTGATCCGGTGAGGGGGTGTAGGTGGTGTCTTTCATCCCCAATGGACGGGTAATCTGCTCTTCAAACAACTGGGTATAAGGCTTGCCGGACGCATTCGCCAGTGCATCTGCCAGCAGGTCGAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACGATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTTTAACTTCACGGTCCCCTGGTCGAGCAATTTCACCAGCATCTCACTGGTCATTAGCTTGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGC  >  NC_000913/395598‑396061
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
gcggcACTGATCCGGTGAGGGGGTGTAGGTGGTGTCTTTCATCCCCAATGGACGGGTAATCTGCTCTTCAAACAACTGGGTATAAGGCTTGCCGGACGCATTCGCCAGTGCATCTGCCAGCAGGTCGAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCgg                                                                                                                                                                                                                       >  1:179575/1‑251 (MQ=255)
                                                       gTAATNNGCNNNTANNACAACTGGGAATAAGGCTTGCCGNACGCATTCGCCAGTGCATCTGACAGAAGGNCNANCGCANNGTTNGCGTACGCCGNCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCANCNTTGCTAGCGCGTTGGCCAGACAAACACCGGANAATGTGCcg                                                                                                                                                                                                                                >  2:369758/1‑187 (MQ=255)
                                                       gTAATCTGCTCTTCAAACAACTGGGTATAAGGCTTGCCGGACGCAGTCGCCATTGCATCTGCCAGCAGTTCTAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTTGCGCGTCTTGCAGCTTCGCCGTAGAAAGGTACTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCcg                                                                                                                                                                                                                                <  1:369758/187‑1 (MQ=255)
                                                                                          gCCGGACGCATTCGCCAGTGCATCTGCCAGCAGGTCGAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAggg                                                                                                                                                                                                     <  1:29485/179‑1 (MQ=255)
                                                                                                  cATTCGCCAGTGCATCTGCCAGCAGGTCGAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCatat                                                                                                                     >  1:337238/1‑251 (MQ=255)
                                                                                                                   gCCAGCAGGTCGAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGGCACCGGGCTGTTCACGGGGCAGGGCGCGGGTATTGGTTGCCAGATTGGCCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATAGTTACGTAGCGGATCGtt                                                                                                    <  1:411238/251‑1 (MQ=255)
                                                                                                                    ccAGCAGGTCGAACGCCAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGtt                                                                                                                                                >  1:158429/1‑206 (MQ=255)
                                                                                                                                   cAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGc                                                                                                                                                                        <  1:463221/167‑1 (MQ=255)
                                                                                                                                   cAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGACGCGCCACCGGGCTGTTCACgggg                                                                                                                                                                                                          >  1:78450/1‑133 (MQ=255)
                                                                                                                                               aCGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTGTTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGGTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGc                                                                                                                              <  1:508096/197‑1 (MQ=255)
                                                                                                                                                     ccTGACAGCCTGGCGCGTCTTTCAGCTTCTCCGGAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACTCCGGACAATGTGCCGCGCCACCGGGCTGGTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGAcc                                                                                                                                                                           <  1:334818/146‑1 (MQ=255)
                                                                                                                                                      cTGACNNCCNNNCGNNGCTTTCAGCTTCGCCGTAGAAAGNTATTTCCAGCGTTGCTCGCGCGTTGGCCANANANACACNNGACNATGTGCCGCGNCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGANCNGCGTAATGGGCGTCCCATTGTAGGTTGGCACNCGTGCGCCTGGCGGAGCATATTTANTTAGNNGANcnnnnnnnn                                                                                               >  2:366253/1‑213 (MQ=255)
                                                                                                                                                                      gCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACtt                                                                                                               <  1:231874/189‑1 (MQ=255)
                                                                                                                                                                      gCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGACAGACAAACACCGGACACTGTGCCGCGCCAACggg                                                                                                                                                                                                                      >  1:245257/1‑86 (MQ=255)
                                                                                                                                                                                  gCCGGAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGAGTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGcn                                                                                                           <  1:55610/181‑2 (MQ=255)
                                                                                                                                                                                    cGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCTCGAAAAACACCGAACACTGTGCCGCGCACCCGGGCTGTCCACGGGGCAGGGCGCTGGTATGGATTGCCAGACTGACCAg                                                                                                                                                                         >  1:249332/1‑117 (MQ=255)
                                                                                                                                                                                     gTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGACTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTTTAACTTCACGGTCCCCTGGTCGAGCAATTTCACCAGCATCTCACTGGt                                                    <  1:479767/233‑1 (MQ=255)
                                                                                                                                                                                                    cagcgtgtcccgcgcgTTGGCCAAACAAACCCTGGACAATGGGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAAGGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTTTAACTTCACGGCCCACTGGGCGAGCAATTTCACCAGCATCTCACTGGTCATTAGCTTGGGGAGGGAAGCAATACGCACGAc                   <  1:366253/241‑1 (MQ=255)
                                                                                                                                                                                                    cAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTCGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTTTAACTTCACGGTCCCCTGGTCGAGCAATTTCACCAGCATCTCACTGGTCATTAGCTTGGTGAGGGAAGCAATACGCACGAc                   >  1:205818/1‑251 (MQ=255)
                                                                                                                                                                                                           gCTCGCGCGTTGGCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTNNCCAGCg                                                                                                                                                                       >  2:8836/1‑96 (MQ=255)
                                                                                                                                                                                                                       gCCAGACAAACACCGGACAATGTGCCNCGCCACNNGGNTGTTCACGGGGCAGGGCGCNGGNNNGNNNTNCCAGANTGNNCAGCGTAATNGGCGTCCCATTGTAGGTTgg                                                                                                                                              >  2:110072/1‑109 (MQ=255)
                                                                                                                                                                                                                       gCCAGACAAACACCGGACAATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTgg                                                                                                                                              <  1:110072/109‑1 (MQ=255)
                                                                                                                                                                                                                                        cGATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGGTGGCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACTCGTGCGCCTGGCTGAGCATATTTACTTAGCGGATCGTTTATCTTCACGGTCCGCTGGGCGAGCAATTTCACCAGCATCTCACTGTGCATTAGCTTGTTGAGGGAAGCAAGACGCACGACCGAATCCAGCTGCGGGc  <  1:350302/232‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
GCGGCACTGATCCGGTGAGGGGGTGTAGGTGGTGTCTTTCATCCCCAATGGACGGGTAATCTGCTCTTCAAACAACTGGGTATAAGGCTTGCCGGACGCATTCGCCAGTGCATCTGCCAGCAGGTCGAACGCAAGGTTAGAGTACGCCGCCTGACTGCCTGGCGCGGCTTTCAGCTTCGCCGTAGAAAGGTATTTCCAGCGTTGCTCGCGCGTTGGCCAGACAAACACCGGACGATGTGCCGCGCCACCGGGCTGTTCACGGGGCAGGGCGCTGGTATGGGTTGCCAGATTGACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTTTAACTTCACGGTCCCCTGGTCGAGCAATTTCACCAGCATCTCACTGGTCATTAGCTTGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGC  >  NC_000913/395598‑396061

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