Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,558,076 C→T 8.1% G27D (GGC→GAC)  PP_4045 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,558,0760CT 8.1% 138.5 / 2.3 49G27D (GGC→GAC) PP_4045hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base T (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGAATGTCCTCGACA  >  NC_002947/4557930‑4558217
                                                                                                                                                  |                                                                                                                                             
cgccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCg                                                                                                                                             >  1:117505/1‑149 (MQ=255)
  ccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgc                                                                                                                                            <  6:334195/148‑1 (MQ=255)
  ccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgc                                                                                                                                            >  5:334195/1‑148 (MQ=255)
     tGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgctgct                                                                                                                                        <  8:143022/149‑1 (MQ=255)
     tGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgctgct                                                                                                                                        >  3:59674/1‑149 (MQ=255)
           gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                                                                                                                  <  5:128824/149‑1 (MQ=255)
           gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                                                                                                                  >  1:61226/1‑149 (MQ=255)
            tAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgc                                                                                                                                            >  8:327593/1‑138 (MQ=255)
            tAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgc                                                                                                                                            <  7:327593/138‑1 (MQ=255)
            tAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtg                                                                                                                                 >  7:124975/1‑149 (MQ=255)
              gCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggt                                                                                                                               <  2:117505/149‑1 (MQ=255)
              gCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggt                                                                                                                               >  7:68727/1‑149 (MQ=255)
               cATCGCCTACCCTAGCAAGAGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTTATCGCCGCCCCACAGCCGAGCCTGTCGCCAAGCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggtg                                                                                                                              <  4:298126/149‑1 (MQ=255)
                   gCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTgggg                                                                                                                           <  2:61226/148‑1 (MQ=255)
                         ccTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTgggg                                                                                                                           >  2:306573/1‑142 (MQ=255)
                         ccTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTgggg                                                                                                                           <  1:306573/142‑1 (MQ=255)
                           tAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGTCGCTGCTGACGCTGGTGGGGCCGCACAAt                                                                                                                  >  4:279424/1‑149 (MQ=255)
                             gCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATg                                                                                                                 <  8:124975/148‑1 (MQ=255)
                             gCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATg                                                                                                                 <  4:59674/148‑1 (MQ=255)
                                     cTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGTCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGt                                                                                                        <  3:279424/149‑1 (MQ=255)
                                          cTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAgg                                                                                                    <  1:299882/148‑1 (MQ=255)
                                          cTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAgg                                                                                                    <  8:68727/148‑1 (MQ=255)
                                             catccatcATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGa                                                                                                >  6:93911/1‑149 (MQ=255)
                                              atccatcATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATg                                                                                                                 <  6:135287/131‑1 (MQ=255)
                                              atccatcATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATg                                                                                                                 >  5:135287/1‑131 (MQ=255)
                                                  atcatcAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATg                                                                                                          <  4:73335/134‑1 (MQ=255)
                                                  atcatcAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATg                                                                                                          >  3:73335/1‑134 (MQ=255)
                                                        aGCGCCGTGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACTAAGAACAGTCTGGAGCTGCTGTTGACGCTGTTGGGGCAGCACAATGGGCAATGGTGAGGGTGACCAACCAGaa                                                                                      <  5:246670/148‑1 (MQ=255)
                                                                                           gAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGt                                                   >  2:159471/1‑148 (MQ=255)
                                                                                           gAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGt                                                   >  6:49075/1‑148 (MQ=255)
                                                                                             gCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTaaa                                                <  5:49075/149‑1 (MQ=255)
                                                                                               cTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAgc                                              <  5:33975/149‑1 (MQ=255)
                                                                                                 gTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAgcgc                                            >  8:296720/1‑149 (MQ=255)
                                                                                                       aaCCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCa                                      >  6:299199/1‑149 (MQ=255)
                                                                                                         cccTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCAc                                     >  8:10955/1‑148 (MQ=255)
                                                                                                             ggcggcCGCCCTGCCGTTGACGAAGAACAGGCTGGAGTCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATg                                >  3:308281/1‑149 (MQ=255)
                                                                                                                ggcCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCg                             >  6:116803/1‑149 (MQ=255)
                                                                                                                            gTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAAt                  >  1:209852/1‑148 (MQ=255)
                                                                                                                             ttGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCa                <  5:299199/149‑1 (MQ=255)
                                                                                                                             ttGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCa                <  7:10955/149‑1 (MQ=255)
                                                                                                                               gACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAg               >  8:267996/1‑148 (MQ=255)
                                                                                                                                aCGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGa              <  5:93911/148‑1 (MQ=255)
                                                                                                                                  gaagaaCAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGAATg           >  2:294941/1‑149 (MQ=255)
                                                                                                                                  gaagaaCAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGAATg           >  2:108370/1‑149 (MQ=255)
                                                                                                                                      aaCAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGAATGTCCt       <  5:116803/149‑1 (MQ=255)
                                                                                                                                         aGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAACCAGAATGTCCTCGa    <  1:159471/149‑1 (MQ=255)
                                                                                                                                           gCTGGAGTCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGAATGTCCTCGACa  <  4:308281/149‑1 (MQ=255)
                                                                                                                                           gCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAgcg                                             >  6:118542/1‑106 (MQ=255)
                                                                                                                                           gCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAgcg                                             <  5:118542/106‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAGGCATTGTCGCCACTCCTGAATCATGTAAAGCGCCACGCACTCATGCCGGCAGCGCCAATCAGAATGTCCTCGACA  >  NC_002947/4557930‑4558217

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