Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,808,460 A→T 10.4% intergenic (+86/+66) 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,4600AT10.4% 131.6 / 5.3 48intergenic (+86/+66)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.81e-01

AGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCC  >  NC_002947/2808340‑2808607
                                                                                                                        |                                                                                                                                                   
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                                                                                          >  1:158741/1‑121 (MQ=255)
                                                           gcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAgatgat                                                                                          <  2:158741/121‑1 (MQ=255)
                                                            cAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAg                                                             <  8:299103/149‑1 (MQ=255)
                                                             aaaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                            >  6:232216/1‑149 (MQ=255)
                                                                   tgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACGAGGGCggg                                                       <  6:48989/148‑1 (MQ=255)
                                                                     aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCTTTGCCCGt                                                                                                                                             >  1:291277/1‑60 (MQ=21)
                                                                     aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCTTTGCCCGt                                                                                                                                             <  2:291277/60‑1 (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                                                                            >  4:105487/1‑104 (MQ=255)
                                                                                          gTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATaaa                                                                            <  3:105487/104‑1 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTcc                            <  6:250586/149‑1 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCACGCAATACCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTcc                            >  6:207378/1‑149 (MQ=255)
                                                                                                 cGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGAc                        <  8:222661/149‑1 (MQ=255)
                                                                                                          ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATc                                       >  5:159121/1‑125 (MQ=255)
                                                                                                          ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATc                                       <  6:159121/125‑1 (MQ=255)
                                                                                                           cgcgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAt                                  >  1:294567/1‑129 (MQ=255)
                                                                                                           cgcgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAt                                  <  2:294567/129‑1 (MQ=255)
                                                                                                             cgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAg            >  7:88714/1‑149 (MQ=255)
                                                                                                              tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                     <  4:168638/120‑1 (MQ=255)
                                                                                                              tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                     >  3:168638/4‑123 (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)
                                                                                                                        |                                                                                                                                                   
AGCAGGGCAGGGTAACCCGCGTCAACTGCGGTTGATGCACGCGGGGGCAACTGCCTTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCC  >  NC_002947/2808340‑2808607

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