Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,604 G→T 10.0% intergenic (+59/‑237) PP_2661 → / → PP_2662 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,6040GT10.0% 142.7 / 5.9 50intergenic (+59/‑237)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/28);  new base T (3/2);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAG  >  NC_002947/3048457‑3048749
                                                                                                                                                   |                                                                                                                                                 
tggtggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGa                                                                                                                                                  <  1:148927/149‑1 (MQ=255)
  gtggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGc                                                                                                                                                >  7:195253/1‑149 (MQ=255)
  gtggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGc                                                                                                                                                <  6:104640/149‑1 (MQ=255)
        gtgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGa                                                                                                                                          <  1:136128/149‑1 (MQ=255)
         tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCcgcg                                                                                                                                           >  6:304206/1‑147 (MQ=255)
         tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTGTGCAGCCCTTCGCGGGCGAGCCCGCGaa                                                                                                                                         <  7:86795/149‑1 (MQ=255)
          gCAGGCACGGCTGGCCTTGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCTGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTAAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGaa                                                                                                                                         <  6:74574/148‑1 (MQ=255)
           cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAgg                                                                                                                                       <  8:2100/149‑1 (MQ=255)
           cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAgg                                                                                                                                       <  8:195253/149‑1 (MQ=255)
           cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAgg                                                                                                                                       >  5:367343/1‑149 (MQ=39)
              gCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGcc                                                                                                                                    >  1:112665/1‑149 (MQ=255)
              gCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGcc                                                                                                                                    >  3:119680/1‑149 (MQ=255)
                aCGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGctg                                                                                                                                   >  4:63618/1‑146 (MQ=38)
                   gCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCGTCGCGGGCGAGCCCGCGAAGGGCCGCAAg                                                                                                                               <  7:41009/149‑1 (MQ=255)
                       gcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggc                                                                                                                             <  6:88091/147‑1 (MQ=255)
                        cctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcggc                                                                                                                          <  4:119680/149‑1 (MQ=255)
                         ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcc                                                                                                                         >  3:36273/1‑149 (MQ=255)
                         ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCAGcc                                                                                                                         <  3:63618/149‑1 (MQ=37)
                            gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCg                                                                                                                      >  6:64805/1‑149 (MQ=255)
                               gcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCAGccccaatt                                                                                                                   <  6:367343/149‑5 (MQ=25)
                                   tcgcaatgTTATGAAATGCACGCTTGCCAGCGTTCGCAGGCATTTGCTTCCTGAAGCCAGGAGATTCCGCTCATATAATGTAATTGGGGCTGGTTTGCAGCCCTTCGCGGGCGAGCCGGCGAAGGGCCGCAAGGCGTCCCCGGttcatt                                                                                                               <  8:154304/142‑1 (MQ=255)
                                      cACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCATTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCCGCCCCGGttcattca                                                                                                             >  8:77168/1‑148 (MQ=255)
                                        ccGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATc                                                                                                          <  8:275695/149‑1 (MQ=255)
                                           ttATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                                                                                                       <  7:203823/149‑1 (MQ=255)
                                           ttATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                                                                                                       <  3:118072/149‑1 (MQ=255)
                                              tGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                                                                                     >  7:29107/1‑148 (MQ=255)
                                              tGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                                                                                     >  3:48618/1‑148 (MQ=255)
                                              tGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                                                                                     <  2:112665/148‑1 (MQ=255)
                                                              cAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCGCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACt                                                                                    >  8:223171/1‑149 (MQ=255)
                                                              cAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACt                                                                                    >  6:136717/1‑149 (MQ=255)
                                                                   ttCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcctgc                                                                               <  5:304206/149‑1 (MQ=255)
                                                                      gatgggaCTTGCTTCCTGATGCCATGAGATTCCGCTCAGATAATGTAATTGTGTCTGCTTTGCAGCCCTTAGCTGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcctgcct                                                                             <  5:359716/142‑1 (MQ=255)
                                                                         ggCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGCGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGc                                                                         <  1:120378/149‑1 (MQ=255)
                                                                             cTGGCTTCAGGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGTGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTAATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCAt                                                                      >  7:234151/1‑148 (MQ=255)
                                                                                gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGt                                                                  <  5:136717/149‑1 (MQ=255)
                                                                                 cTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtc                                                                 <  1:21871/149‑1 (MQ=255)
                                                                                       gAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCAGCGAAGGGCCGCAAGGCGGCACCGGTTCATTCAATCGTATTGATACCCGACTGCCACTCCTGACTGAGACATAAGTCTAGAAc                                                           >  5:358204/1‑149 (MQ=255)
                                                                                        aaGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTACTGCCTGAGCCATAAGTCTCGAAc                                                           >  1:360750/1‑148 (MQ=255)
                                                                                                 gaTTCCGCTCACATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCaa                                                 <  8:286669/149‑1 (MQ=255)
                                                                                                           aGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCg                                       <  4:335683/149‑1 (MQ=255)
                                                                                                                 tGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                         >  8:105394/1‑125 (MQ=255)
                                                                                                                 tGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                         <  7:105394/125‑1 (MQ=255)
                                                                                                                   tAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCaaa                               >  7:286768/1‑149 (MQ=255)
                                                                                                                         gggCTGCTTTGCAGCGCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCAcc                          <  7:223171/148‑1 (MQ=255)
                                                                                                                          ggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCgg                        <  4:36273/149‑1 (MQ=255)
                                                                                                                                         cTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGc         <  8:286768/149‑1 (MQ=255)
                                                                                                                                          ttCGCGGGCTCGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAATCAGGTACCGGGACCACATCACCCGCCCGTACCCTCCCGcc        >  4:221383/1‑149 (MQ=255)
                                                                                                                                            cgcgGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGccc       >  3:290027/1‑148 (MQ=255)
                                                                                                                                            cgcgGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGccc       <  8:29107/148‑1 (MQ=255)
                                                                                                                                                ggCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAg  <  5:64805/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAG  >  NC_002947/3048457‑3048749

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