Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,101,078 G→A 7.3% V127I (GTC→ATC) ‡ PP_2710 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,101,0780GA 7.3% 160.1 / 2.8 55V127I (GTC→ATC) ‡PP_2710hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (2/2);  total (27/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 = 7.52e-01

GAGGAGGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAG  >  NC_002947/3100931‑3101221
                                                                                                                                                   |                                                                                                                                               
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                                  >  3:168841/5‑123 (MQ=255)
                                                                                                                                        cgaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCt                                  <  4:168841/119‑1 (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                                                                    >  2:151837/1‑84 (MQ=255)
                                                                                                                                             aaCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAg                                                                    <  1:151837/84‑1 (MQ=255)
                                                                                                                                               cccGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAg  >  6:136764/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GAGGAGGAATACGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAG  >  NC_002947/3100931‑3101221

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