Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,724,453 T→A 10.2% intergenic (+20/+29) PP_1517 → / ← PP_1518 RND family transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,724,4530TA10.2% 134.0 / 6.3 49intergenic (+20/+29)PP_1517/PP_1518RND family transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/19);  new base A (3/2);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

GCGTAGTGTGTTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCTTCTCGGTC  >  NC_002947/1724310‑1724599
                                                                                                                                               |                                                                                                                                                  
gcgTAGTGTGTTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGccc                                                                                                                                               >  6:53575/1‑149 (MQ=255)
   tAGTGTGTTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGccctcc                                                                                                                                            >  3:186122/1‑149 (MQ=255)
    aGTGTGTTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCAGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGccctccc                                                                                                                                           >  3:52166/1‑149 (MQ=255)
     gtgtgtTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGccctccct                                                                                                                                          >  7:113889/1‑149 (MQ=255)
      tgtgtTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCtt                                                                                                                                         >  4:79464/1‑149 (MQ=255)
          ttCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTTGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTcgcg                                                                                                                                     <  5:257121/149‑1 (MQ=255)
            cTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCgg                                                                                                                                    <  2:19003/148‑1 (MQ=255)
             tACGGGCCGATGGCGGTGGCGATCATGGACGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCggg                                                                                                                                   >  1:21725/1‑148 (MQ=255)
                gggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGGTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTaaa                                                                                                                               >  4:4237/1‑149 (MQ=255)
                       tGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTc                                                                                                                        <  4:186122/149‑1 (MQ=255)
                           ggTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                    <  6:238274/149‑1 (MQ=255)
                            gTGGCGAGCTTGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGCCTGAAGGGCGCAGGGTGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTGCTACa                                                                                                                   <  3:135109/149‑1 (MQ=255)
                                 gATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGccctcc                                                                                                                                            >  3:81847/1‑119 (MQ=255)
                                 gATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAgcgcg                                                                                                              >  3:275804/1‑149 (MQ=255)
                                   tCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGcc                                                                                                            <  5:278533/149‑1 (MQ=255)
                                     aTGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGa                                                                                                          <  5:53575/149‑1 (MQ=255)
                                      tGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGaa                                                                                                         <  8:113889/149‑1 (MQ=255)
                                                 aTCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTgg                                                                                              <  4:123793/149‑1 (MQ=255)
                                                     tGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAg                                                                                          <  4:275804/149‑1 (MQ=255)
                                                       gCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAgg                                                                                         <  3:79464/148‑1 (MQ=255)
                                                           cGGTGCGTCCACTGCTGTGCCTGCCGGCGTTGTAAGCGGCGTGGGTCAGTGTGACGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGGAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGt                                                                                    <  8:324545/149‑1 (MQ=255)
                                                               gCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAg                                                                                  >  8:337866/1‑147 (MQ=255)
                                                               gCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAg                                                                                  <  7:337866/147‑1 (MQ=255)
                                                                cTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAgccg                                                                               >  8:122405/1‑149 (MQ=255)
                                                                  tACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAgccgcc                                                                             >  3:136526/1‑149 (MQ=255)
                                                                            ttCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGt                                                                   >  1:211699/1‑149 (MQ=255)
                                                                                 gCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTgg                                                              <  3:4237/149‑1 (MQ=255)
                                                                                      cGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAAt                                                         <  4:136526/149‑1 (MQ=255)
                                                                                            aTGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACGTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTg                                                   >  2:167086/1‑149 (MQ=255)
                                                                                               cggcgTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTc                                                <  7:122405/149‑1 (MQ=255)
                                                                                                 gcgTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCt                                              >  3:112572/1‑149 (MQ=255)
                                                                                                  cgTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTg                                             >  2:256263/1‑149 (MQ=255)
                                                                                                      gTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTg                                         >  6:112867/1‑149 (MQ=255)
                                                                                                            ggTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGat                                    >  2:124896/1‑148 (MQ=255)
                                                                                                                 agaagGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCAtctc                              >  4:124696/1‑149 (MQ=255)
                                                                                                                  gaagGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCAtctc                              >  2:511/1‑148 (MQ=255)
                                                                                                                         tGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGc                      >  8:348054/1‑149 (MQ=255)
                                                                                                                               ccgcgaagggAGGGCTAAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATc                  >  4:158165/8‑147 (MQ=255)
                                                                                                                                cGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAg               >  1:297691/1‑149 (MQ=255)
                                                                                                                                   gaagggAGGGCTAAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAgcgc            >  6:61560/4‑149 (MQ=255)
                                                                                                                                    gggggAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGcc           >  5:25496/1‑149 (MQ=255)
                                                                                                                                     ggggAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCt          >  8:20453/1‑149 (MQ=255)
                                                                                                                                      gggAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCtt         <  4:112572/149‑1 (MQ=255)
                                                                                                                                      gggAGGGCTAAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCg                                                                      <  7:29782/88‑1 (MQ=255)
                                                                                                                                      gggAGGGCTAAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCg                                                                      >  8:29782/1‑88 (MQ=255)
                                                                                                                                        gAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATCTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCTtct       >  5:52406/1‑149 (MQ=255)
                                                                                                                                           ggCTAAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCTTCTCgg    <  3:158165/149‑1 (MQ=255)
                                                                                                                                            gCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCTTCTCGGt   <  1:125177/149‑1 (MQ=255)
                                                                                                                                             cTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCTTCTCGGTc  <  1:124896/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCGTAGTGTGTTCTACGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTGGCCACGGTGCTTACACTGCTGTTCCTGCCGGCGTTGTATGCGGCGTGGTTCAGGGTGAAGAAGGGCTGAAGGCCGCAGGGGGAGGGCTTAGCCCTCCCTTCGCGGGTAAACCCGCTCCTACAGCGCGCCGAAGACCTTCTTGGCCAGGCTGGTAGCCGCCGCTGCCGGGTTCTGGCGAATGTTCTGCTCCTGTTTGCCGATCATCTCGAACAGGCCATCCAGCGCCTTCTCGGTC  >  NC_002947/1724310‑1724599

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