Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 708,502 G→T 10.0% E302* (GAG→TAG) ‡ murJ → lipid II flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947708,5020GT10.0% 151.1 / 7.2 50E302* (GAG→TAG) ‡murJlipid II flippase
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base T (2/3);  total (25/25)
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.18e-01

CCTGGTGGCCGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATGCGGCGATGA  >  NC_002947/708366‑708646
                                                                                                                                        |                                                                                                                                                
cctGGTGGCCGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGTACTCCTcgcg                                                                                                                                        <  2:272251/147‑1 (MQ=255)
cctGGGGGCCGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGTACTCCTcgcgg                                                                                                                                       <  5:128386/148‑2 (MQ=255)
  tggtggCCGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAAt                                                                                                                                     >  6:55625/1‑148 (MQ=255)
    gtggCCGGCTCGGTGTCTTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGTTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATtct                                                                                                                                  <  2:201641/149‑1 (MQ=255)
     tggCCGGCTCGGTGTCGTGGATGTACTACGCCGATCGCCTCATGGAGCTGCCTTCCGGTGAGCTGGGCGTTGCCCTGGGCACTATCCTGATGCCGACACTGGCAAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATtcac                                                                                                                                 >  8:164648/1‑147 (MQ=255)
     tggCCGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATtctc                                                                                                                                 <  4:84937/149‑1 (MQ=255)
         cGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAAtt                                                                                                                                    >  7:118984/1‑142 (MQ=255)
         cGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAAtt                                                                                                                                    <  8:118984/142‑1 (MQ=255)
            cTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTggg                                                                                                                          <  5:55625/149‑1 (MQ=255)
               ggTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGcc                                                                                                                       >  3:346674/1‑149 (MQ=255)
                 tgtCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTg                                                                                                                     >  2:73066/1‑149 (MQ=255)
                    cGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTgcgc                                                                                                                  >  2:52970/1‑149 (MQ=255)
                        gATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCtgt                                                                                                              <  5:349868/149‑1 (MQ=255)
                          tGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGc                                                                                                            <  3:240478/149‑1 (MQ=255)
                           gTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTggg                                                                                                                          <  7:157385/134‑1 (MQ=255)
                           gTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTggg                                                                                                                          >  8:157385/1‑134 (MQ=255)
                              ctacGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTTCCGACACTTGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTc                                                                                                         <  2:334968/148‑1 (MQ=255)
                                          ccTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGc                                                                                            >  5:182226/1‑149 (MQ=255)
                                               tGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGc                                                                                       <  1:52970/149‑1 (MQ=255)
                                                           ccGGTGTGCTGGGCGTTGCCCTGGGCAATATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAATGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGtt                                                                            <  7:164648/148‑1 (MQ=255)
                                                            cGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTg                                                                           <  7:44548/148‑1 (MQ=255)
                                                            cGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTCGGGCCTGCGCCTGTCCTTCCTGCTGGTAATGCCCTGCACCCTGGCGTTgg                                                                          >  3:190649/1‑149 (MQ=255)
                                                             ggTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGc                                                                         >  6:120612/1‑149 (MQ=255)
                                                              gtgtGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTggg                                                                                                                          >  5:304772/1‑99 (MQ=255)
                                                              gtgtGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTggg                                                                                                                          <  6:304772/99‑1 (MQ=255)
                                                                     gggCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTATAGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTAGGCCATTCTTg                                                                 >  2:329117/1‑149 (MQ=255)
                                                                          ttGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGaa                                                            >  3:255296/1‑149 (MQ=255)
                                                                               cTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCgagga                                                                                                                                                 <  4:335712/59‑1 (MQ=255)
                                                                               cTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCgagga                                                                                                                                                 >  3:335712/1‑59 (MQ=255)
                                                                               cTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGAGTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCt                                                       >  1:119915/1‑149 (MQ=255)
                                                                               cTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTg                                                                                                                     <  4:26031/87‑1 (MQ=255)
                                                                                  ggCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTg                                                      <  4:182512/147‑1 (MQ=255)
                                                                                       tATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGtag                                               >  8:317778/1‑147 (MQ=255)
                                                                                        aTCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTc                                                                    <  8:113022/127‑1 (MQ=255)
                                                                                        aTCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTc                                                                    >  7:113022/1‑127 (MQ=255)
                                                                                              ctgcCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTgg                                                                          <  4:246445/115‑1 (MQ=255)
                                                                                              ctgcCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTgg                                                                          >  3:246445/1‑115 (MQ=255)
                                                                                              ctgcCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGtt                                        <  2:119915/149‑1 (MQ=255)
                                                                                                 ccGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGttt                                       <  4:346674/147‑1 (MQ=255)
                                                                                                  cgtctcTGTCCAAGACCTACTCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCaa                                    <  4:190649/144‑1 (MQ=255)
                                                                                                  cGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTCGCTTTGTTTCaa                                    <  7:317778/149‑1 (MQ=255)
                                                                                                      acTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCa                                                                       <  6:238699/110‑1 (MQ=255)
                                                                                                      acTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCa                                                                       >  5:238699/1‑110 (MQ=255)
                                                                                                           ccAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAg                           >  2:69577/1‑149 (MQ=255)
                                                                                                                   tACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAgcgc                   >  4:243623/1‑149 (MQ=255)
                                                                                                                              ggACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATgcgg        >  2:332921/1‑149 (MQ=255)
                                                                                                                                ttcgcgAGTACTCCTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATgcggc       >  1:314073/3‑148 (MQ=255)
                                                                                                                                 ccGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATGCGGCGa     >  5:344091/1‑149 (MQ=255)
                                                                                                                                   gcgAGTACTCCTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCt                                             <  4:326891/107‑1 (MQ=255)
                                                                                                                                   gcgAGTACTCCTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCt                                             >  3:326891/1‑107 (MQ=255)
                                                                                                                                    cgAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATGCGGCGATGa  <  5:120612/149‑1 (MQ=255)
                                                                                                                                    cgAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATGCGGCGATGa  >  2:99301/1‑149 (MQ=255)
                                                                                                                                        gagTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGt                                         <  7:257965/106‑1 (MQ=255)
                                                                                                                                        gagTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGt                                         >  8:257965/1‑106 (MQ=255)
                                                                                                                                        |                                                                                                                                                
CCTGGTGGCCGGCTCGGTGTCGTGGATGTACTACGCCGACCGCCTCATGGAGCTGCCTTCCGGTGTGCTGGGCGTTGCCCTGGGCACTATCCTGCTGCCGACACTGGCCAAGACCTACGCCAACAAGGACCGCGAGGAGTACTCGCGAATTCTCGACTGGGGCCTGCGCCTGTGCTTCCTGCTGGTACTGCCCTGCACCCTGGCGTTGGCCATTCTTGCCGAACCGCTGACCGTTGCTTTGTTTCAATATGGCAAGTTCAGCGCGTTTGATGCGGCGATGA  >  NC_002947/708366‑708646

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