Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,885,730 C→A 10.0% R50L (CGC→CTC)  PP_2541 ← transcriptional factor‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,885,7300CA10.0% 115.5 / 3.1 40R50L (CGC→CTC) PP_2541transcriptional factor‑like protein
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base A (2/2);  total (18/22)
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.51e-01

GTACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGG  >  NC_002947/2885584‑2885874
                                                                                                                                                  |                                                                                                                                                
gTACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCgcgc                                                                                                                                                <  4:319965/149‑1 (MQ=255)
 tACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCgcgc                                                                                                                                                >  8:50963/1‑148 (MQ=255)
 tACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCgcgc                                                                                                                                                <  7:50963/148‑1 (MQ=255)
   cgacgTGTGTCGGTGCGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTc                                                                                                                                              <  7:304642/148‑1 (MQ=255)
      cgTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCGCGACCTGGTCGAGCGCGCCg                                                                                                                                          <  1:409516/149‑1 (MQ=255)
                  gggTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtg                                                                                                                              <  3:301026/149‑1 (MQ=255)
                    gTAAGTCCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggt                                                                                                                            <  6:373398/149‑1 (MQ=255)
                     tAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggtg                                                                                                                           >  1:234900/1‑149 (MQ=255)
                       aGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggtgg                                                                                                                          <  5:105451/148‑1 (MQ=255)
                         tGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTTCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCgcgc                                                                                                                                                <  5:166143/124‑1 (MQ=255)
                         tGCGTCGGCTTGAAGCTGCTGCGCTTCGCCAAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGAATTCCTTTACTGTTGCAATGAAGCGGGAAAAGTGAGGCGAGCTCGACCAGGTCgcgc                                                                                                                                                >  6:166143/1‑124 (MQ=255)
                          gCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATg                                                                                                                      >  4:182736/1‑149 (MQ=255)
                           cGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGt                                                                                                                                    <  4:42731/134‑1 (MQ=255)
                               gggTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTAt                                                                                                                       >  3:330420/1‑143 (MQ=255)
                               gggTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTAt                                                                                                                       <  4:330420/143‑1 (MQ=255)
                                    gCAGCGGCGGCGCTTCGGCCAGGCTCTAGCCCGGTTGCCAGTACCAGCCATACGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGGCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGt                                                                                                             <  6:31055/148‑1 (MQ=255)
                                     cAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGTTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGtt                                                                                                            <  4:233255/148‑1 (MQ=255)
                                           ctgcGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAgg                                                                                                      <  7:31452/148‑1 (MQ=255)
                                                    gCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGc                                                                                             <  3:182736/148‑1 (MQ=255)
                                                              aaGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTgg                                                                                   <  2:234900/148‑1 (MQ=255)
                                                              aaGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGa                                                                                  <  3:305148/149‑1 (MQ=255)
                                                                          cAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAgg                                                                                                      <  8:375712/117‑1 (MQ=255)
                                                                          cAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAgg                                                                                                      >  7:375712/1‑117 (MQ=255)
                                                                          cAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCgg                                                                      >  7:295859/1‑149 (MQ=255)
                                                                           aGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCggg                                                                     <  4:263797/149‑1 (MQ=255)
                                                                               ccAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTggc                                                                 >  3:130554/1‑149 (MQ=255)
                                                                                            aTGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCaa                                                    >  8:238770/1‑149 (MQ=255)
                                                                                                   tACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAg                                              <  7:238770/148‑1 (MQ=255)
                                                                                                        ttGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGa                                               >  2:139985/1‑142 (MQ=255)
                                                                                                        ttGAACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGa                                               <  1:139985/142‑1 (MQ=255)
                                                                                                           cACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCa                                     >  5:16445/1‑149 (MQ=255)
                                                                                                              tGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTg                                  >  2:148567/1‑149 (MQ=255)
                                                                                                                          aGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGgt                      >  3:292644/1‑149 (MQ=255)
                                                                                                                              cggcgAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTgg                  >  5:70136/1‑149 (MQ=255)
                                                                                                                               gtcgAGCGCGACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGc                 >  6:98241/3‑149 (MQ=255)
                                                                                                                                   aGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggt             >  5:6964/1‑149 (MQ=255)
                                                                                                                                   aGCGCGACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGTGCAGCCTGCAACTGGGTGGAGATCATTCACGGGTGGCGGCCGGGCAGCAACTAGATCAGTTGCCACTGGCGGGCGATGGTGAGGcgggt             >  4:177608/1‑149 (MQ=255)
                                                                                                                                      tCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgg          >  4:342055/1‑149 (MQ=255)
                                                                                                                                        gACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCGGGGTGGCGGCCGGGCAGCATCTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgggt        <  3:177608/149‑1 (MQ=255)
                                                                                                                                        gACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgggt        <  6:70136/149‑1 (MQ=255)
                                                                                                                                              gTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCggg                                                          >  7:192730/1‑93 (MQ=255)
                                                                                                                                              gTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCggg                                                          <  8:192730/93‑1 (MQ=255)
                                                                                                                                              gTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTgg  <  6:6964/149‑1 (MQ=255)
                                                                                                                                              gTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTgg  <  6:16445/149‑1 (MQ=255)
                                                                                                                                              gTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCAGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTgg  >  6:128054/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGG  >  NC_002947/2885584‑2885874

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