Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,869,612 T→C 7.1% P40P (CCT→CCC PP_5144 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,869,6120TC 7.1% 163.4 / 2.6 56P40P (CCT→CCCPP_5144hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base C (2/2);  total (30/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 = 2.98e-01

ACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCT  >  NC_002947/5869468‑5869758
                                                                                                                                                 |                                                                                                                                                  
aCCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTggg                                                                                                                                                 <  2:266939/148‑1 (MQ=255)
  cAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTg                                                                                                                                              <  4:116498/149‑1 (MQ=255)
      aaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCCAGGTGCAtc                                                                                                                                          <  4:203422/149‑2 (MQ=255)
        acacAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGcc                                                                                                                                        <  5:153230/149‑1 (MQ=255)
          acaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCa                                                                                                                                            <  8:46043/143‑1 (MQ=255)
          acaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCa                                                                                                                                            >  7:46043/1‑143 (MQ=255)
          acaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCCAGGTGCAtcctcc                                                                                                                                      <  1:244320/149‑6 (MQ=255)
           caGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCggg                                                                                                                                     >  1:216608/1‑149 (MQ=255)
             gagGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGgcg                                                                                                                                   <  7:109229/149‑1 (MQ=255)
                 cccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAc                                                                                                                               <  3:168541/149‑1 (MQ=255)
                  ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                              >  5:68956/1‑149 (MQ=255)
                       cTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAc                                                                                                                               <  8:212603/143‑1 (MQ=255)
                       cTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAc                                                                                                                               >  7:212603/1‑143 (MQ=255)
                              ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                  <  5:104255/149‑1 (MQ=255)
                                  aTCGTATTTGCGTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCGGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTa                                                                                                               <  2:66453/148‑1 (MQ=255)
                                       aTTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggc                                                                                                         >  1:145855/1‑149 (MQ=255)
                                           gCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTggg                                                                                                     >  3:35184/1‑149 (MQ=255)
                                            ccTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTggg                                                                                                     <  2:216608/148‑1 (MQ=255)
                                             cTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                              >  4:119888/1‑122 (MQ=255)
                                                    ccGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGTATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGccc                                                                                            <  6:118201/149‑1 (MQ=255)
                                                     cGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCca                                                                                           >  2:227836/1‑149 (MQ=255)
                                                                     gCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCcg                                                                           >  2:249660/1‑149 (MQ=255)
                                                                        gATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCgcgc                                                                        >  4:47285/1‑149 (MQ=255)
                                                                         aTCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCt                                                                       >  2:240653/1‑149 (MQ=255)
                                                                             tGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGAcc                                                                   <  4:35184/149‑1 (MQ=255)
                                                                                  gcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACAt                                                              <  1:227836/149‑1 (MQ=255)
                                                                                   ccAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATc                                                             >  8:105336/1‑149 (MQ=255)
                                                                                   ccAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATc                                                             <  7:105336/149‑1 (MQ=255)
                                                                                            acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCg                                                    >  8:127907/1‑149 (MQ=255)
                                                                                              gaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCggg                                                  >  3:203662/1‑149 (MQ=255)
                                                                                                    tATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAgg                                             >  1:48146/1‑149 (MQ=255)
                                                                                                        cccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGTGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAg                                        >  6:159327/1‑149 (MQ=255)
                                                                                                        cccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAg                                        >  5:1461/1‑149 (MQ=255)
                                                                                                         ccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                       >  8:184812/1‑149 (MQ=255)
                                                                                                               cACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctc                                 <  4:203662/149‑1 (MQ=255)
                                                                                                                    aGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctcgctc                             <  2:48146/148‑1 (MQ=255)
                                                                                                                      ccccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTcgcg                          >  8:242583/1‑149 (MQ=255)
                                                                                                                       cccTGGCGGCGTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCggg                                                  <  7:217442/124‑1 (MQ=255)
                                                                                                                       cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCTCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCgg                         >  4:58963/1‑149 (MQ=255)
                                                                                                                       cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCggg                                                  >  8:217442/1‑124 (MQ=255)
                                                                                                                        ccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctc                                 <  4:133451/140‑1 (MQ=255)
                                                                                                                        ccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctc                                 >  3:133451/1‑140 (MQ=255)
                                                                                                                            gcggcCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGa                    <  3:47285/149‑1 (MQ=255)
                                                                                                                               gcCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCt                 <  1:240653/149‑1 (MQ=255)
                                                                                                                                  tGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggt              >  2:41797/1‑149 (MQ=255)
                                                                                                                                    ggaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtgg            <  1:41797/149‑1 (MQ=255)
                                                                                                                                     gaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctc                                 <  6:113054/127‑1 (MQ=255)
                                                                                                                                     gaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctc                                 >  5:113054/1‑127 (MQ=255)
                                                                                                                                     ccggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGTGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtgg            >  2:11243/6‑149 (MQ=255)
                                                                                                                                        gcaTGCACCCCAGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAg         >  8:52338/3‑149 (MQ=255)
                                                                                                                                         gaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGc        >  7:43423/1‑149 (MQ=255)
                                                                                                                                              aCCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGttt    >  1:91558/1‑148 (MQ=255)
                                                                                                                                              aCCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGttt    <  2:91558/148‑1 (MQ=255)
                                                                                                                                              aCCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTc   >  8:40235/1‑149 (MQ=255)
                                                                                                                                              aCCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGCTGGCAGCGTTTc   >  2:141329/1‑149 (MQ=255)
                                                                                                                                               ccTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCt  <  7:242583/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCT  >  NC_002947/5869468‑5869758

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: