Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,238,932 T→A 8.0% intergenic (+240/‑41) PP_5701 → / → PP_4616 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,238,9320TA 8.0% 139.4 / 3.4 50intergenic (+240/‑41)PP_5701/PP_4616hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base A (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACTGGCCGAGCGCCTG  >  NC_002947/5238791‑5239080
                                                                                                                                             |                                                                                                                                                      
cATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCg                                                                                                                                                  >  4:160260/1‑148 (MQ=255)
  ttCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCaa                                                                                                                                               <  8:58965/149‑1 (MQ=255)
  ttCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCaa                                                                                                                                               <  8:54425/149‑1 (MQ=255)
             gCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGc                                                                                                                                     >  2:310100/1‑148 (MQ=255)
               cTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCaa                                                                                                                                   >  6:320003/1‑148 (MQ=255)
               cTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCaa                                                                                                                                   >  7:20648/1‑148 (MQ=255)
                 gtgCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCGCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCaaa                                                                                                                                  >  1:239791/1‑147 (MQ=255)
                  tgCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCaaaaa                                                                                                                                <  6:97159/148‑1 (MQ=255)
                    ccTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAgg                                                                                                                              <  1:310100/148‑1 (MQ=255)
                     cTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCa                                                                                                                            <  5:320003/149‑1 (MQ=255)
                         aaCTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATaga                                                                                                                        <  7:5900/149‑1 (MQ=255)
                         aaCTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATa                                                                                                                          >  4:113711/1‑147 (MQ=255)
                         aaCTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATa                                                                                                                          >  2:141139/1‑147 (MQ=255)
                         aaCTGTTTCGGGCAAATGGTTCATGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATa                                                                                                                          >  7:51128/1‑147 (MQ=255)
                          aCTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATagag                                                                                                                       <  8:240721/149‑1 (MQ=255)
                              tttCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAGCCAGCAAAAAGGCATAGAGGtga                                                                                                                   >  2:176982/1‑149 (MQ=255)
                                tCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGtgat                                                                                                                  <  2:243845/148‑1 (MQ=255)
                                               cTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAgt                                                                                                >  2:102712/1‑149 (MQ=255)
                                                tGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAgtg                                                                                               >  5:159738/1‑149 (MQ=255)
                                                tGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGAGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑CAATGGCCGTCGAAgtt                                                                                               >  8:76598/1‑148 (MQ=255)
                                                  gCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAgtg                                                                                               <  8:20648/147‑1 (MQ=255)
                                                    aaaGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCaaa                                                                                                                                  >  2:82540/1‑112 (MQ=255)
                                                          gcTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCg                                                                                       >  5:234734/1‑147 (MQ=255)
                                                              ccGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCgcagc                                                                                   <  6:152311/147‑1 (MQ=255)
                                                               cGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCagcagc                                                                                >  1:247895/1‑149 (MQ=255)
                                                                gCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGcc                                                                               >  7:285252/1‑149 (MQ=255)
                                                                     gcTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGAc                                                                          <  1:100912/149‑1 (MQ=255)
                                                                               tgCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAgcg                                                                  <  6:330798/147‑1 (MQ=255)
                                                                                 ccACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTg                                                              >  8:52134/1‑149 (MQ=255)
                                                                                 ccACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTg                                                              >  7:186455/1‑149 (MQ=255)
                                                                                     gCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCa                                                          <  2:239791/149‑1 (MQ=255)
                                                                                       tACAAGTTCTGCTTGCTCACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATg                                                        <  6:234734/149‑1 (MQ=255)
                                                                                         cAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATg                                                        <  2:247895/147‑1 (MQ=255)
                                                                                          aaTTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGa                                                      >  3:47606/1‑148 (MQ=255)
                                                                                               cTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAgg                                                 <  6:159738/148‑1 (MQ=255)
                                                                                                     gCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACAGCAACACAACCAGCAAAAAGGCAGAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGATCATGGATAAGGCCGAAg                                           >  6:196155/1‑148 (MQ=255)
                                                                                                      cGAACCAACCAGTATTTCCATTACTGTATATATATAAATTAAACCGCAACACAACCCCCAAAAAGGCATATATGTGTTTCCAATGGCCGTCGAAGTGGTGTACTGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGaa                                            <  1:102712/148‑1 (MQ=255)
                                                                                                              ccaGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCa                                 <  8:51128/149‑1 (MQ=255)
                                                                                                               caGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCATCAAAAAGGCATAGAGGTTAT‑‑TAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCAc                                <  7:76598/149‑1 (MQ=255)
                                                                                                               caGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCAc                                <  1:141139/149‑1 (MQ=255)
                                                                                                                       ccATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGAt                        >  8:108450/1‑149 (MQ=255)
                                                                                                                            aCTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTc                    <  6:25911/148‑1 (MQ=255)
                                                                                                                            aCTGTATATAGATACAGAAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCt                     <  5:66159/147‑1 (MQ=255)
                                                                                                                            aCTGTATATAGATACAGAAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCt                     >  6:71086/1‑147 (MQ=255)
                                                                                                                            aCTGTATATAGATACAGAAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCt                     >  6:66159/1‑147 (MQ=255)
                                                                                                                            aCTGTATATAGATACAGAAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCt                     <  5:71086/147‑1 (MQ=255)
                                                                                                                                  tataGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACTgg             <  7:52134/149‑1 (MQ=255)
                                                                                                                                   ataGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACTGGc            >  3:84870/1‑149 (MQ=255)
                                                                                                                                      gATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACTGGCCGa         <  8:285252/149‑1 (MQ=255)
                                                                                                                                      gATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGGGAT‑‑GAATCGCCGTCGAACTCGTGTACAGCAGCAGCCGCGACCCGGCGCTCTTGATCATGGAGAAGGCCGCAGCAGACCGTCACGACCAGATCCTCGCACAGGcc           >  5:94211/1‑147 (MQ=255)
                                                                                                                                             tAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACTGGCCGAGCGCCTg  >  1:55619/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
CATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTAAACCGCAACACAACCAGCAAAAAGGCATAGAGGTGAT‑‑GAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACTGGCCGAGCGCCTG  >  NC_002947/5238791‑5239080

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