Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,721,945 T→A 7.4% L179Q (CTG→CAG)  PP_1517 → RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,721,9450TA 7.4% 175.1 / 4.2 54L179Q (CTG→CAG) PP_1517RND family transporter
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base A (2/2);  total (30/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.29e-01

TTCAACGACGAGTTCGGCACCACCTTCGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGTTTCTTCGAGACGCCCAGCGAGCGCT  >  NC_002947/1721806‑1722094
                                                                                                                                           |                                                                                                                                                     
ttcAACGACGAGTTCGGCACCACCTTCGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTg                                                                                                                                             >  2:1202652/1‑150 (MQ=255)
   aaCGACGAGTTCGGCACCACCTTCGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAg                                                                                                                                          >  2:1208729/1‑150 (MQ=255)
    acgaccAGTTCGGCAACACCTCAGGCCCTATCTATCGAAAAACTGGCGAATGCTTGGACTATGCTGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAgg                                                                                                                                         <  2:1084564/150‑1 (MQ=255)
                   ccaccTTCGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGa                                                                                                                          <  1:362753/150‑1 (MQ=255)
                       cTTCGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGa                                                                                                                      <  1:1208729/150‑1 (MQ=255)
                          cGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTATGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCt                                                                                                                   <  1:355433/150‑1 (MQ=255)
                           ggCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTg                                                                                                                  >  1:1229173/1‑150 (MQ=255)
                            gCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGt                                                                                                                 >  1:1158167/1‑150 (MQ=255)
                             cAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTc                                                                                                                >  1:463450/1‑150 (MQ=255)
                                         gcTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCt                                                                                                    >  2:597925/1‑150 (MQ=255)
                                                 gCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCAccc                                                                                            >  2:285695/1‑150 (MQ=255)
                                                       gCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCg                                                                                      >  1:504600/1‑150 (MQ=255)
                                                         tttGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTg                                                                                    >  1:1210406/1‑150 (MQ=255)
                                                            gACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCg                                                                                 >  2:733756/1‑150 (MQ=255)
                                                            gACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGCCGATCAGATCGACCttgc                                                                                                                                           >  1:344030/1‑89 (MQ=255)
                                                            gACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGCCGATCAGATCGACCttgc                                                                                                                                           <  2:344030/92‑4 (MQ=255)
                                                             aCTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGc                                                                                                                                            >  2:1202199/1‑90 (MQ=255)
                                                             aCTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGc                                                                                                                                            <  1:1202199/90‑1 (MQ=255)
                                                              cTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCt                                                                               <  2:1158167/150‑1 (MQ=255)
                                                               tATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTg                                                                              >  1:706676/1‑150 (MQ=255)
                                                                 tGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGa                                                                            >  2:99026/1‑150 (MQ=255)
                                                                    ggTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCa                                                                                                                                           >  1:519032/1‑84 (MQ=255)
                                                                    ggTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCa                                                                                                                                           <  2:519032/84‑1 (MQ=255)
                                                                       gCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGt                                                                      <  2:463450/150‑1 (MQ=255)
                                                                            aGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGa                                                                 >  2:461666/1‑150 (MQ=255)
                                                                              gACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAg                                                               <  2:504600/150‑1 (MQ=255)
                                                                                 tACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGATGCCGTGCAGAAGGCg                                                            >  2:989479/1‑150 (MQ=255)
                                                                                    gCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTg                                                         >  1:998766/1‑150 (MQ=255)
                                                                                     ccGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGc                                                        >  1:439431/1‑150 (MQ=255)
                                                                                         ccGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAaga                                                    <  2:1037164/150‑1 (MQ=255)
                                                                                            tATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGCCGATCAGATCGACCtt                                                                                                                                             >  1:1205655/1‑57 (MQ=37)
                                                                                            tATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGCCGATCAGATCGACCtt                                                                                                                                             <  2:1205655/58‑2 (MQ=37)
                                                                                              tccagatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGa                                               >  2:1292769/1‑150 (MQ=255)
                                                                                                cagatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCAGTGCAGaa                                                                >  1:1211262/1‑131 (MQ=255)
                                                                                                cagatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCAGTGCAGaa                                                                <  2:1211262/131‑1 (MQ=255)
                                                                                                 agatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAAcg                                            <  2:998766/150‑1 (MQ=255)
                                                                                                 agatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAAcg                                            >  2:816147/1‑150 (MQ=255)
                                                                                                  gatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGc                                                                                       >  1:1048214/1‑106 (MQ=255)
                                                                                                  gatccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGc                                                                                       <  2:1048214/106‑1 (MQ=255)
                                                                                                    tccagCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCgg                                         <  2:1424852/150‑1 (MQ=255)
                                                                                                       agCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGc                                                             <  2:93378/127‑1 (MQ=255)
                                                                                                       agCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGc                                                             >  1:93378/1‑127 (MQ=255)
                                                                                                       agCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGa                                      <  1:285695/150‑1 (MQ=255)
                                                                                                       agCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGa                                      <  1:989479/150‑1 (MQ=255)
                                                                                                           gCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCAc                                  >  1:1025124/1‑150 (MQ=255)
                                                                                                           gCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCAc                                  <  1:99026/150‑1 (MQ=255)
                                                                                                                  gTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGt                           >  1:320747/1‑150 (MQ=255)
                                                                                                                   tCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGtt                          >  1:619826/1‑150 (MQ=255)
                                                                                                                   tCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGtt                          <  1:1202652/150‑1 (MQ=255)
                                                                                                                   tCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGtt                          >  2:1193617/1‑150 (MQ=255)
                                                                                                                        gATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGTTTCTTC                     <  2:619826‑M1/150‑1 (MQ=255)
                                                                                                                         aTGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGTTTCTTCG                    <  2:706676‑M1/150‑1 (MQ=255)
                                                                                                                          tgtgGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGTTTCTTCGA                   >  2:1167322‑M1/1‑150 (MQ=255)
                                                                                                                                    gTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGAACAACCTCAAGCTGGTCATCCTTGGCGTTCCGCTTGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGTTTCTTCGAGACGCCCagc         <  2:1229173/150‑1 (MQ=255)
                                                                                                                                           tGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCCAGAGAAGAAACAGACGAACACTGCCGCTTGCGCCAAGACGCCCAGCGCGCGCt  >  2:373925/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
TTCAACGACGAGTTCGGCACCACCTTCGGCAATATCTATGCGCTGACTGGCGAAGGCTTTGACTATGCGGTGCTCAAGGACTACGCCGACCGTATCCAGATCCAGCTGCAGCGGGTCAAGGATGTGGGCAAGGTCGATCTGATCGGCCTGCAGGACGAGAAAGTCTGGATCGAGCTGTCCAACCTCAAGCTGGCCACCCTCGGCGTGCCGCTGGAAGCCGTGCAGAAGGCGTTGCAAGAGCAGAACGCGGTGAGCACCGCCGGTTTCTTCGAGACGCCCAGCGAGCGCT  >  NC_002947/1721806‑1722094

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