Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 959,850 A→T 9.1% D40V (GAC→GTC)  PP_0821 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947959,8500AT 9.1% 121.3 / 3.5 44D40V (GAC→GTC) PP_0821hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/18);  new base T (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GATGCCTCCGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGG  >  NC_002947/959705‑959998
                                                                                                                                                 |                                                                                                                                                    
gATGCCTCCGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGATTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGAc                                                                                                                                                     >  1:190947/1‑147 (MQ=255)
 aTGCCTCCGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAg                                                                                                                                                  <  4:61598/149‑1 (MQ=255)
     cTCCGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCt                                                                                                                                              >  4:254471/1‑149 (MQ=255)
      tCCGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTc                                                                                                                                             >  5:78210/1‑149 (MQ=255)
        cGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTcc                                                                                                                                            <  7:233787/148‑1 (MQ=255)
                  ccAGACGCTCTGATACCCACAAGGAGGCTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGAGGc                                                                                                                                 >  4:37899/2‑149 (MQ=255)
                            tGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAgcgc                                                                                                                        <  6:78210/148‑1 (MQ=255)
                              aTACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGt                                                                                                                      <  7:19499/148‑1 (MQ=255)
                                  ccACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATc                                                                                                                 >  8:161133/1‑149 (MQ=255)
                                            ggTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGAAATGGCGGACGATTACACGCCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGaaa                                                                                                       >  5:216863/1‑149 (MQ=255)
                                              ttAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGTCAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAaca                                                                                                     >  5:33659/1‑149 (MQ=255)
                                              ttAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGTCAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAac                                                                                                      >  5:29661/1‑148 (MQ=255)
                                                 gCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACaa                                                                                                  <  3:254471/149‑1 (MQ=255)
                                                   aTGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGt                                                                                                >  3:6075/1‑149 (MQ=255)
                                                      ttCAAGGATTACCTGGCAGGGCTTGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGc                                                                                              <  6:216863/147‑1 (MQ=255)
                                                       tCAAGCATTACCTGGCAGGACTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAg                                                                                            >  2:78095/1‑149 (MQ=255)
                                                             aTTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAg                                                                                      <  7:161133/149‑1 (MQ=255)
                                                              ttACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAgc                                                                                     >  1:64663/1‑149 (MQ=255)
                                                                aCCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAACTGCAGGACCAgcgc                                                                                   >  5:210598/1‑149 (MQ=255)
                                                                 ccTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCt                                                                                  <  4:6075/149‑1 (MQ=255)
                                                                     gCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGaa                                                                              >  5:8787/1‑149 (MQ=255)
                                                                                    gCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAag                                                                 >  2:256450/1‑147 (MQ=255)
                                                                                            gCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATgg                                                       <  6:210598/149‑1 (MQ=255)
                                                                                             cTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGTCAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATgg                                                       <  6:29661/148‑1 (MQ=255)
                                                                                             cTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGTCAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATgg                                                       <  6:33659/148‑1 (MQ=255)
                                                                                                            gCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCa                                                    <  7:12324/136‑1 (MQ=255)
                                                                                                            gCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCa                                                    >  8:12324/1‑136 (MQ=255)
                                                                                                                 cGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGgcag                                   >  5:77548/1‑148 (MQ=255)
                                                                                                                  gATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCgg                                                                          <  4:72619/108‑1 (MQ=255)
                                                                                                                  gATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCgg                                                                          >  3:72619/1‑108 (MQ=255)
                                                                                                                    ttACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCaga                               <  6:77548/149‑1 (MQ=255)
                                                                                                                      acacGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCagaag                             >  3:142955/1‑149 (MQ=255)
                                                                                                                       cacGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCagaaga                            <  6:8787/149‑1 (MQ=255)
                                                                                                                         cggcggCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGc                          >  8:60521/1‑149 (MQ=255)
                                                                                                                         cggcggCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGc                          <  7:60521/149‑1 (MQ=255)
                                                                                                                          ggcggcCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCa                                                                                             >  2:251513/1‑81 (MQ=255)
                                                                                                                          ggcggcCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCa                                                                                             <  1:251513/81‑1 (MQ=255)
                                                                                                                                cTACGGGCAGTGCATGTACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGAAGAAGAGCCTGCggg                   <  3:37899/149‑1 (MQ=255)
                                                                                                                                         gTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGttagat          >  3:197226/1‑145 (MQ=255)
                                                                                                                                           gCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAgcg         >  1:196629/1‑148 (MQ=255)
                                                                                                                                              tGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAg     <  2:196629/149‑1 (MQ=255)
                                                                                                                                               ggACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCt                                                            <  2:205669/93‑1 (MQ=255)
                                                                                                                                               ggACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCt                                                            >  1:205669/1‑93 (MQ=255)
                                                                                                                                               ggACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGt    >  1:36121/1‑149 (MQ=255)
                                                                                                                                                 aCAAGGCCTCCAGCACCGTGGCCATGAGCGCTTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGtgg  <  2:190947/149‑1 (MQ=255)
                                                                                                                                                 aCAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGtgg  <  1:256450/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATGCCTCCGGAGCAGCCGCAGACGCTTTGATACCCACAAGGAGGGTTAGCATGATCAAGCATTACCTGGCAGGCCTGGCACTGGCCTGCGTGCTACCGGTAGCAATGGCGGACGATTACACGGCGGCCTACGGGCAGTGCATGGACAAGGCCTCCAGCACCGTGGCCATGAGCGCGTGCATCGGGGCCGAAACACAAGTGCAGGACCAGCGCTTGAACCGGGTGTACAAGCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGG  >  NC_002947/959705‑959998

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