Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,257,480 G→C 9.4% A158G (GCC→GGC)  PP_2854 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,257,4800GC 9.4% 151.6 / 7.4 53A158G (GCC→GGC) PP_2854hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base C (3/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

CGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGCCGTTG  >  NC_002947/3257332‑3257623
                                                                                                                                                     |                                                                                                                                               
cGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGGTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcgg                                                                                                                                                 <  7:110618/149‑1 (MQ=255)
      gTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                           <  2:210642/149‑1 (MQ=255)
      gTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                           <  4:324354/149‑1 (MQ=255)
         cTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                               >  1:217997/1‑143 (MQ=255)
         cTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                               <  2:217997/142‑1 (MQ=255)
            cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAgaccc                                                                                                                                      <  6:517/148‑5 (MQ=255)
            cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGCCGGCCAgacc                                                                                                                                       <  7:214786/147‑4 (MQ=255)
            cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGATGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCt                                                                                                                                     <  3:73048/149‑1 (MQ=255)
                agGCCACCTTGCGCCAGGCCTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACACACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGcc                                                                                                                                 >  4:242980/1‑149 (MQ=255)
                       cTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgt                                                                                                                          >  6:167798/1‑149 (MQ=255)
                          gggCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                         <  3:131407/147‑1 (MQ=255)
                          gggCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                         <  2:160514/147‑1 (MQ=255)
                                ggCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCCGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTAcggcg                                                                                                                 >  3:284386/1‑149 (MQ=255)
                                      gCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGcc                                                                                                           >  8:173553/1‑149 (MQ=255)
                                           gctgctGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGt                                                                                                       <  7:197278/149‑1 (MQ=255)
                                            ctgctgCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGtt                                                                                                      <  1:183973/149‑1 (MQ=255)
                                               ctgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                                    <  6:137292/148‑1 (MQ=255)
                                                   ccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                               <  4:327247/149‑1 (MQ=255)
                                                    cTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTgg                                                                                              >  8:281557/1‑149 (MQ=255)
                                                        ttgttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCt                                                                                                          <  2:302398/133‑1 (MQ=255)
                                                        ttgttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCt                                                                                                          >  1:302398/1‑133 (MQ=255)
                                                           ttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTg                                                                                                         >  7:27032/1‑131 (MQ=255)
                                                           ttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTg                                                                                                         <  8:27032/131‑1 (MQ=255)
                                                             gCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                               <  1:135220/139‑1 (MQ=255)
                                                                 gCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATc                                                                                 >  6:161264/1‑149 (MQ=255)
                                                                       aTGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgtt                                                                           >  1:157334/1‑149 (MQ=255)
                                                                        tGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttg                                                                          <  5:167798/149‑1 (MQ=255)
                                                                          ttGAGCTGCACACCCTCACCGCTGGTAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                        <  8:314098/149‑1 (MQ=255)
                                                                          ttGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCTCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                        <  8:309834/149‑1 (MQ=255)
                                                                          ttGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                        >  6:208238/1‑149 (MQ=255)
                                                                             aGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTtgttgt                                                                      <  3:245631/148‑1 (MQ=255)
                                                                                 gCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                                    >  7:59225/1‑114 (MQ=255)
                                                                                 gCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcg                                                                 <  4:284386/149‑1 (MQ=255)
                                                                                   acacccacacgGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcggc                                                               <  7:164784/149‑1 (MQ=255)
                                                                                       ccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgttg                                                           >  6:263471/1‑149 (MQ=255)
                                                                                             cGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCa                                                     <  7:281557/149‑1 (MQ=255)
                                                                                              gCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCAc                                                    <  7:173553/149‑1 (MQ=255)
                                                                                                        acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACg                                          >  3:218727/1‑149 (MQ=255)
                                                                                                        acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACg                                          >  8:191567/1‑149 (MQ=255)
                                                                                                          acTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGa                                        >  6:127982/1‑149 (MQ=255)
                                                                                                            tGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGac                                       >  3:152968/1‑148 (MQ=255)
                                                                                                               gcgcgcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGaccacc                                   >  8:11354/1‑149 (MQ=255)
                                                                                                                 gcgcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACTTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTg                                 >  4:239061/1‑149 (MQ=255)
                                                                                                                      ggCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTa                            <  5:263471/149‑1 (MQ=255)
                                                                                                                       gcgcCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                        >  8:116259/1‑104 (MQ=255)
                                                                                                                         gcCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                         >  3:265964/1‑149 (MQ=255)
                                                                                                                                tCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCgg                   <  3:242980/148‑1 (MQ=255)
                                                                                                                                 ccccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggt                  <  4:218727/148‑1 (MQ=255)
                                                                                                                                   ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggt               >  8:84729/1‑149 (MQ=255)
                                                                                                                                     ggccAGGTACTGGCCGCCGGCCAGACCCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtgg             >  4:112697/3‑149 (MQ=255)
                                                                                                                                      aCCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtggt            <  5:208238/149‑1 (MQ=255)
                                                                                                                                        caggtacTGGCCGCCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggt                  >  3:317868/7‑141 (MQ=255)
                                                                                                                                        caggtacTGGCCGCCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggt                  <  4:317868/135‑1 (MQ=255)
                                                                                                                                         aggtacTGGCCGCCGGCCAGACCCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAg         >  4:204545/6‑149 (MQ=255)
                                                                                                                                          gggTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAgg        <  5:161264/149‑1 (MQ=255)
                                                                                                                                                ggCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGCCGTTg  >  7:213629/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGCCGTTG  >  NC_002947/3257332‑3257623

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: