Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,518,500 G→A 13.8% L55L (CTG→CTA murF → D‑alanyl‑D‑alanine‑adding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,518,5000GA13.8% 62.1 / 3.0 29L55L (CTG→CTAmurFD‑alanyl‑D‑alanine‑adding enzyme
Reads supporting (aligned to +/- strand):  ref base G (13/12);  new base A (2/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

ACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCAGCTCGGTGCATTGAACCGCGCCGCTTTCGACAAGCCAGTG  >  NC_002947/1518353‑1518647
                                                                                                                                                   |                                                                                                                                                   
aCTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGc                                                                                                                                                   >  2:268791/1‑150 (MQ=255)
    agccagcTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGAt                                                                                                                                               <  1:268791/150‑1 (MQ=255)
    agccagcTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTAGCAGAt                                                                                                                                               >  2:179501/1‑150 (MQ=255)
    acccagcTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTAGCAGAt                                                                                                                                               >  2:179390/3‑150 (MQ=255)
                      cTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGcgcg                                                                                                                             >  1:507834/1‑150 (MQ=255)
                        gAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGt                                                                                                                           <  1:445866/150‑1 (MQ=255)
                             gTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAgccg                                                                                                                      >  2:560315/1‑150 (MQ=255)
                                aaCACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGccc                                                                                                                   <  1:173456/150‑1 (MQ=255)
                                     ataggTGGCGCCGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGGGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTc                                                                                                              >  2:134070/4‑150 (MQ=255)
                                       cGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGa                                                                                                            <  2:507834/150‑1 (MQ=255)
                                              gACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTg                                                                                                       >  1:390116/1‑148 (MQ=255)
                                              gACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTg                                                                                                       <  2:390116/148‑1 (MQ=255)
                                                                          aTAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGc                                                                         >  1:34630/1‑150 (MQ=255)
                                                                                 cAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTAGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGc                                                                  <  1:179390/150‑1 (MQ=255)
                                                                                 cAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTAGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGc                                                                  <  1:179501/150‑1 (MQ=255)
                                                                                      tCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTg                                                                                                       >  1:479390/1‑108 (MQ=255)
                                                                                      tCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTg                                                                                                       <  2:479390/108‑1 (MQ=255)
                                                                                      tCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGAtt                                                             <  2:3477/150‑1 (MQ=255)
                                                                                               gtggaCTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTgg                                                    >  1:21554/5‑150 (MQ=255)
                                                                                               ggCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCTGTTGTAGGTGGCCGATTGCCGCCTgg                                                    <  1:124828/150‑1 (MQ=255)
                                                                                                               ccTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCa                                          <  2:568221/144‑1 (MQ=255)
                                                                                                               ccTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCa                                          >  1:568221/1‑144 (MQ=255)
                                                                                                                 tGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCTAGCTCGGTg                                  <  2:2527/150‑1 (MQ=255)
                                                                                                                              gttCGATGGCCATGACTACCTGGCAGATGTGAGAGCCAAAGGGGCGGTAGCAGCACTGGTGGAGCGTGACGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCTCTCGGCACGCTAGGTGCATTGAACGgcgc                     >  2:59149/2‑150 (MQ=255)
                                                                                                                                            aCTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCAGCTCGGTGCATTGAACCGCGCCGCTTTCGACAAGc       <  1:59149/150‑1 (MQ=255)
                                                                                                                                             cTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGcc                                                  >  2:133594/1‑106 (MQ=255)
                                                                                                                                             cTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGcc                                                  <  1:133594/106‑1 (MQ=255)
                                                                                                                                              tACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCAGCTCGGTGCATTGAACCGCGCCGCTTTCGACAAGCCa     >  1:536780/1‑150 (MQ=255)
                                                                                                                                                 cTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCAGCTCGGTGCATTGAACCGCGCCGCTTTCGACAAGCCAgtg  >  2:364007/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCAGCATCGATAGCCGCAGTGTCAGCGCCGGGCAACTTTTCGTCGCCCTGGCCGGCCCGCGTTTCGATGGCCATGACTACCTGGCAGATGTGAAAGCCAAAGGCGCGGTAGCCGCCCTGGTCGAGCGTGAAGTAGCCGATGTCGACCTGCCGCAGTTGCAGGTGGCCGATTGCCGCCTGGCCCTCGGCCAGCTCGGTGCATTGAACCGCGCCGCTTTCGACAAGCCAGTG  >  NC_002947/1518353‑1518647

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