Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,352,025 G→T 9.8% intergenic (‑230/+34) 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,0250GT 9.8% 156.7 / 6.0 51intergenic (‑230/+34)PP_4707/pnpOsmY‑like protein/polyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base T (2/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCG  >  NC_002947/5351883‑5352173
                                                                                                                                              |                                                                                                                                                    
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                                                                                                                     <  2:687170/150‑1 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                     >  1:647502/1‑150 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                     <  1:1129366/150‑1 (MQ=255)
                          aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                     <  1:233769/150‑1 (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                                                                                          <  2:650472/150‑1 (MQ=255)
                                                     gtggtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAgg                                                                                          <  1:1139717/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                                                                                                        <  2:205169/71‑1 (MQ=255)
                                                                                                                      ccGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAg                                                                                                        >  1:205169/1‑71 (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)
                                                                                                                                              |                                                                                                                                                    
TGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCG  >  NC_002947/5351883‑5352173

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