Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 989,635 A→G 7.5% Q53R (CAG→CGG)  PP_0855 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947989,6350AG 7.5% 151.3 / 2.7 53Q53R (CAG→CGG) PP_0855hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base G (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

CCGATGATGATGACCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTGGCCAAGGT  >  NC_002947/989488‑989779
                                                                                                                                                   |                                                                                                                                                
ccGATGATGATGACCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAg                                                                                                                                                 >  4:215805/1‑149 (MQ=255)
 cGATGATGATGACCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGt                                                                                                                                                >  1:218019/1‑149 (MQ=255)
 cGATGATGATGACCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGt                                                                                                                                                >  5:56337/1‑149 (MQ=255)
  gatgatgatGACCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGtt                                                                                                                                               <  5:76897/149‑1 (MQ=255)
    tgatgatgaCCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGt                                                                                                                                             <  3:215805/149‑1 (MQ=255)
    tgatgatgaCCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGt                                                                                                                                             <  7:163997/149‑1 (MQ=255)
         atgcCCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTTGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCa                                                                                                                                         <  5:282404/148‑1 (MQ=255)
            aCCTTGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTTGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGgc                                                                                                                                      <  4:328621/148‑1 (MQ=255)
            aCCTGGCAGGCGTCAAGGACTGGTCGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGgcg                                                                                                                                     >  4:206403/1‑149 (MQ=255)
             ccTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGgcgc                                                                                                                                    <  6:72568/149‑1 (MQ=255)
             ccTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGgcgc                                                                                                                                    <  3:206403/149‑1 (MQ=255)
               tgtgaGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACGGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCa                                                                                                                                                  <  6:204236/129‑1 (MQ=255)
               tGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCa                                                                                                                                                  >  5:204236/1‑133 (MQ=255)
                 gCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGAt                                                                                                                                <  8:44827/149‑1 (MQ=255)
                    ggCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGa                                                                                                                             <  8:94809/149‑1 (MQ=255)
                        tCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAAc                                                                                                                          >  7:96150/1‑148 (MQ=255)
                                   tggAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACacc                                                                                                              >  7:8743/1‑149 (MQ=255)
                                   tggAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACacc                                                                                                              <  6:56337/149‑1 (MQ=255)
                                   tggAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACacc                                                                                                              >  8:196421/1‑149 (MQ=255)
                                       aCCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACaccgacc                                                                                                          >  4:333364/1‑149 (MQ=255)
                                         cGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGTCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCgg                                                                                                        >  8:217520/1‑149 (MQ=255)
                                          gCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGc                                                                                                       >  5:79950/1‑149 (MQ=255)
                                           cAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGc                                                                                                       >  5:53248/1‑148 (MQ=255)
                                             aCGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAg                                                                                                    <  8:96150/149‑1 (MQ=255)
                                              cGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGc                                                                                                   <  2:218019/149‑1 (MQ=255)
                                                          gctgACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGccacc                                                                                        >  5:223129/1‑148 (MQ=255)
                                                          gctgACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGccacc                                                                                        >  5:218066/1‑148 (MQ=255)
                                                           ctgACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGccacca                                                                                       <  8:8743/148‑1 (MQ=255)
                                                                  gTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGGAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTcgcg                                                                               >  5:255654/1‑149 (MQ=255)
                                                                  gTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAg                                                                                                    <  8:299013/128‑1 (MQ=255)
                                                                  gTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAg                                                                                                    >  7:299013/1‑128 (MQ=255)
                                                                  gTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTcgcg                                                                               <  7:196421/149‑1 (MQ=255)
                                                                       ttGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACt                                                                          >  7:185528/1‑149 (MQ=255)
                                                                             gctggcGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGTCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCgg                                                                    <  7:217520/149‑1 (MQ=255)
                                                                                           gtTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGc                                                       >  1:213814/1‑148 (MQ=255)
                                                                                              gggCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAg                                                    <  6:79950/148‑1 (MQ=255)
                                                                                                 cTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTg                                                <  3:333364/149‑1 (MQ=255)
                                                                                                 cTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTg                                                >  2:174453/1‑149 (MQ=255)
                                                                                                          cTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCg                                       >  5:107306/1‑149 (MQ=255)
                                                                                                               acaAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTAc                                  <  6:223129/149‑1 (MQ=255)
                                                                                                               acaAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTAc                                  <  8:250182/149‑1 (MQ=255)
                                                                                                               acaAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTAc                                  <  6:218066/149‑1 (MQ=255)
                                                                                                                caAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCGGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAg                                                                                     >  3:349972/1‑97 (MQ=255)
                                                                                                                caAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCGGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAg                                                                                     >  3:355549/1‑97 (MQ=255)
                                                                                                                caAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCGGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAg                                                                                     <  4:349972/97‑1 (MQ=255)
                                                                                                                caAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCGGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAg                                                                                     <  4:355549/97‑1 (MQ=255)
                                                                                                                   gTACCAGACCAAACAGTAGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACCCCACCAAGGTCGCGGAAATGCCCGGCAAGCTCAAGAGCGAGTATGGCGGCACCGCTTACGcc                               >  2:110764/1‑148 (MQ=255)
                                                                                                                       cAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGaa                                                                                                                            >  2:320811/1‑51 (MQ=255)
                                                                                                                       cAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGaa                                                                                                                            <  1:320811/51‑1 (MQ=255)
                                                                                                                        aGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTAc                         >  1:31333/1‑149 (MQ=255)
                                                                                                                                       aaGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTg          <  6:107306/149‑1 (MQ=255)
                                                                                                                                        aGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTgg         <  6:53248/149‑1 (MQ=255)
                                                                                                                                             ccTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTGGCCAAg    >  3:102889/1‑149 (MQ=255)
                                                                                                                                               tCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTGGCCAAGGt  >  2:108080/1‑149 (MQ=255)
                                                                                                                                               tCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTGGCCAAGGt  >  8:224757/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGATGATGATGACCTGGCAGGCGTCAAGGACTGGTGGAACCGCAACGGCAAGCCGCTGCTGACCGGTGCATTGCTGGCTGGCGTGGTGGTGTTGGGCTGGAATACCTGGCACAAGTACCAGACCAACCAGTCGCAAGGTGCCTCGCAGTTGTATCAGGCGCTGATTGAAACCAGCCTGACACCGACCGGCCAGCCGGACGCCACCAAGGTCGCGGAACTGGCCGGCAAGCTCAAGAGCGAGTTTGGCGGCACCGCTTACGCCCAGTACGGCAGCCTGTTCGTGGCCAAGGT  >  NC_002947/989488‑989779

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