Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,623,447 A→C 7.8% P114P (CCT→CCG)  PP_4941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,623,4470AC 7.8% 138.1 / 4.2 51P114P (CCT→CCG) PP_4941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/28);  new base C (2/2);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAA  >  NC_002947/5623300‑5623584
                                                                                                                                                   |                                                                                                                                         
cTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGACCTGTTCGCCAGCCAGCAGCTCGACCCGAgg                                                                                                                                         >  2:1041161/1‑150 (MQ=255)
     gCCGGCATTGCCGGCTGCAGGGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga                                                                                                                                    <  1:753995/150‑1 (MQ=255)
     gCCGGCATTGCCGGCTGCAGGGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga                                                                                                                                    <  1:753025/150‑1 (MQ=255)
     gCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga                                                                                                                                    <  1:8236/150‑1 (MQ=255)
     gCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga                                                                                                                                    <  1:567602/150‑1 (MQ=255)
      ccGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgag                                                                                                                                   <  2:1095003/150‑1 (MQ=255)
        ggCAATGCCGGCCGCGCAGTCGTAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgag                                                                                                                                   >  2:658562/1‑148 (MQ=255)
        ggCAATGCCGGCCGCGCAGTCGTAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgag                                                                                                                                   <  1:658562/148‑1 (MQ=255)
                        cAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAgct                                                                                                                 <  1:278651/150‑1 (MQ=255)
                         aGTCGCAACGTACCGGCCTGCAGTGGCGCAGCCAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAgctg                                                                                                                >  2:1042987/1‑150 (MQ=255)
                           tCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAgctgct                                                                                                              >  1:1019168/1‑150 (MQ=255)
                                aCGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTa                                                                                                         <  2:1019168/150‑1 (MQ=255)
                                  gTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGc                                                                                                       <  2:966957/150‑1 (MQ=255)
                                  gTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGc                                                                                                       <  1:1202577/150‑1 (MQ=255)
                                       ggCCTGCAGTGGCGCAGCAAGACCCGTCGGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAg                                                                                                  >  1:1163205/1‑150 (MQ=255)
                                        gCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGTAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGt                                                                                                 >  2:76828/1‑150 (MQ=255)
                                                    gcagcaAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAgg                                                                                     <  1:1005185/150‑1 (MQ=255)
                                                    gcagcaAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAgg                                                                                     <  1:580637/150‑1 (MQ=255)
                                                    gcaccaAGACCCGTCATGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAgg                                                                                     <  1:1104998/150‑1 (MQ=255)
                                                      agcaAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTa                                                                                                         <  1:1011710/128‑1 (MQ=255)
                                                      agcaAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTa                                                                                                         >  2:1011710/1‑128 (MQ=255)
                                                         aaGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGCTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTtc                                                                                <  2:783325/150‑1 (MQ=255)
                                                          aGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTc                                                                                                                                      <  2:895429/95‑1 (MQ=255)
                                                          aGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTc                                                                                                                                      >  1:895429/1‑95 (MQ=255)
                                                              ccGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGc                                                                           <  1:1042987/150‑1 (MQ=255)
                                                                 tCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCa                                                                                                                    <  1:775275/106‑1 (MQ=255)
                                                                 tCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCa                                                                                                                    >  2:775275/1‑106 (MQ=255)
                                                                      ttAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAAc                                                                   >  2:887951/1‑150 (MQ=255)
                                                                      ttAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAAc                                                                   <  1:887951/150‑1 (MQ=255)
                                                                          ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGc                                                               >  2:81901/1‑150 (MQ=255)
                                                                            aaGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCaa                                                             <  1:81901/150‑1 (MQ=255)
                                                                              gcgcACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTGTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATc                                                           >  2:810068/1‑150 (MQ=255)
                                                                                       aTCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGACGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGag                                                  >  1:1186998/1‑150 (MQ=255)
                                                                                          gACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagcca                                               <  2:636942/150‑1 (MQ=255)
                                                                                           aCGCCACGCTCGAACGCTTGCCGATCGACCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagccag                                              <  1:1041161/150‑1 (MQ=255)
                                                                                                aCGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATc                                         >  2:1137259/1‑150 (MQ=255)
                                                                                                aCGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCACCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAAGCAGCCTTGAGCCAGCCATc                                         <  1:211038/150‑1 (MQ=255)
                                                                                                      gAACGCTTGCCGATCGGCCTGTTCTCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAg                                   >  2:965922/1‑150 (MQ=255)
                                                                                                             tGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTggaa                            <  1:1137259/150‑3 (MQ=255)
                                                                                                                               ccagccagCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGt                                       >  1:815912/1‑121 (MQ=255)
                                                                                                                               ccagccagCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGt                                       <  2:815912/121‑1 (MQ=255)
                                                                                                                                cagccagcAGCTCGACCCGAGGGTCGAGAAAGCGACCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTg                                                                       >  2:667953/1‑88 (MQ=255)
                                                                                                                                cagccagcAGCTCGACCCGAGGGTCGAGAAAGCGACCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTg                                                                       <  1:667953/88‑1 (MQ=255)
                                                                                                                                    cagcagCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGcc     <  1:810068/150‑1 (MQ=255)
                                                                                                                                     agcagcTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAAc                       <  2:736676/131‑1 (MQ=255)
                                                                                                                                     agcagcTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAAc                       >  1:736676/1‑131 (MQ=255)
                                                                                                                                     agcagcTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCg    <  2:613778/150‑1 (MQ=255)
                                                                                                                                     agcagcTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCg    >  2:833363/1‑150 (MQ=255)
                                                                                                                                       tttcTCGACCCTCGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGaa  >  1:230173/4‑150 (MQ=255)
                                                                                                                                             gACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACctgctg                 >  2:437543/1‑129 (MQ=255)
                                                                                                                                             gACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACctgctg                 <  1:437543/129‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
CTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCAGCTGCTTTTTAGCCGCAGTTCAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAA  >  NC_002947/5623300‑5623584

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