Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,101,079 T→C 7.4% V127A (GTC→GCC) ‡ PP_2710 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,101,0790TC 7.4% 154.8 / 2.7 54V127A (GTC→GCC) ‡PP_2710hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base C (2/2);  total (26/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.54e-01

GAGGAGGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCACATG  >  NC_002947/3100931‑3101227
                                                                                                                                                    |                                                                                                                                                    
gaggagGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGt                                                                                                                                                      >  5:32418/1‑149 (MQ=255)
 aggaggAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTc                                                                                                                                                     >  6:328474/1‑149 (MQ=255)
 aggaggAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCT‑ATCGCCCGTGGCAACCCGGTc                                                                                                                                                     <  7:173649/148‑1 (MQ=255)
   gaggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCgg                                                                                                                                                   <  6:2095/149‑1 (MQ=255)
    aggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCggg                                                                                                                                                  <  5:328474/149‑1 (MQ=255)
    aggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCggg                                                                                                                                                  <  7:55481/149‑1 (MQ=255)
     ggaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgt                                                                                                                                                 >  5:5055/1‑149 (MQ=255)
      gaATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgtt                                                                                                                                                >  6:234941/1‑149 (MQ=255)
           cGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTc                                                                                                                                            <  1:10749/148‑1 (MQ=255)
           cGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCg                                                                                                                                           <  3:196752/149‑1 (MQ=255)
                  gAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGt                                                                                                                                    >  6:312696/1‑149 (MQ=255)
                           aTCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCt                                                                                                                           <  6:32418/149‑1 (MQ=255)
                            tCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTc                                                                                                                          <  7:75613/149‑1 (MQ=255)
                               tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCg                                                                                                                                           >  6:228850/1‑129 (MQ=255)
                               tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCg                                                                                                                                           >  6:228851/1‑129 (MQ=255)
                               tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCg                                                                                                                                           <  5:228850/129‑1 (MQ=255)
                               tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCg                                                                                                                                           <  5:228851/129‑1 (MQ=255)
                               tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCg                                                                                                                         >  7:264150/1‑147 (MQ=255)
                                 aaggACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAg                                                                                                                      <  6:5055/145‑1 (MQ=255)
                                 aaCGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCATCGTGCTCGaa                                                                                                                       <  6:63190/147‑1 (MQ=255)
                                     aCATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCAc                                                                                                                  >  3:51613/1‑148 (MQ=255)
                                                   ggcggTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCggg                                                                                                                                                  >  1:132618/1‑102 (MQ=255)
                                                   ggcggTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCggg                                                                                                                                                  <  2:132618/102‑1 (MQ=255)
                                                               gAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAggcggc                                                                                       >  6:312585/1‑149 (MQ=255)
                                                                   ccGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCggcgggc                                                                                   <  5:234941/149‑1 (MQ=255)
                                                                         gcgATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATc                                                                             <  5:312696/149‑1 (MQ=255)
                                                                                       gTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGc                                                               <  8:264150/149‑1 (MQ=255)
                                                                                            aTCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCAc                                                                                                                  >  7:155862/1‑93 (MQ=255)
                                                                                            aTCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCAc                                                                                                                  <  8:155862/93‑1 (MQ=255)
                                                                                            aTCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgc                                                          <  6:258647/149‑1 (MQ=255)
                                                                                               cATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCAc                                                                                                                  >  4:120208/1‑90 (MQ=255)
                                                                                               cATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCAc                                                                                                                  <  3:120208/90‑1 (MQ=255)
                                                                                                         aTGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGtt                                             >  5:61445/1‑149 (MQ=255)
                                                                                                          tGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGtt                                             >  2:211270/1‑148 (MQ=255)
                                                                                                           ggTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGtt                                             >  7:39754/1‑147 (MQ=255)
                                                                                                            gTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAg                                          >  6:260586/1‑149 (MQ=255)
                                                                                                             tCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGc                                         >  8:51025/1‑149 (MQ=255)
                                                                                                             tCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGc                                         >  5:144254/1‑149 (MQ=255)
                                                                                                               gCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGt                                                                                 <  6:243998/107‑1 (MQ=255)
                                                                                                                    ccTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGAcc                                                    >  2:111634/1‑131 (MQ=255)
                                                                                                                     cTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAg                                 <  5:312585/149‑1 (MQ=255)
                                                                                                                       gTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCt                               >  7:296569/1‑149 (MQ=255)
                                                                                                                                 aTCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTg                      <  8:39754/148‑1 (MQ=255)
                                                                                                                                     cccGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTc                 <  7:51025/149‑1 (MQ=255)
                                                                                                                                      ccGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCACTGGCCACCCAGGCGGCGGGCGTTATCCAGAAAACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTc                 >  2:247061/1‑148 (MQ=255)
                                                                                                                                       cGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGAGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGa               >  3:213437/1‑149 (MQ=255)
                                                                                                                                        cgaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCt                                        <  4:168841/119‑1 (MQ=255)
                                                                                                                                        cgaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCt                                        >  3:168841/5‑123 (MQ=255)
                                                                                                                                        cgaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGt                       <  2:233411/136‑1 (MQ=255)
                                                                                                                                        cgaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGt                       >  1:233411/5‑140 (MQ=255)
                                                                                                                                           gCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAg           <  6:61445/149‑1 (MQ=255)
                                                                                                                                            cAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAgg          <  1:211270/149‑1 (MQ=255)
                                                                                                                                             aaCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAg                                                                          <  1:151837/84‑1 (MQ=255)
                                                                                                                                             aaCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAg                                                                          >  2:151837/1‑84 (MQ=255)
                                                                                                                                               cccGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAg        >  6:136764/1‑148 (MQ=255)
                                                                                                                                                    tCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCACATg  <  8:296569/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGGAGGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCACATG  >  NC_002947/3100931‑3101227

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: ‑