Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,633,500 T→G 17.4% 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,5000TG17.4% 40.3 / 5.9 23intergenic (+100/‑82)hupB/PP_2304DNA binding regulator subunit beta B/peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base T (9/10);  new base G (2/2);  total (11/12)
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.48e-01

AGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATC  >  NC_002947/2633358‑2633635
                                                                                                                                              |                                                                                                                                       
aGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATgc                                                                                                                                  <  2:678915/150‑1 (MQ=255)
                       tGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTa                                                                                                           <  1:379805/150‑1 (MQ=255)
                                     aaCTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCgg                                                                                             >  1:582062/1‑150 (MQ=255)
                                            tcgtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGa                                                                                      >  2:286776/1‑150 (MQ=255)
                                              gtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATc                                                                                                                    <  1:292033/118‑1 (MQ=255)
                                              gtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATc                                                                                                                    >  2:292033/1‑118 (MQ=255)
                                                    cAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTgg                                                                              <  2:582062/150‑1 (MQ=255)
                                                      ctcGGACCAGGCCAGGCCGGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCa                                                                                                                   >  1:217267/3‑111 (MQ=255)
                                                      ctcGGACCAGGCCAGGCCGGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCa                                                                                                                   <  2:217267/109‑1 (MQ=255)
                                                               ggcctggcctGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGtgt                                                                   >  2:610924/1‑148 (MQ=255)
                                                                       ctggTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCGAGGATgcgc                                                                                                                                >  1:467436/1‑81 (MQ=255)
                                                                       ctggTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCGAGGATgcgc                                                                                                                                <  2:467436/81‑1 (MQ=255)
                                                                                                  gcAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCa                                        <  2:118596/142‑1 (MQ=255)
                                                                                                  gcAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCa                                        >  1:118596/1‑142 (MQ=255)
                                                                                                    aaaCGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGGCTCAGTGGTACAGCCGTTCCTGGGGGACGCATGCTGCAAAATATCAGGGACAAttt                              >  1:384544/1‑149 (MQ=255)
                                                                                                     aaCGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACTCTCCGCGGTTGTCGCTGCACTGGTACCGCCGTCTTTGGGCGACGCACTCTGCTAAATATCAGGGACAATTca                             <  2:650069/150‑1 (MQ=255)
                                                                                                         ttCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAg                         >  2:280827/1‑150 (MQ=255)
                                                                                                                    caggtgtGCTAAGAAAAGGCGCATCCTCGGATGGGCCTTTCTTTTATCAGGATTCTACCCGCGCTCCGCGGTTGTCCACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGcc              <  2:299891/144‑1 (MQ=255)
                                                                                                                                gAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTc                                                                                        <  2:142525/64‑1 (MQ=255)
                                                                                                                                gAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTc                                                                                        >  1:142525/1‑64 (MQ=255)
                                                                                                                                gAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCatcatc  <  1:280827/150‑1 (MQ=255)
                                                                                                                                      gcgcATCCGAGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTgg                                                               >  2:262157/1‑83 (MQ=255)
                                                                                                                                      gcgcATCCGAGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTgg                                                               <  1:262157/83‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                       
AGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATC  >  NC_002947/2633358‑2633635

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