Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,515,933 C→A 8.1% P276T (CCG→ACG)  ftsI → peptidoglycan DD transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,515,9330CA 8.1% 150.8 / 2.6 49P276T (CCG→ACG) ftsIpeptidoglycan DD transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base A (2/2);  total (25/24)
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.44e-01

CTGGCTACCCGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGCCTG  >  NC_002947/1515789‑1516070
                                                                                                                                                |                                                                                                                                          
cTGGCTACCCGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCcggcg                                                                                                                                       >  1:526162/1‑150 (MQ=255)
     tACCCGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCAt                                                                                                                                  <  1:14383/150‑1 (MQ=255)
       cccGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCTGATCTGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGc                                                                                                                                <  1:437975/150‑1 (MQ=255)
       cccGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGc                                                                                                                                <  1:67307/150‑1 (MQ=255)
        ccGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCg                                                                                                                               >  2:354254/1‑150 (MQ=255)
        ccGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCg                                                                                                                               >  1:210011/1‑150 (MQ=255)
          gcgAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAAACGACTGACCAGGTTCCCCGCGGCCATGCGTa                                                                                                                             >  2:323444/1‑150 (MQ=255)
                 gcgcAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGc                                                                                                                      <  1:354254/150‑1 (MQ=255)
                     aaCGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGc                                                                                                                                >  2:349640/1‑136 (MQ=255)
                     aaCGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGc                                                                                                                                <  1:349640/136‑1 (MQ=255)
                         ccATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCg                                                                                                              >  1:848854/1‑150 (MQ=255)
                           aTTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGAc                                                                                                            <  2:436367/150‑1 (MQ=255)
                                 gagcAGGGCGCCAAGGGCGGCAGGCGGGTGATTATGGATGTCAAGCGCGGTGAGGTACTGGCGATGGTCAAGCAAGCGACCTACAACCCCCACGACCGACGTAGCATGTTCCCGGCGGCCCTGCGTAACCGGGCGATCATCGACGTGTTc                                                                                                      <  2:137319/147‑1 (MQ=255)
                                 gAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTc                                                                                                      >  1:333348/1‑150 (MQ=255)
                                  aaCATGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCggg                                                                                                                       <  2:112646/132‑1 (MQ=255)
                                  aaCATGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCggg                                                                                                                       >  1:112646/1‑132 (MQ=255)
                                     aGGACGCCAAGGCAAGCAGTCCGAAGCTAATAAATGTAAAGAACGGGTATGTCCTGGCCATTTTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGc                                                                                                  >  2:650313/1‑150 (MQ=255)
                                         cGCCAATGCCGGCAGCCTGTTGATTATGGATGTCAAGACCGGTGAGGTCCTGGCCAGGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCggg                                                                                             <  2:280931/150‑1 (MQ=255)
                                                          tGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCATCAGCCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGTTCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGa                                                                            <  1:323444/150‑1 (MQ=255)
                                                                  aTGGATGTCAAGACCGGTGGGGTCCTGGCGATGGTCAACCAACCGACCTACAACCCCAACAACCGACATAGCATGTTACCGGCGGCCATGCGTAACCGAGCGATCATCGACGTGTTCGAGCCGTGCTGCACGGTCAAGCCGATCTCCATg                                                                    >  1:298717/1‑150 (MQ=255)
                                                                  aTGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATg                                                                    >  2:737675/1‑150 (MQ=255)
                                                                   tGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGa                                                                   <  1:632853/150‑1 (MQ=255)
                                                                   tGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGa                                                                   >  1:406967/1‑150 (MQ=255)
                                                                   tGGATCTCAAGCCCGCTCAAGTCCTGGCAAAGGTGCAACAAGCGAGCTATATCCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGa                                                                   <  2:225429/150‑1 (MQ=255)
                                                                           aaGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCACCGGTGGCTCCACGGTCAAGCCGATCTCCATGAGTgcggcg                                                           >  1:437637/1‑150 (MQ=255)
                                                                                  gTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGc                                                                                                                                <  1:480876/75‑1 (MQ=255)
                                                                                  gTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGc                                                                                                                                >  2:480876/1‑75 (MQ=255)
                                                                                   tGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCagag                                                   >  2:305275/1‑150 (MQ=255)
                                                                                   tGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTAGAGCCGGACTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCagca                                                   >  1:91459/1‑148 (MQ=255)
                                                                                                   gTCAACCAACCGATCTATAACCCCCACCACCGCGGTAGTATGTTCCCGGCGGCCACGCGTAACCGGGCGATCATGCACGTCGGCGAGCCGGGCTCTGGGGTCCACCCGTTCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGccg                                   <  2:210011/150‑1 (MQ=255)
                                                                                                      aaccaaccGACCTACAACCCCAACAACCGACGTAGCATGTTCACGGCGGCCATGc                                                                                                                                <  1:471878/55‑1 (MQ=255)
                                                                                                      aaccaaccGACCTACAACCCCAACAACCGACGTAGCATGTTCACGGCGGCCATGc                                                                                                                                >  2:471690/1‑55 (MQ=255)
                                                                                                      aaccaaccGACCTACAACCCCAACAACCGACGTAGCATGTTCACGGCGGCCATGc                                                                                                                                >  2:471878/1‑55 (MQ=255)
                                                                                                      aaccaaccGACCTACAACCCCAACAACCGACGTAGCATGTTCACGGCGGCCATGc                                                                                                                                <  1:471690/55‑1 (MQ=255)
                                                                                                       accaaccGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGccgaccg                               >  1:340538/1‑150 (MQ=255)
                                                                                                        ccaaccGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGccgaccga                              >  1:833046/1‑150 (MQ=255)
                                                                                                           accGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACaa                           >  1:283075/1‑150 (MQ=255)
                                                                                                              gACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGt                        >  1:292226/1‑150 (MQ=255)
                                                                                                                 cTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGa                     <  2:283075/150‑1 (MQ=255)
                                                                                                                  tACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAg                    <  2:848854/150‑1 (MQ=255)
                                                                                                                    cAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGgt                  <  2:406967/150‑1 (MQ=255)
                                                                                                                             caacCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGGGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTgg         <  2:333348/150‑1 (MQ=255)
                                                                                                                              aacCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGaa                                       <  2:257623/119‑1 (MQ=255)
                                                                                                                              aacCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGaa                                       >  1:257623/1‑119 (MQ=255)
                                                                                                                                 cGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGc     >  1:547207/1‑150 (MQ=255)
                                                                                                                                 cGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGc     <  1:305275/150‑1 (MQ=255)
                                                                                                                                 cGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGc     >  1:139217/1‑150 (MQ=255)
                                                                                                                                 cGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTAGCGCGCGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGc     <  1:650313/150‑1 (MQ=255)
                                                                                                                                    cGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGCCTg  <  1:737675/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                          
CTGGCTACCCGCGAGTTGCGCAACGCCATTGCCGAACAGGACGCCAAGGCCGGCAGCCTGGTGATCATGGATGTCAAGACCGGTGAGGTCCTGGCCATGGTCAACCAACCGACCTACAACCCCAACAACCGACGTAGCATGTTCCCGGCGGCCATGCGTAACCGGGCGATCATCGACGTGTTCGAGCCGGGCTCCACGGTCAAGCCGATCTCCATGAGTGCGGCGCTGCAGAGTGGCCGCTGGAAGCCGACCGACAAGGTCGAGGTGTACCCTGGCAGCCTG  >  NC_002947/1515789‑1516070

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