Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,801:1 +T 9.3% intergenic (+67/+60) PP_2445 → / ← PP_2446 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,793,8011.T 9.3% 162.6 / 11.3 43intergenic (+67/+60)PP_2445/PP_2446membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (16/23);  new base T (2/2);  total (18/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.88e-01

GTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGC  >  NC_002947/2793659‑2793945
                                                                                                                                               |                                                                                                                                                
gTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGCTGTTAACAGGGGTGCCACCCGCCCCTCGTCCTGGCAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAAcgc                                                                                                                                           <  1:647175/150‑1 (MQ=255)
gTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAAcgc                                                                                                                                           <  1:2425816/150‑1 (MQ=255)
  tctcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCa                                                                                                                                               <  2:2416847/144‑1 (MQ=255)
  tctcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCa                                                                                                                                               >  1:2416847/1‑144 (MQ=255)
  tctcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGa                                                                                                                                         <  2:614817/150‑1 (MQ=255)
   ctcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                                                                        <  1:1072148/150‑1 (MQ=255)
   ctcGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                                                                        <  2:31772/150‑1 (MQ=255)
      gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATAc                                                                                                                                     >  2:1435623/1‑150 (MQ=255)
      gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATAc                                                                                                                                     <  2:2382452/150‑1 (MQ=255)
      gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATAc                                                                                                                                     <  1:2307887/150‑1 (MQ=255)
      gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGAACAGGTATTCGCGTGGCACTGGCAACGCGAATAc                                                                                                                                     <  1:1677189/150‑1 (MQ=255)
        cGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCt                                                                                                                                   <  2:208290/150‑1 (MQ=255)
        cGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCt                                                                                                                                   <  2:2384973/150‑1 (MQ=255)
        cGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCt                                                                                                                                   <  1:2010946/150‑1 (MQ=255)
            tATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAAt                                                                                                                                       >  2:2229238/1‑142 (MQ=38)
                  ggTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCggg                                                                                                                         <  1:1435623/150‑1 (MQ=255)
                         aaCATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAAt                                                                                                                                       >  1:375809/1‑129 (MQ=255)
                         aaCATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAAt                                                                                                                                       <  2:375809/129‑1 (MQ=255)
                                    ccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCtg                                                                                                                                  >  1:1297088/1‑123 (MQ=25)
                                    ccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCtg                                                                                                                                  >  1:1285057/1‑123 (MQ=25)
                                    ccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCtg                                                                                                                                  <  2:1285057/123‑1 (MQ=25)
                                    ccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCtg                                                                                                                                  <  2:1297088/123‑1 (MQ=25)
                                                                              tAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTgg                                                                                                                                                 >  1:1711305/1‑66 (MQ=255)
                                                                              tAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTgg                                                                                                                                                 <  2:1711305/66‑1 (MQ=255)
                                                                               aaaaaGGGGCGACTCGCGACCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCgggg                                                            >  2:585922/1‑150 (MQ=255)
                                                                                 aaaGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTACAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGcc                                                          >  2:932170/1‑150 (MQ=255)
                                                                                    ggggCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCg                                                                                                                           >  2:1560029/1‑82 (MQ=255)
                                                                                    ggggCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCg                                                                                                                           <  1:1560029/82‑1 (MQ=255)
                                                                                                       ttGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTACAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGc                                    <  1:932170/150‑1 (MQ=255)
                                                                                                                   ccGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCgg                        >  2:2052317/1‑150 (MQ=255)
                                                                                                                     gCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTg                      <  1:2052317/150‑1 (MQ=255)
                                                                                                                       tGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAg                                                                                   >  2:1557713/1‑87 (MQ=255)
                                                                                                                       tGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAg                                                                                   <  1:1557713/87‑1 (MQ=255)
                                                                                                                           caGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggt                >  2:699147/1‑150 (MQ=255)
                                                                                                                              gTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggtggt             <  1:699147/150‑1 (MQ=255)
                                                                                                                                aTTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGgtgt           >  2:2015720/1‑150 (MQ=255)
                                                                                                                                 ttCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGtgtg          >  2:1373785/1‑150 (MQ=255)
                                                                                                                                      gTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGa     >  1:1937924/1‑150 (MQ=255)
                                                                                                                                       tGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTgg                                                <  2:1501421/106‑1 (MQ=255)
                                                                                                                                       tGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTgg                                                >  1:1501421/1‑106 (MQ=255)
                                                                                                                                         gaaCTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAcgc  <  1:585922/148‑1 (MQ=255)
                                                                                                                                         gCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCggg                                                             <  2:1272538/91‑1 (MQ=255)
                                                                                                                                         gCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCggg                                                             >  1:1272538/1‑91 (MQ=255)
                                                                                                                                           aCTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGcc                                                          >  2:1710575/1‑92 (MQ=255)
                                                                                                                                           aCTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGcc                                                          <  1:1710575/92‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGC  >  NC_002947/2793659‑2793945

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