Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 900,536 T→G 11.8% intergenic (+28/+17) PP_0784 → / ← cysZ hypothetical protein/sulfate transport protein CysZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947900,5360TG11.8% 137.6 / 10.7 51intergenic (+28/+17)PP_0784/cysZhypothetical protein/sulfate transport protein CysZ
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base G (3/3);  total (28/23)
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.54e-01

CTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCTT  >  NC_002947/900391‑900681
                                                                                                                                                  |                                                                                                                                                 
cTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTg                                                                                                                                                  >  2:229683/1‑147 (MQ=255)
   ggCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAg                                                                                                                                               <  1:171432/147‑1 (MQ=255)
    gCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCAGGGCAGcc                                                                                                                                             <  4:164883/148‑1 (MQ=255)
      gTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCt                                                                                                                                          <  6:4536/149‑1 (MQ=255)
          ggcggGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCATT‑CCGCAAAAAGGGGCTGCAGGGCAGCCCCttttt                                                                                                                                      <  1:207344/149‑5 (MQ=255)
               gTCAAGCTACCTTTCTAACCTCGAGAAGTTTATGGAAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACACTACAGTTTATCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCAGGGCAGCCCCtttttgcg                                                                                                                                   <  4:314324/147‑8 (MQ=255)
               gTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTc                                                                                                                                 <  8:303176/149‑1 (MQ=255)
               gTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTc                                                                                                                                 <  1:68329/149‑1 (MQ=255)
                         cTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAgctccgctc                                                                                                                       >  4:266145/1‑149 (MQ=255)
                         cTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAgctccgctc                                                                                                                       >  8:116884/1‑149 (MQ=255)
                           tGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCAACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTccc                                                                                                                     >  1:316624/1‑149 (MQ=255)
                            gCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGGAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAgctccgctcc                                                                                                                      <  1:229683/147‑1 (MQ=255)
                                 ccTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCAccc                                                                                                               <  7:116884/149‑1 (MQ=255)
                                     gagaAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGa                                                                                                            <  4:281949/148‑1 (MQ=255)
                                           gTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACg                                                                                                      >  7:233793/1‑148 (MQ=255)
                                                ggCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATATAATTATAAGGGATAGCCATGCAGCCCATCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGcc                                                                                                 >  1:105946/1‑149 (MQ=255)
                                                    gTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCTCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGc                                                                                            <  6:82287/149‑1 (MQ=255)
                                                       gAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCAc                                                                                         >  1:151373/1‑149 (MQ=255)
                                                       gAACATCTAATACTGCTTTCTCAGATTGCCAAACTCGAACACTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCCCTCCTTCAGCTCCGCACCCGCACCCAGAACAACCTAGACCCGGCCa                                                                                          >  2:123617/1‑148 (MQ=255)
                                                             cTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTgccgc                                                                                   >  1:91276/1‑149 (MQ=255)
                                                                aTACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCgg                                                                                >  4:260445/1‑149 (MQ=255)
                                                                    tGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcat                                                                             >  2:189017/1‑148 (MQ=255)
                                                                      cTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGaa                                                                                                           >  2:34149/1‑116 (MQ=255)
                                                                      cTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGaa                                                                                                           <  1:34149/116‑1 (MQ=255)
                                                                      cTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatcat                                                                          >  6:208204/1‑149 (MQ=255)
                                                                      cTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatc                                                                            <  1:139381/147‑1 (MQ=255)
                                                                           tcAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCa                                                                     >  7:112885/1‑149 (MQ=255)
                                                                            cAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCAc                                                                    <  1:269969/149‑1 (MQ=255)
                                                                               aTTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGtt                                                                 <  2:105946/149‑1 (MQ=255)
                                                                               aTTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGtt                                                                 >  7:294105/1‑149 (MQ=255)
                                                                                   ccACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACg                                                             >  8:188780/1‑149 (MQ=255)
                                                                                   ccACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAACCAAT‑CCCCAAAAAGGGGCTGCCGTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCAc                                                              <  2:316624/148‑1 (MQ=255)
                                                                                     aCCCTCGAAACCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGATCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGcc                                                           >  5:314981/1‑149 (MQ=255)
                                                                                            aaCCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAat                                                     >  3:811/1‑148 (MQ=255)
                                                                                                cTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTCCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAAtcatca                                                <  2:155681/149‑1 (MQ=255)
                                                                                                      tttCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAgcgc                                            >  4:76811/1‑147 (MQ=255)
                                                                                                            cGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCaggtag                                     <  6:140599/148‑1 (MQ=255)
                                                                                                              aCTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCaggtaggt                                   >  4:98855/1‑148 (MQ=255)
                                                                                                              aCTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTg                                  <  2:91276/149‑1 (MQ=255)
                                                                                                                  cTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgc                              >  5:203052/1‑149 (MQ=255)
                                                                                                                   tAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAAGCATCAGCGCCAGGTAGGTGATgcc                             >  2:285774/1‑149 (MQ=255)
                                                                                                                    aaaatcaaaT‑CCGAAAAAAGGGGCTGCCCTGCCGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgccg                            <  2:151373/142‑1 (MQ=255)
                                                                                                                         gCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCaaa                       <  7:188780/149‑1 (MQ=255)
                                                                                                                            aT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGcc                    >  3:220102/1‑149 (MQ=255)
                                                                                                                               ccGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCa                  <  1:189017/149‑1 (MQ=255)
                                                                                                                                      gaGGGGCTGCAGGGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTg           >  3:247005/2‑149 (MQ=255)
                                                                                                                                      aaGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTg           <  8:112885/149‑1 (MQ=255)
                                                                                                                                       aGGGGCTGCAGGGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGc          >  3:306817/1‑149 (MQ=255)
                                                                                                                                       aGGGGCTGCAGGGCAGCCCCTCGCTCCTTCAGCTCCCCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGc          >  6:216423/1‑149 (MQ=255)
                                                                                                                                             tGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAgcgc    >  6:105880/1‑149 (MQ=255)
                                                                                                                                               cccTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCta  >  4:66089/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAAT‑CCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCTT  >  NC_002947/900391‑900681

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