Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,352,028 A→C 9.6% intergenic (‑233/+31) PP_4707 ← / ← pnp OsmY‑like protein/polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,352,0280AC 9.6% 129.8 / 8.1 52intergenic (‑233/+31)PP_4707/pnpOsmY‑like protein/polyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base C (2/3);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

TGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGGTGA  >  NC_002947/5351883‑5352177
                                                                                                                                                 |                                                                                                                                                     
tgctTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAggg                                                                                                                                                   >  2:488676/1‑150 (MQ=255)
tgctTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAggg                                                                                                                                                   <  1:91223/150‑1 (MQ=255)
     cTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCtt                                                                                                                                              <  2:280198/150‑1 (MQ=255)
      tGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCttt                                                                                                                                             >  1:1090147/1‑150 (MQ=255)
      tGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCttt                                                                                                                                             >  2:233769/1‑150 (MQ=255)
        ttATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTtct                                                                                                                                           <  2:1090147/150‑1 (MQ=255)
           tttCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTa                                                                                                                                        >  2:456317/1‑150 (MQ=255)
               cTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGc                                                                                                                                    <  1:456317/150‑1 (MQ=255)
                ttttAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCg                                                                                                                                   >  2:1129366/1‑150 (MQ=255)
                 tttAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGt                                                                                                                                  <  1:488676/150‑1 (MQ=255)
                      aTCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCGGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGt                                                                                                                             >  1:694413/1‑150 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                         >  2:1119016/1‑150 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                         <  1:233769/150‑1 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                         <  2:687170/150‑1 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                         <  1:1129366/150‑1 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                         >  1:647502/1‑150 (MQ=255)
                              gTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAg                                                                                                                     >  1:179648/1‑150 (MQ=255)
                              gTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAg                                                                                                                     >  2:856034/1‑150 (MQ=255)
                                      atatCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGa                                                                                                             <  2:179648/150‑1 (MQ=255)
                                            atcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCtt                                                                                                       >  1:534100/1‑150 (MQ=255)
                                             tcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTg                                                                                                      >  1:428731/1‑150 (MQ=255)
                                              cACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                     <  2:647502/150‑1 (MQ=255)
                                                    ggtggtCGTGCATTTTGCAATCAGCCCGCACGTACTGATTTGAAGATGCAAGTCGGGAAAAGGTTCCCGTCCCCATGAAAAAAGGGCTCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAg                                                                                               <  2:96393/150‑1 (MQ=255)
                                                     gtggtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAgg                                                                                              <  1:1139717/150‑1 (MQ=255)
                                                     gtggtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAgg                                                                                              <  2:650472/150‑1 (MQ=255)
                                                           gTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTc                                                                                        <  2:534100/150‑1 (MQ=255)
                                                           gTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCCTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTTGCCGCAGGGACGTc                                                                                        >  2:523478/1‑150 (MQ=255)
                                                               aTTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCC‑‑ATCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGAt                                                                                  <  2:190100/150‑1 (MQ=255)
                                                                                gTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACg                                                                   >  2:1145736/1‑150 (MQ=255)
                                                                                    aCCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCcg                                                               >  2:939302/1‑150 (MQ=255)
                                                                                      cTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCAGAGGCCTTGGCCGCAGGGGCGTCGTTTATCGACATCGTGAGACGGCcgcg                                                             >  2:1164090/1‑150 (MQ=255)
                                                                                              aagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCACAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCa                                                     >  2:359060/1‑150 (MQ=255)
                                                                                                  tgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTc                                                 <  1:160942/150‑1 (MQ=255)
                                                                                                            cAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGTCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGc                                       >  2:781641/1‑150 (MQ=255)
                                                                                                               aaGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACt                                    >  2:393256/1‑150 (MQ=255)
                                                                                                                      ccGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAg                                                                                                            >  1:205169/1‑71 (MQ=255)
                                                                                                                      ccGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAg                                                                                                            <  2:205169/71‑1 (MQ=255)
                                                                                                                       cGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTAtt                                                                                                                                      >  1:421609/1‑44 (MQ=255)
                                                                                                                       cGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTAtt                                                                                                                                      <  2:421609/44‑1 (MQ=255)
                                                                                                                        gTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCACAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTcc                           <  1:359060/150‑1 (MQ=255)
                                                                                                                                   aaaaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGccttcctt                <  1:393256/150‑1 (MQ=255)
                                                                                                                                    aaaGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTg               >  1:985444/1‑150 (MQ=255)
                                                                                                                                    aaaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTg               >  2:1154624/1‑150 (MQ=255)
                                                                                                                                    aaaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTg               >  2:1032112/1‑150 (MQ=255)
                                                                                                                                     aaGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGa              <  2:985444/150‑1 (MQ=255)
                                                                                                                                      aGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCtt                                                                                                       >  2:236882/1‑60 (MQ=38)
                                                                                                                                      aGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCtt                                                                                                       <  1:236882/60‑1 (MQ=38)
                                                                                                                                       gggCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAgg            <  1:1119016/150‑1 (MQ=255)
                                                                                                                                         gCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGAt          >  1:531101/1‑150 (MQ=255)
                                                                                                                                          cccTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATg         >  1:1024945/1‑150 (MQ=255)
                                                                                                                                             tGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCg      <  1:856034/150‑1 (MQ=255)
                                                                                                                                               cGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGGt    <  1:523478/150‑1 (MQ=255)
                                                                                                                                                 aaGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCTCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGGTGa  >  1:336103/1‑150 (MQ=255)
                                                                                                                                                 aaGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGGTGa  <  1:180923/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGGTGA  >  NC_002947/5351883‑5352177

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