Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,405,314 C→T 9.1% G83R (GGG→AGG)  PP_3886 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,405,3140CT 9.1% 118.2 / 2.3 44G83R (GGG→AGG) PP_3886hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base T (2/2);  total (22/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 = 1.62e-01

ACAACCGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGG  >  NC_002947/4405168‑4405462
                                                                                                                                                  |                                                                                                                                                     
aCAACCGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCg                                                                                                                                                     <  2:279577/149‑1 (MQ=255)
    ccGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCAt                                                                                                                                                 >  8:12094/1‑149 (MQ=255)
           ccAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGa                                                                                                                                          <  7:65650/149‑1 (MQ=255)
            cAGCGATCTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCTTGACGTAATTTGGAGGGACTGGGCATAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGTCATGTGCCGa                                                                                                                                          <  7:121851/148‑1 (MQ=255)
                gAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCtt                                                                                                                                    <  5:268768/149‑1 (MQ=255)
                gAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCtt                                                                                                                                    <  7:12094/149‑1 (MQ=255)
                      ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGG‑GCtttttcg                                                                                                                               <  2:179499/148‑6 (MQ=255)
                       cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACTGAATTTGGAGGGACTGGGCTGAAATCTGAGCATCTAACGTCTCGACGCTTTCGATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATgg                                                                                                                             <  5:252214/149‑1 (MQ=255)
                       cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATgg                                                                                                                             >  7:37541/1‑149 (MQ=255)
                       cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATgg                                                                                                                             <  2:72538/149‑1 (MQ=255)
                          ggAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTg                                                                                                                          <  6:66997/149‑1 (MQ=255)
                               aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGAt                                                                                                                     <  6:37263/149‑1 (MQ=255)
                               aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGAt                                                                                                                     >  3:286273/1‑149 (MQ=255)
                                gggCGACCCGACTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATg                                                                                                                    <  2:159251/149‑1 (MQ=255)
                                 ggCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGTGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGG‑Gc                                                                                                                                      >  6:330913/1‑130 (MQ=255)
                                 ggCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGTGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGG‑Gc                                                                                                                                      <  5:330913/130‑1 (MQ=255)
                                    gACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCaa                                                                                                                >  8:9514/1‑149 (MQ=255)
                                     aCCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAAc                                                                                                               <  7:9514/149‑1 (MQ=255)
                                        cGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCt                                                                                                            >  3:300980/1‑149 (MQ=255)
                                                  cgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGc                                                                                                  <  3:114096/149‑1 (MQ=255)
                                                  cgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGc                                                                                                  >  4:114096/1‑149 (MQ=255)
                                                       aCCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGTGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTc                                                                                              <  1:61854/149‑1 (MQ=255)
                                                        cccGTGACGGAATTTGCAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGt                                                                                            >  1:302510/1‑149 (MQ=255)
                                                        accGTGACGGTATTTGGAGGGACTGGGCAGAAATCTTAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGt                                                                                            <  4:286273/148‑1 (MQ=255)
                                                           gTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgc                                                                                         >  2:3387/1‑149 (MQ=255)
                                                                ggAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCACCCCTGTCACGTCGCTGTCGAGGCGCTTc                                                                                    >  8:237639/1‑149 (MQ=255)
                                                                ggAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                    >  2:24787/1‑149 (MQ=255)
                                                                ggAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCATGGAATGGGTGGAGATGGCAAACCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                    >  5:69736/1‑149 (MQ=255)
                                                                    tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                <  8:37541/149‑1 (MQ=255)
                                                                                           tGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTAcc                                                         >  2:57223/1‑149 (MQ=255)
                                                                                                    aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTc                                                                           <  8:177239/122‑1 (MQ=255)
                                                                                                    aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTc                                                                           >  7:177239/1‑122 (MQ=255)
                                                                                                    aaCGTCTCGAAGCTATCCATGACTGGCTACATTTAACGAAAAATCACCGGCATGTGCCGAG‑GCTTGGAATGTGTGGAGATGGAAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCATCCAGTCCCCACGCTGACTGTTAAACCTGCACGt                                                >  6:181725/1‑149 (MQ=255)
                                                                                                     aCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGtt                                               >  3:255281/1‑149 (MQ=255)
                                                                                                                   tCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                                  <  1:24787/148‑1 (MQ=255)
                                                                                                                      aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatc                              <  4:300980/149‑1 (MQ=255)
                                                                                                                      aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                               >  6:266208/1‑148 (MQ=255)
                                                                                                                        gACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatcat                            >  5:44699/1‑149 (MQ=255)
                                                                                                                               tCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaa                      >  1:63715/1‑148 (MQ=255)
                                                                                                                               tCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaa                      >  4:100498/1‑148 (MQ=255)
                                                                                                                                ccATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaaa                     >  2:204599/1‑148 (MQ=255)
                                                                                                                                  aTTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAAc                  <  4:255281/149‑1 (MQ=255)
                                                                                                                                         gAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg           <  2:63715/149‑1 (MQ=255)
                                                                                                                                             aaGCCTCGGCACATGCCGGG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGcc         >  4:161714/1‑147 (MQ=255)
                                                                                                                                                 cccGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCGGcca   >  2:267606/1‑147 (MQ=255)
                                                                                                                                                 cccGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCACCAAAAACGACAAAGCCCACCTg   >  7:156451/1‑149 (MQ=255)
                                                                                                                                                  ccGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTgg  <  1:57223/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACAACCGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGG  >  NC_002947/4405168‑4405462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑