Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,101,079 T→C 8.0% 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 8.0% 133.0 / 3.2 51V127A (GTC→GCC) PP_2710hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/21);  new base C (2/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

CGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCA  >  NC_002947/3100942‑3101223
                                                                                                                                         |                                                                                                                                                
cGGCCACGAGCGGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGACCTGgttgccac                                                                                                                                       <  3:20157/149‑5 (MQ=255)
cGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTc                                                                                                                                        <  4:317090/148‑1 (MQ=255)
 ggCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCgg                                                                                                                                      >  8:226776/1‑149 (MQ=255)
   ccACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCa                                                                                                                                    <  3:74431/149‑1 (MQ=255)
    cACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCAt                                                                                                                                   >  5:239192/1‑149 (MQ=255)
           aGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAAc                                                                                                                            <  5:216512/149‑1 (MQ=255)
                    tCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCg                                                                                                                     >  6:97705/1‑147 (MQ=255)
                     cAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGTGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAg                                                                                                                  >  6:56587/1‑149 (MQ=255)
                           cATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCa                                                                                                            <  1:358763/149‑1 (MQ=255)
                              cAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCa                                                                                                         >  2:145932/1‑149 (MQ=255)
                                gagCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgttgtt                                                                                                                                         >  3:315628/1‑115 (MQ=255)
                                gagCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGgttgtt                                                                                                                                         <  4:315628/115‑1 (MQ=255)
                                  gCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTTCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGc                                                                                                     >  7:88262/1‑149 (MQ=255)
                                        ggcggTGACTGGGAGGCCGTGTGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGc                                                                                               <  5:56587/149‑1 (MQ=255)
                                           ggTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCa                                                                                                               >  1:9494/1‑130 (MQ=255)
                                           ggTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCa                                                                                                               <  2:9494/130‑1 (MQ=255)
                                                 tGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAggc                                                                                      >  8:82716/1‑149 (MQ=255)
                                                   ggAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGACGGGTTGTTCGGCATGGTCAACGTGATCGAAGGCACCAGCATCGCCCTGGCCACCCAgccgg                                                                                    >  5:234033/1‑149 (MQ=255)
                                                            gcgcGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGACCGGgttgcca                                                                                                                                        <  6:271781/88‑4 (MQ=255)
                                                            gcgcGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGACCGGgttgcca                                                                                                                                        >  5:271781/1‑85 (MQ=255)
                                                              gcgATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCAGCGCCCTGGCCACCCAGGCGGCGGGCGTTATc                                                                         <  5:339942/149‑1 (MQ=255)
                                                                    gCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGaa                                                                    <  6:239192/148‑1 (MQ=255)
                                                                      cAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAaca                                                                  >  6:163892/1‑148 (MQ=255)
                                                                      cAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAaca                                                                  >  7:143434/1‑148 (MQ=255)
                                                                        gCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACAt                                                               >  5:354253/1‑149 (MQ=255)
                                                                          cGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTg                                                             <  5:97705/149‑1 (MQ=255)
                                                                             tGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGcc                                                          >  2:208729/1‑149 (MQ=255)
                                                                              gCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGcggcgg                                                                                 >  2:272412/1‑125 (MQ=255)
                                                                              gCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGcggcgg                                                                                 <  1:272412/125‑1 (MQ=255)
                                                                                 aTCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGGCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgc                                                      >  2:59806/1‑149 (MQ=255)
                                                                                 aTCCATCGAGCTGATGGTCGCCTACCTGTACGACCTAATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgc                                                      >  2:178580/1‑149 (MQ=255)
                                                                                    cATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgcct                                                    <  7:226776/148‑1 (MQ=255)
                                                                                      tCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGAc                                                 >  3:26569/1‑149 (MQ=255)
                                                                                        gAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCa                                               <  8:143434/149‑1 (MQ=255)
                                                                                              aTGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGtt                                         <  7:82716/149‑1 (MQ=255)
                                                                                                ggTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCggcgggc                                                                               >  4:121388/1‑109 (MQ=255)
                                                                                                ggTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCggcgggc                                                                               <  3:121388/109‑1 (MQ=255)
                                                                                                      cTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgcct                                                    >  2:227105/1‑130 (MQ=255)
                                                                                                      cTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCcctgcct                                                    <  1:227105/130‑1 (MQ=255)
                                                                                                        aCCTGTACGACCTGATCGCGCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCggcggg                                                                                >  8:93404/1‑100 (MQ=255)
                                                                                                         ccTGTACGACCTAATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTc                               <  1:178580/148‑1 (MQ=255)
                                                                                                          cTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAg                             >  7:8047/1‑149 (MQ=255)
                                                                                                                 aCCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGa                                                                     >  7:315611/1‑102 (MQ=255)
                                                                                                                 aCCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGa                                                                     <  8:315611/102‑1 (MQ=255)
                                                                                                                 aCCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGc                            >  7:129046/1‑143 (MQ=255)
                                                                                                                 aCCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGc                            <  8:129046/143‑1 (MQ=255)
                                                                                                                              gaacAACCCGACCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGAc          >  3:345986/4‑148 (MQ=255)
                                                                                                                                gCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCa        <  8:88262/148‑1 (MQ=255)
                                                                                                                                   aCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGt                   >  2:156876/1‑134 (MQ=255)
                                                                                                                                   aCCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGt                   <  1:156876/134‑1 (MQ=255)
                                                                                                                                     ccGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGca  >  4:37665/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                
CGGCCACGAGCAGTGGATCCTCAACGACATCAGAGCGTGCGGCGGTGACTGGGAGGCCGTGCGCGATGGCCAGCCGGTGCAATCCATCGAGCTGATGGTCGCCTACCTGTACGACCTGATCGCCCGTGGCAACCCGGTCGGGTTGTTCGGCATGGTCAACGTGCTCGAAGGCACCAGCATCGCCCTGGCCACCCAGGCGGCGGGCGTTATCCAGAACACATTGGCCCTGCCTGACCAGGCGTTCAGCTACCTCAGCTCGCACGGTGCGCTCGACCAGGAGCA  >  NC_002947/3100942‑3101223

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