Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,633,500 T→G 12.8% intergenic (+100/‑82) hupB → / → PP_2304 DNA binding regulator subunit beta B/peptidyl‑prolyl cis‑trans isomerase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,633,5000TG12.8% 109.9 / 10.4 47intergenic (+100/‑82)hupB/PP_2304DNA binding regulator subunit beta B/peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base T (16/25);  new base G (3/3);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTA  >  NC_002947/2633360‑2633642
                                                                                                                                            |                                                                                                                                              
gTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCGAGGATgcgc                                                                                                                                       <  1:319378/150‑1 (MQ=255)
      aGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGCCACGAAATGGATTTATTGCCGCAGCGTGAAGTCGGTGCGCAATTTTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCGAGGATGCGCcttttc                                                                                                                                 <  2:739436/150‑3 (MQ=255)
      aGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCctttct                                                                                                                                 <  2:437871/150‑1 (MQ=255)
        gCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCctttcttt                                                                                                                               >  2:965242/1‑150 (MQ=255)
         ccGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCtttt                                                                                                                              <  1:818095/150‑1 (MQ=255)
          cGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTa                                                                                                                             <  2:894400/150‑1 (MQ=255)
                aGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGa                                                                                                                       <  1:530402/150‑1 (MQ=255)
                  gCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGAtt                                                                                                                     <  2:985619/150‑1 (MQ=255)
                     tGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTa                                                                                                                  <  2:817654/150‑1 (MQ=255)
                         aGACGCCGTCAACTAAGTCGTTTTTCAGCCTGACCCGCCCTGGCCTGGTCCGAGCACGGTGATTGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCa                                                                                                              <  1:9735/150‑1 (MQ=255)
                           aCGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACg                                                                                                            >  2:987020/1‑150 (MQ=255)
                               cGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTcc                                                                                                        >  2:87026/1‑150 (MQ=255)
                                     cTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGtt                                                                                                  <  2:551017/150‑1 (MQ=255)
                                          tcgtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGa                                                                                             >  2:134014/1‑150 (MQ=255)
                                           cgtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGAc                                                                                            <  1:987020/150‑1 (MQ=255)
                                            gtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACg                                                                                           >  1:747595/1‑150 (MQ=255)
                                            gtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACg                                                                                           <  1:965242/150‑1 (MQ=255)
                                                   aGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAATGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCAc                                                                                                             <  2:824001/125‑1 (MQ=255)
                                                   aGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAATGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCAc                                                                                                             >  1:824001/1‑125 (MQ=255)
                                                    gCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTa                                                                                   >  1:476524/1‑150 (MQ=255)
                                                     ccGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTAc                                                                                  >  2:1012990/1‑150 (MQ=255)
                                                      cGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACa                                                                                 <  1:134014/150‑1 (MQ=255)
                                                       ggACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAg                                                                                >  2:575935/1‑150 (MQ=255)
                                                                 tggcctggTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTgg                                                                      <  1:1012990/150‑1 (MQ=255)
                                                                 tggcctggTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTgg                                                                      <  1:372195/150‑1 (MQ=255)
                                                                            gAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCAtgc                                                           <  2:289819/150‑1 (MQ=255)
                                                                                         gCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAg                                              <  2:323824/150‑1 (MQ=255)
                                                                                          cGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTTTTGGACGCATGCTGCAAAATATCAgg                                             <  1:338595/150‑1 (MQ=255)
                                                                                                    aCGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTcac                                   <  1:575935/150‑1 (MQ=255)
                                                                                                                 acacGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGc                      <  1:87026/150‑1 (MQ=255)
                                                                                                                          cTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCAtgct                                                          <  2:284996/105‑1 (MQ=255)
                                                                                                                          cTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCAtgct                                                          >  1:284996/1‑105 (MQ=255)
                                                                                                                           tAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGtt                                                                                                  <  2:405170/64‑1 (MQ=255)
                                                                                                                           tAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGtt                                                                                                  >  1:405170/1‑64 (MQ=255)
                                                                                                                            aataaaaGGCGCCTCCTCGGATTCTCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAAAATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatca           <  1:1152009/147‑1 (MQ=255)
                                                                                                                            aaGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatca           >  1:784018/1‑150 (MQ=255)
                                                                                                                               aaaaGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTg                            <  2:149592/130‑1 (MQ=255)
                                                                                                                               aaaaGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTg                            >  1:149592/1‑130 (MQ=255)
                                                                                                                               aaaaGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatcatag        >  1:202590/1‑148 (MQ=255)
                                                                                                                                  aGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGa                                           <  2:394828/112‑1 (MQ=255)
                                                                                                                                  aGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGa                                           >  1:394828/1‑112 (MQ=255)
                                                                                                                                  aGGCGCATCCGAGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGt                                                                            >  2:1049508/1‑79 (MQ=39)
                                                                                                                                  aGGCGCATCCGAGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGt                                                                            <  1:1049508/79‑1 (MQ=39)
                                                                                                                                  aGGCGCATCCGAGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGcc     >  1:490542/1‑150 (MQ=255)
                                                                                                                                  aGGCGCATCCGAGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGcc     >  2:950159/1‑150 (MQ=255)
                                                                                                                                    gcgcATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTACGCGGCAGTCGTCGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCtt   >  2:296284/1‑150 (MQ=255)
                                                                                                                                     cgcATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTa  <  2:747595/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTA  >  NC_002947/2633360‑2633642

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