Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,351,007 Δ1 bp 35.3% coding (2019/2337 nt) arcB ← sensory histidine kinase ArcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,351,0070A.35.3% 22.9 / 17.3 17coding (2019/2337 nt)arcBsensory histidine kinase ArcB
Reads supporting (aligned to +/- strand):  ref base A (8/3);  new base . (3/3);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

ATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTC  >  NC_000913/3350870‑3351148
                                                                                                                                          |                                                                                                                                             
aTTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATAcgg                                                                                                                                          >  1:229940/1‑141 (MQ=255)
 ttttATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGt                                                                                                                                         >  1:11841/1‑143 (MQ=255)
        tccttcctCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCa                                                                                                                                  >  2:128861/1‑143 (MQ=255)
                            ttttATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCCTCATCTTCTCAAACACAGAATATGCCGTGGTGATCTGCTCCGGTCCTACACGTTCGGGATCCTGTTCCAGCATGGGCATCTCGAGCACTGCt                                                                                                              >  2:140647/1‑143 (MQ=255)
                             ttttTCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCTGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCtt                                                                                                             <  2:118191/143‑1 (MQ=255)
                                                       ccAGCCCGCCGACACGGCCCGGCAGCATCTTCTCAAGCACCGCTGACCCGTCGGCGAGCAGCTTCGGTACTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCa                                                                                  <  2:192245/143‑1 (MQ=255)
                                                         gccACCCTGACATAGCCCTCCACCATCTTCTCGACCACATGGACCTCTTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCACTGCTTCTGATTTACTGTTCTCTTCTGTCGTCAcc                                                                                <  2:174416/141‑1 (MQ=255)
                                                          gcgtgCTGTCATAGCCCGGCTGCTTCTGCTCAACCCCCGCTAGCCCGTCGGTCGTCAGCTTCGGTCCTAGAAGTGCGTGATACTGTTCCAGCATTGGAATATCGAGCAATGCTTCTGTTTTACTGTTCTCTTCTGTCGTCAcc                                                                                <  2:381458/139‑1 (MQ=255)
                                                          gCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTCCTCTTCGGGCGTCACCg                                                                               >  1:53021/1‑143 (MQ=255)
                                                           cACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCg                                                                               >  1:4914/1‑143 (MQ=255)
                                                                              atcatcTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTgg                                                           <  1:61764/143‑1 (MQ=255)
                                                                                 atcTTCTCAAACACCACTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCGTGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGGAATCTCCTCATCATCCTGGGTa                                                        >  1:397972/1‑143 (MQ=255)
                                                                                  tcttctCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTa                                                        >  1:170664/1‑143 (MQ=255)
                                                                                         aaaCACCGCTAACCCGTCGGTGATCAGCTGCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGa                                                 <  1:128861/143‑1 (MQ=255)
                                                                                                                    gCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGGATTTCTTGATCATCGCGGCTCGCGcc                       >  2:42593/1‑143 (MQ=255)
                                                                                                                             cTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGaaa             >  1:234914/1‑143 (MQ=255)
                                                                                                                                cAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGc           >  2:160929/1‑143 (MQ=255)
                                                                                                                                  acttctggaTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGGTTTACTGTCCTCTTCTGTCGTCACCGTACTCTCCTCACCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCgg         <  2:25050/136‑1 (MQ=255)
                                                                                                                                         gaTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCCCCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTc  >  1:257931/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                             
ATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTC  >  NC_000913/3350870‑3351148

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