Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,808,462 A→T 10.4% intergenic (+88/+64) PP_2462 → / ← PP_mr35 hypothetical protein/cspA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,808,4620AT10.4% 132.1 / 5.7 48intergenic (+88/+64)PP_2462/PP_mr35hypothetical protein/cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base T (3/2);  total (27/21)
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.57e-01

AGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTT  >  NC_002947/2808340‑2808610
                                                                                                                          |                                                                                                                                                    
aGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGc                                                                                                                                                    <  1:247171/125‑1 (MQ=255)
aGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGc                                                                                                                                                    >  2:247171/1‑125 (MQ=255)
aGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCgg                                                                                                                            >  7:299103/1‑149 (MQ=255)
 gcagggcaggGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGt                                                                                                                           >  7:222661/1‑149 (MQ=255)
     ggcagggTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCg                                                                                                                       <  3:101707/149‑1 (MQ=255)
        agggTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAg                                                                                                                    <  6:301212/149‑1 (MQ=255)
            tAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCCGGCaa                                                                                                                >  2:163711/1‑149 (MQ=255)
            tAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCaa                                                                                                                >  5:48989/1‑149 (MQ=255)
                ccGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGccc                                                                                                                                                  <  4:237646/111‑1 (MQ=255)
                ccGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGccc                                                                                                                                                  >  3:237646/1‑111 (MQ=255)
                ccGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATccc                                                                                                            <  4:91681/149‑1 (MQ=255)
                                    gCACGCGGGGGCAACTGCCTTGAGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCAt                                                                                         >  2:320638/1‑148 (MQ=255)
                                                aaCTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGc                                                                                                                                                    >  8:50318/1‑77 (MQ=255)
                                                aaCTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGc                                                                                                                                                    <  7:50318/77‑1 (MQ=255)
                                                aaCTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATaaaaa                                                                             >  5:18266/1‑148 (MQ=255)
                                                   tGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAAcc                                                                         >  1:288561/1‑149 (MQ=255)
                                                     ccTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAg                                                                                                                                                     >  6:18358/1‑71 (MQ=255)
                                                     ccTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAg                                                                                                                                                     <  5:18358/71‑1 (MQ=255)
                                                           gcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAgatgat                                                                                             <  2:158741/121‑1 (MQ=255)
                                                           gcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAgatgat                                                                                             >  1:158741/1‑121 (MQ=255)
                                                            cAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAg                                                                <  8:299103/149‑1 (MQ=255)
                                                             aaaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                               >  6:232216/1‑149 (MQ=255)
                                                                   tgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACGAGGGCggg                                                          <  6:48989/148‑1 (MQ=255)
                                                                     aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCTTTGCCCGt                                                                                                                                                <  2:291277/60‑1 (MQ=21)
                                                                     aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCTTTGCCCGt                                                                                                                                                >  1:291277/1‑60 (MQ=21)
                                                                      aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttt                                                       <  2:301750/148‑1 (MQ=255)
                                                                       aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                     <  1:320638/149‑1 (MQ=255)
                                                                             gcgcGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                               >  5:250586/1‑149 (MQ=255)
                                                                              cgcgGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGTTTTTTTatg                                                <  4:205607/147‑3 (MQ=39)
                                                                                         tGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGc                                                                                                          >  3:103150/1‑78 (MQ=255)
                                                                                         tGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGc                                                                                                          <  4:103150/78‑1 (MQ=255)
                                                                                          gTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATaaa                                                                               <  3:105487/104‑1 (MQ=255)
                                                                                          gTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATaaa                                                                               >  4:105487/1‑104 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTcc                               <  6:250586/149‑1 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCACGCAATACCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTcc                               >  6:207378/1‑149 (MQ=255)
                                                                                                 cGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGAc                           <  8:222661/149‑1 (MQ=255)
                                                                                                          ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATc                                          <  6:159121/125‑1 (MQ=255)
                                                                                                          ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATc                                          >  5:159121/1‑125 (MQ=255)
                                                                                                           cgcgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAt                                     >  1:294567/1‑129 (MQ=255)
                                                                                                           cgcgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAt                                     <  2:294567/129‑1 (MQ=255)
                                                                                                             cgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAg               >  7:88714/1‑149 (MQ=255)
                                                                                                              tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                        >  3:168638/4‑123 (MQ=255)
                                                                                                              tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                        <  4:168638/120‑1 (MQ=255)
                                                                                                                gcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCt            >  4:184226/2‑149 (MQ=255)
                                                                                                                 cacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGAAGGCAATCCCGCCCATTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTAAAGCTATCCTAATTACTCCGTACGAACCTGTGCAGCCTg           >  5:3876/1‑149 (MQ=255)
                                                                                                                  acGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTg           >  1:209108/1‑148 (MQ=255)
                                                                                                                    gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCAt        >  4:115301/1‑149 (MQ=255)
                                                                                                                     ggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATg       >  8:16523/1‑149 (MQ=255)
                                                                                                                        aaaGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGcc     >  3:4373/1‑148 (MQ=255)
                                                                                                                          aGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCtt  <  8:15334/149‑1 (MQ=255)
                                                                                                                          |                                                                                                                                                    
AGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTT  >  NC_002947/2808340‑2808610

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