Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,351,007 Δ1 bp 55.2% 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.55.2% ‑2.0 / 53.5 29coding (2019/2337 nt)arcBsensory histidine kinase ArcB
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base . (9/7);  total (15/14)
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.07e-01

TTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCA  >  NC_000913/3350873‑3351149
                                                                                                                                      |                                                                                                                                              
gtctGTCCTTCCTCAACAACGCCTGTTTGATCCGGCGCCGTCAGATTCGGCTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTGACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTc                                                                                                                                        <  2:71055/140‑1 (MQ=255)
        ttcctCAACAATGCCTTTTTTATACTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAACTACTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACACGTTCGAGATACAGTCCAAGAATgg                                                                                                                                >  2:354363/1‑143 (MQ=255)
            tCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAata                                                                                                                           >  1:576522/1‑143 (MQ=255)
                    gCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGGTTCGGACTTACAATTTCGAGTCCTGTACCAGCATAGGAAGAGCGAGCaa                                                                                                                   >  2:80189/1‑143 (MQ=255)
                       ttgtttATCCTGCGCCGTCATATTCGCCTCCAGCACGCTGACATCGCCCGGCATCATCTTCTCAAACACCGCTATCCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGCGTACTGTGCCAGCATGGGAATATCGAGCAATGc                                                                                                                <  1:757010/140‑1 (MQ=255)
                            tATCCTGCGCCCTCAGATTCGCCTCCAGCCCGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTg                                                                                                           <  2:2254/143‑1 (MQ=255)
                               ccTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGAt                                                                                                         <  1:429717/143‑1 (MQ=255)
                                  gcgcCGTCAGATTCGACTCCAGCCCGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTAc                                                                                                     <  1:18368/143‑1 (MQ=255)
                                             tgcgccTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCtcttc                                                                                           <  1:271452/138‑1 (MQ=255)
                                              tCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCtcttct                                                                                          <  1:65111/143‑1 (MQ=255)
                                                 aCTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACGCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTGCTCTtgtggc                                                                                       >  1:361613/1‑141 (MQ=255)
                                                   tCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCCATGCTTCTGAGTTACCGATCTCTTCTgtcgtc                                                                                    >  2:534011/1‑143 (MQ=255)
                                                             tGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTActc                                                                           >  2:631463/1‑143 (MQ=255)
                                                             tGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGGGATCAGCTTCGGTCCTACAAGATCGAGATACTGTTCCAGCATGTGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTActc                                                                           >  2:330705/1‑143 (MQ=255)
                                                               aCATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATAGCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTctcc                                                                        <  2:481869/143‑1 (MQ=255)
                                                                 aTAGCCCGGCATCGTCTGCTCCAACACCGCTACCCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTCCTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTctcctc                                                                      <  2:498417/143‑1 (MQ=255)
                                                                  tAGCCCGGCATCATCTTCTCATACCCCGCTAGCCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGCTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTctcctc                                                                      <  2:70295/143‑1 (MQ=255)
                                                                  catcccGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCtca                                                                     <  2:528230/140‑1 (MQ=255)
                                                                          catcatCTTCTCACACACCCCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCt                                                              <  2:615639/143‑1 (MQ=255)
                                                                                       aaCACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATAATCCTGGGTATCCCAGAAt                                                >  1:223423/1‑143 (MQ=255)
                                                                                            cGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCt                                            >  1:678504/1‑143 (MQ=255)
                                                                                            cGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGACTTTCtt                                           >  2:154763/1‑143 (MQ=255)
                                                                                                    cGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATcg                                   >  1:660144/1‑143 (MQ=255)
                                                                                                     gTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTGATGATTTACAGTCCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGAATCCCAGAATTTAATGAACATcgc                                  >  2:54685/1‑143 (MQ=255)
                                                                                                                         cTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGaa               >  2:234055/1‑143 (MQ=255)
                                                                                                                          tACAAGTTCGAGGACTGCTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAg             <  2:421435/143‑1 (MQ=255)
                                                                                                                             aaGTTCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCgg          >  2:809934/1‑143 (MQ=255)
                                                                                                                              aGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCgg          <  1:575312/143‑1 (MQ=255)
                                                                                                                                ttCGAGTACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGCGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGcct       >  2:188172/1‑141 (MQ=255)
                                                                                                                                      aTACTGTTCCCGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCGTCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCa  <  2:816650/143‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                              
TTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTCGAGATACTGTTCCAGCATGGGAATATCGAGCAATGCTTCTGATTTACTGTTCTCTTCTGTCGTCACCGTACTCTCCTCATCATCCTGGGTATCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCA  >  NC_000913/3350873‑3351149

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: