Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,486,216 A→C 6.8% S30R (AGC→CGC)  PP_3978 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,486,2160AC 6.8% 181.8 / 2.6 58S30R (AGC→CGC) PP_3978hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base C (2/2);  total (29/29)
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.99e-01

GGCAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGC  >  NC_002947/4486069‑4486364
                                                                                                                                                   |                                                                                                                                                    
ggCAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGTAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAg                                                                                                                                                     >  2:282477/1‑149 (MQ=255)
  cAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCt                                                                                                                                                   >  3:160468/1‑149 (MQ=255)
    aGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGc                                                                                                                                                 <  6:152958/149‑1 (MQ=255)
    aGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGc                                                                                                                                                 <  6:157302/149‑1 (MQ=255)
      aTCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCa                                                                                                                                               <  1:115817/149‑1 (MQ=255)
         aaTGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                              >  6:306881/1‑145 (MQ=255)
         aaTGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                              <  5:306881/147‑3 (MQ=255)
              tgcGCCACTACAATTACATCGTTGCCGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                              <  3:193202/142‑3 (MQ=255)
              tgcGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                              >  4:193202/1‑140 (MQ=255)
                     cTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGTAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGAt                                                                                                                                <  1:282477/149‑1 (MQ=255)
                       aCAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATgg                                                                                                                              <  4:160468/149‑1 (MQ=255)
                           ttACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGc                                                                                                                                                 >  1:146772/1‑126 (MQ=255)
                           ttACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGc                                                                                                                                                 <  2:146772/126‑1 (MQ=255)
                                        gacCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGc                                                                                                             >  5:260864/1‑149 (MQ=255)
                                        gacCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATCAGCTCCACAGGCACACAAAGCATGACCCGCCAGCGAGCTGCCATGAACCAGCCCTGATGGAGCAGTTTCACCGTGGc                                                                                                             >  4:129984/1‑149 (MQ=255)
                                           ccGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCttt                                                                                                          <  6:205350/149‑1 (MQ=255)
                                           ccGATCATTCCAAACGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCttt                                                                                                          <  3:129984/149‑1 (MQ=255)
                                            cGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTa                                                                                                         <  1:128334/149‑1 (MQ=255)
                                             gATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGAAGCATTTGAACCGTGGCttt                                                                                                          >  6:242006/1‑147 (MQ=255)
                                                     cATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGAt                                                                                                <  7:231375/149‑1 (MQ=255)
                                                     cATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGAt                                                                                                >  5:91159/1‑149 (MQ=255)
                                                      aTAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATc                                                                                               >  6:103804/1‑149 (MQ=255)
                                                           tATGGGCCACCCGTGGTTTAGCCACAAACGCAACATCGATGCCAGCACCTATGAGCTCCACAGGCACACAAAGCAAGACCCGTCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTGTAAGACGGGATCGATAc                                                                                          <  5:242006/149‑1 (MQ=255)
                                                               tgCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGa                                                                                       >  8:85680/1‑148 (MQ=255)
                                                                gCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGa                                                                                       >  8:65207/1‑147 (MQ=255)
                                                                 ccACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAg                                                                                    >  7:109398/1‑149 (MQ=255)
                                                                     ccGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTggg                                                                                >  2:200977/1‑149 (MQ=255)
                                                                          gtgtAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACa                                                                           <  6:260864/149‑1 (MQ=255)
                                                                             tAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATgg                                                                        >  8:127100/1‑149 (MQ=255)
                                                                                 cacaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCATGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGt                                                                     >  6:175336/1‑148 (MQ=255)
                                                                                  acaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTaa                                                                   <  7:127100/149‑1 (MQ=255)
                                                                                     aaCGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCg                                                                >  8:234622/1‑149 (MQ=255)
                                                                                      aCGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAg                                                                                    <  4:216072/128‑1 (MQ=255)
                                                                                      aCGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAg                                                                                    >  3:216072/1‑128 (MQ=255)
                                                                                           aCATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGc                                                          >  4:101092/1‑149 (MQ=255)
                                                                                           aCATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGc                                                          <  1:200977/149‑1 (MQ=255)
                                                                                              tCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                                       <  8:109398/149‑1 (MQ=255)
                                                                                              tCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                                       <  5:103804/149‑1 (MQ=255)
                                                                                                      agcagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCtt                                               >  6:250988/1‑149 (MQ=255)
                                                                                                       gcagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                                       <  6:32470/140‑1 (MQ=255)
                                                                                                       gcagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAAcagatcgc                                               >  5:32470/1‑141 (MQ=255)
                                                                                                        cagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTcc                                             <  7:65207/149‑1 (MQ=255)
                                                                                                              tGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAAc                                         <  6:91159/147‑1 (MQ=255)
                                                                                                                     cacaGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTc                                       >  7:88899/1‑142 (MQ=255)
                                                                                                                     cacaGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTc                                       <  8:88899/142‑1 (MQ=255)
                                                                                                                     cacaGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAt                                 >  7:253783/1‑148 (MQ=255)
                                                                                                                     cacaGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACAGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAt                                 >  7:200224/1‑148 (MQ=255)
                                                                                                                      acaGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGc                                                                       >  7:23391/1‑109 (MQ=255)
                                                                                                                      acaGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGc                                                                       <  8:23391/109‑1 (MQ=255)
                                                                                                                          gCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAggg                           >  5:247067/1‑149 (MQ=255)
                                                                                                                            acacaAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGc                          >  1:132772/1‑148 (MQ=255)
                                                                                                                              acaAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGAt                       <  7:234622/149‑1 (MQ=255)
                                                                                                                               caAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATg                      <  7:85680/149‑1 (MQ=255)
                                                                                                                               caAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATg                      <  2:132772/149‑1 (MQ=255)
                                                                                                                                          ccGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTg                                                                                  >  6:140489/1‑78 (MQ=255)
                                                                                                                                          ccGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTg                                                                                  <  5:140489/78‑1 (MQ=255)
                                                                                                                                          ccGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATg                      >  7:221999/1‑138 (MQ=255)
                                                                                                                                          ccGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATg                      <  8:221999/138‑1 (MQ=255)
                                                                                                                                                   aGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGc  <  5:175336/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGC  >  NC_002947/4486069‑4486364

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