Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 981,097 C→A 7.1% C17* (TGC→TGA)  ldtD → murein L,D‑transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913981,0970CA 7.1% 194.5 / 2.2 56C17* (TGC→TGA) ldtDmurein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base A (2/2);  total (28/28)
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.06e-01

CATACTTTCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGAT  >  NC_000913/980959‑981234
                                                                                                                                          |                                                                                                                                         
cATACTTTCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGACTGGCCGt                                                                                                                                   >  2:310819/1‑147 (MQ=255)
cATACTTTCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGACTGGCCGt                                                                                                                                   >  2:150937/1‑147 (MQ=255)
   aCTTTCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACa                                                                                                                               >  1:99750/1‑148 (MQ=255)
      ttCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTAGTCAGCTGTCGGCAATCAGTTTTTGCCTGGCCGTAACATTc                                                                                                                            >  8:349836/1‑148 (MQ=255)
      ttCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGCGGCAAGGGATGTTGCTTAATATGAGGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTc                                                                                                                            >  4:235128/1‑148 (MQ=255)
            aGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGACTGGCCGTAACATTCGCTcc                                                                                                                       <  1:150937/147‑1 (MQ=255)
            aGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGACTGGCCGTAACATTCGCTcc                                                                                                                       <  1:310819/147‑1 (MQ=255)
               aTTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTg                                                                                                                   >  8:32361/1‑148 (MQ=255)
                ttGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGt                                                                                                                  <  2:99750/148‑1 (MQ=255)
                   gCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCa                                                                                                               <  7:194606/148‑1 (MQ=255)
                   gCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCa                                                                                                               <  4:344215/148‑1 (MQ=255)
                     gTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAAt                                                                                                             >  1:294566/1‑148 (MQ=255)
                         tACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATgcgc                                                                                                         <  7:32361/148‑1 (MQ=255)
                                aaGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGa                                                                                                  <  3:6221/148‑1 (MQ=255)
                                  gTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGAt                                                                                                 >  8:272682/1‑147 (MQ=255)
                                       aaaaTAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCt                                                                                           >  7:201358/1‑148 (MQ=255)
                                         aaTAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCaa                                                                                                              <  1:330642/127‑1 (MQ=255)
                                         aaTAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCaa                                                                                                              >  2:330642/1‑127 (MQ=255)
                                         aaTAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGa                                                                                         <  3:95598/148‑1 (MQ=255)
                                             aCAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGt                                                                                      >  8:81173/1‑147 (MQ=255)
                                               aGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATgcg                                                                                                          <  6:68863/125‑1 (MQ=255)
                                               aGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATgcg                                                                                                          >  5:68863/1‑125 (MQ=255)
                                                 ccTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATcc                                                                                 >  2:10244/1‑148 (MQ=255)
                                                          tCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGa                                                                          >  7:197857/1‑146 (MQ=255)
                                                                 aTGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGt                                                                 >  1:75536/1‑148 (MQ=255)
                                                                   gATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGAc                                                                         >  4:248871/1‑138 (MQ=255)
                                                                   gATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGAc                                                                         <  3:248871/138‑1 (MQ=255)
                                                                     tAAAAACAGGGGGAAAGGGATGTTGCTTAATATGATTTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTTGCCGTAGCATTCGCTCAAATTTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGAGACAGCCCGGTGTCt                                                             >  5:353918/1‑148 (MQ=255)
                                                                       aaaaCAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGt                                                           <  2:294566/148‑1 (MQ=255)
                                                                        aaaCAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTc                                                          <  7:81173/148‑1 (MQ=255)
                                                                              ggggCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGaa                                                    <  7:272682/148‑1 (MQ=255)
                                                                               gggCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAAc                                                   <  2:75536/148‑1 (MQ=255)
                                                                                  cAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAgg                                                >  8:338284/1‑148 (MQ=255)
                                                                                     gggATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCg                                             <  7:338284/148‑1 (MQ=255)
                                                                                     gggATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCg                                             <  8:197857/148‑1 (MQ=255)
                                                                                      ggATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGc                                              >  7:357560/1‑146 (MQ=255)
                                                                                         tGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGc                                         <  3:235128/148‑1 (MQ=255)
                                                                                          gTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTc                                                                                                                            >  2:37814/1‑64 (MQ=255)
                                                                                          gTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTc                                                                                                                            <  1:37814/64‑1 (MQ=255)
                                                                                          gTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAAc                                                   <  6:108045/137‑1 (MQ=255)
                                                                                          gTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCACTGGCGACAGCCCGGTGGCTGTCAGTGAAc                                                   >  5:108045/1‑137 (MQ=255)
                                                                                           ttGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCAc                                       <  1:10244/148‑1 (MQ=255)
                                                                                                        tgtgtgGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCa                          <  8:357560/148‑1 (MQ=255)
                                                                                                           gtgGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAAc                                                   >  1:97544/1‑120 (MQ=255)
                                                                                                           gtgGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAAc                                                   <  2:97544/120‑1 (MQ=255)
                                                                                                             ggTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACa                                                  <  2:188235/119‑1 (MQ=255)
                                                                                                             ggTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACa                                                  >  1:188235/1‑119 (MQ=255)
                                                                                                                cgtcAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGc                                 >  7:396778/1‑133 (MQ=255)
                                                                                                                cgtcAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGc                                 <  8:396778/133‑1 (MQ=255)
                                                                                                                    aGCTGTAGGCAAACAGTTTTGTCCTGTCCGTAACCGTCGCTCCACTGTTCCATGCGAAGGCCGATGAGCCAGAAGTAATCCCTGGCGACAGCCCTGTGGGTTTCAGTGAACAGTGCGAGTAACTGCCGCAGGCGCAAGCCACGGCaat              <  6:353918/148‑1 (MQ=255)
                                                                                                                             cAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCggg     >  6:14349/1‑148 (MQ=255)
                                                                                                                                tCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGAt  >  3:36552/1‑148 (MQ=255)
                                                                                                                                 cAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGt                                                                 >  4:205477/1‑84 (MQ=255)
                                                                                                                                 cAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGt                                                                 <  3:205477/84‑1 (MQ=255)
                                                                                                                                     ttGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCaata             >  6:361272/1‑132 (MQ=255)
                                                                                                                                     ttGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCaata             <  5:361272/132‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CATACTTTCGTGAGTATTGGCGTTGTACAGGCAAGTCGTAAAATAACAGCCTGGCTATTCAGAGTATGATAAAAACAGGGGGCAAGGGATGTTGCTTAATATGATGTGTGGTCGTCAGCTGTCGGCAATCAGTTTGTGCCTGGCCGTAACATTCGCTCCACTGTTCAATGCGCAGGCCGATGAGCCTGAAGTAATCCCTGGCGACAGCCCGGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGAT  >  NC_000913/980959‑981234

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