Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,878,603 C→T 8.0% A242V (GCC→GTC)  PP_1684 → transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,878,6030CT 8.0% 138.9 / 2.8 50A242V (GCC→GTC) PP_1684transporter
Reads supporting (aligned to +/- strand):  ref base C (26/19);  new base T (2/2);  total (28/22)
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.86e-01

GATCGCCCGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCGGGTTGTACCTGCT  >  NC_002947/1878458‑1878750
                                                                                                                                                  |                                                                                                                                                   
gATCGCCCGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTa                                                                                                                                                  >  7:339531/1‑149 (MQ=255)
  tCGCCCGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTAcc                                                                                                                                                >  3:257521/1‑149 (MQ=255)
     cccGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGTCAACCAgg                                                                                                                                             >  8:303552/1‑149 (MQ=255)
     cccGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTgg                                                                                                                                             <  7:302232/149‑1 (MQ=255)
     cccGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTgg                                                                                                                                             <  8:124772/149‑1 (MQ=255)
       cGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGtt                                                                                                                                           <  6:82152/149‑1 (MQ=255)
           aCCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGAc                                                                                                                                        <  2:147402/148‑1 (MQ=255)
           aCCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGAc                                                                                                                                        <  3:222592/148‑1 (MQ=255)
              cAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGTCAACCAGGTAGGCCAGc                                                                                                                                    <  8:200474/149‑1 (MQ=255)
                gTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGTCAACCAGGTAGGCCAGCag                                                                                                                                  <  7:303552/149‑2 (MQ=255)
                                   ccTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCCCCGGTa                                                                                                               >  6:308879/1‑149 (MQ=255)
                                         ccAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGc                                                                                                                          <  8:306609/132‑1 (MQ=255)
                                         ccAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGc                                                                                                                          >  7:306609/1‑132 (MQ=255)
                                                     cgAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTgg                                                                                             >  6:146604/1‑149 (MQ=255)
                                                         gCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCgg                                                                                         <  4:293958/149‑1 (MQ=255)
                                                         gCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCg                                                                                          <  4:257521/148‑1 (MQ=255)
                                                           gTGCACGCACTCCGGCGTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACTACGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGc                                                                                        <  1:307004/149‑1 (MQ=255)
                                                                    cTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTAGGCCTGTTCAACa                                                                              >  5:153725/1‑149 (MQ=255)
                                                                     tCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACAt                                                                             >  3:156321/1‑149 (MQ=255)
                                                                       cGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACAtcg                                                                           >  4:261360/1‑149 (MQ=255)
                                                                                       gCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCgc                                                           >  6:192867/1‑149 (MQ=255)
                                                                                        cACTGGGCTTTTTTGTCTTTTGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGACTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTTGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCgct                                                          <  5:365219/149‑1 (MQ=255)
                                                                                         aCTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCgctg                                                         <  3:261360/149‑1 (MQ=255)
                                                                                                ttttttGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGtt                                                  <  7:347418/149‑1 (MQ=255)
                                                                                                  ttttGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTgg                                                >  6:233734/1‑149 (MQ=255)
                                                                                                  ttttGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTgg                                                >  1:74325/1‑149 (MQ=255)
                                                                                                   tttGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGtt                                                                                                                                           <  7:339978/57‑1 (MQ=255)
                                                                                                   tttGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGtt                                                                                                                                           >  8:339978/1‑57 (MQ=255)
                                                                                                     tGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTgggc                                              <  5:308879/148‑1 (MQ=255)
                                                                                                          gtgGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGggcgggcg                                         >  4:67108/1‑148 (MQ=255)
                                                                                                           tgGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGACGGTACCTTCACCGCTGGCTGGTTGCGCGGGCGCa                                       >  8:193225/1‑149 (MQ=255)
                                                                                                              cTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGt                                    >  7:21705/1‑149 (MQ=255)
                                                                                                                       cGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCcgc                           >  8:226977/1‑149 (MQ=255)
                                                                                                                        gTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTc                                                                                  <  3:218392/93‑1 (MQ=255)
                                                                                                                        gTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTc                                                                                  >  4:218392/1‑93 (MQ=255)
                                                                                                                        gTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCcgcg                          >  6:12832/1‑149 (MQ=255)
                                                                                                                         tATTGATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCgcgc                         <  4:12503/149‑1 (MQ=255)
                                                                                                                         tATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCgcgc                         <  5:233734/149‑1 (MQ=255)
                                                                                                                           ttCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCct                       >  7:277477/1‑149 (MQ=255)
                                                                                                                              aTTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACAtcg                                                                           >  7:82628/1‑94 (MQ=255)
                                                                                                                              aTTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACAtcg                                                                           <  8:82628/94‑1 (MQ=255)
                                                                                                                               ttGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAAc                                                                               >  4:372682/1‑89 (MQ=255)
                                                                                                                                tGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGa                  >  6:118395/1‑149 (MQ=255)
                                                                                                                                 ggCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGa                  <  5:192867/148‑1 (MQ=255)
                                                                                                                                   cGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCg               <  3:67108/149‑1 (MQ=255)
                                                                                                                                     ttCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCGGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCggg             >  5:226919/1‑149 (MQ=255)
                                                                                                                                        aCCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCGGGTTg          >  6:77771/1‑149 (MQ=255)
                                                                                                                                           ctgctgGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCGGGTTGTa        >  6:320759/1‑149 (MQ=255)
                                                                                                                                           atgctgGTCAACCAGGTAGGCCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCGGGTTGTa        >  3:317676/2‑149 (MQ=255)
                                                                                                                                                 gCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCGGGTTGTACctgct  >  5:78765/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATCGCCCGTTACCCAGTCACGGTGCCGAGCAGACCCTGGGCCAGGCGCTGCGCGAGGCGTGCACGCACTCCGGCTTCTGGTTGCTGGCACTGGGCTTTTTTGTCTGTGGCTTCCAGGTCGTATTCATTGGCGTTCACCTGCTGGCCTACCTGGTTGACCAGCATCTGGCTGCGGCCACCGGTACGACAGTGCTGGCGCTGGTCGGCCTGTTCAACATCGTCGGTACCTTCACCGCTGGCTGGTTGGGCGGGCGCATGTCCAAGCCGCGCCTGCTGACCGGGTTGTACCTGCT  >  NC_002947/1878458‑1878750

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: