Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,286,475 T→C 9.1% intergenic (‑11/+34) PP_3757 ← / ← PP_3758 chemotaxis protein CheY/two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,286,4750TC 9.1% 118.9 / 2.6 44intergenic (‑11/+34)PP_3757/PP_3758chemotaxis protein CheY/two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (18/22);  new base C (2/2);  total (20/24)
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.49e-01

CAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCA  >  NC_002947/4286327‑4286620
                                                                                                                                                    |                                                                                                                                                 
cAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCtt                                                                                                                                                   <  3:298879/149‑1 (MQ=255)
    aTGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGaaa                                                                                                                                               <  7:252474/149‑1 (MQ=255)
       gCAAACTGCTTCTAGCCATCCTCAGACACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATTCTAATCGCAGGTTCATCCTCGACGACAAGAACCACATATTGTGCATCTTCACTCATGTTAACGCCTTGAAATa                                                                                                                                             <  1:350182/148‑1 (MQ=255)
        cAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACt                                                                                                                                           >  6:133015/1‑149 (MQ=255)
        cAAACGCCTCCACGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCACGTAATCTCCCAAAATCATGCGAATCCCCGGTTCATCCTCGACGACCAGAACAACAACTTGTGCATCTTCACTCAAGTAAACGCCTTGAAAAAc                                                                                                                                            >  2:350182/1‑148 (MQ=255)
               cTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGtt                                                                                                                                    <  2:150722/149‑1 (MQ=255)
                 gCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACt                                                                                                                                           >  4:245035/1‑140 (MQ=255)
                 gCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACt                                                                                                                                           <  3:245035/140‑1 (MQ=255)
                     gCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCAt                                                                                                                              <  6:313605/149‑1 (MQ=255)
                     gCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCAt                                                                                                                              <  8:215350/149‑1 (MQ=255)
                       cATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATACTCGACGACCAGAACAACATCTTGAGCATCATCACTCATGTTAACGCCTTGAAATACATACAGATGGCCAtaa                                                                                                                            >  2:294995/1‑147 (MQ=255)
                          ccTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGAc                                                                                                                         >  7:196685/1‑149 (MQ=255)
                              aGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAAGCCTCCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGTAATACTTACAGTTGGCCATCAGAccc                                                                                                                       <  1:294995/147‑1 (MQ=255)
                               gCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAg                                                                                                                    <  8:317865/149‑1 (MQ=255)
                                ccaccaGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAgg                                                                                                                   <  3:79144/149‑1 (MQ=255)
                                              tACCCCTCACCCGCCAGGGAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGtt                                                                                                     >  8:277024/1‑149 (MQ=255)
                                                 cccTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGAcc                                                                                                                        >  2:241246/1‑127 (MQ=255)
                                                 cccTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGAcc                                                                                                                        <  1:241246/127‑1 (MQ=255)
                                                   cTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTgg                                                                                                                                  <  7:219482/115‑1 (MQ=255)
                                                   cTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCcg                                                                                                >  6:165056/1‑149 (MQ=255)
                                                   cTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATATTCACTCATGTTAACGCCTTGAAATACTTACAGTTgg                                                                                                                                  >  8:219482/1‑115 (MQ=255)
                                                      aCCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCcg                                                                                                >  1:126585/1‑146 (MQ=255)
                                                           ccAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGt                                                                                        >  2:340077/1‑149 (MQ=255)
                                                             aGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCg                                                                                      >  7:15361/1‑149 (MQ=255)
                                                                  aTCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGtttt                                                                                  <  5:133015/148‑1 (MQ=255)
                                                                     gcgcAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTc                                                                              <  3:15136/149‑1 (MQ=255)
                                                                      cgcAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCgcgc                                                                                             >  1:187384/1‑133 (MQ=255)
                                                                      cgcAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCgcgc                                                                                             <  2:187384/133‑1 (MQ=255)
                                                                       gcAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCtt                                                                                                                                                   <  7:78120/78‑1 (MQ=255)
                                                                       gcAAAATCATGCGAATCGCCGGTTCATCATCGACGACCAGAACCACATCTTGTGAATCATAACTCATGTTAACGCCtt                                                                                                                                                   >  8:78120/1‑78 (MQ=255)
                                                                          aaaTCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGa                                                                         <  8:196685/149‑1 (MQ=255)
                                                                                   gAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGc                                                                  >  8:79210/1‑147 (MQ=255)
                                                                                                   cTCGACGACGAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGt                                                >  8:266419/1‑149 (MQ=255)
                                                                                                       acgacCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAAt                                            <  7:277024/149‑1 (MQ=255)
                                                                                                       acgacCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAAt                                            <  8:15361/149‑1 (MQ=255)
                                                                                                            cAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTgg                                                                                         >  6:320534/1‑99 (MQ=255)
                                                                                                            cAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTgg                                                                                         <  5:320534/99‑1 (MQ=255)
                                                                                                                   acaTCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTg                                <  5:165056/149‑1 (MQ=255)
                                                                                                                                   cACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCtt                    <  2:126585/145‑1 (MQ=255)
                                                                                                                                    aCTCATGTTAACGCCTCGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAAc               >  2:339686/1‑149 (MQ=255)
                                                                                                                                    aCTCATGTTAACGCCTCGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAAc               >  2:334478/1‑149 (MQ=255)
                                                                                                                                        aTGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTg            >  4:327114/1‑148 (MQ=255)
                                                                                                                                             aaCGCCTCGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATcc       <  1:334478/148‑1 (MQ=255)
                                                                                                                                             aaCGCCTCGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATcc       <  1:339686/148‑1 (MQ=255)
                                                                                                                                                 ccTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGTATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCa  <  1:340077/149‑1 (MQ=255)
                                                                                                                                                  cTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCtt                    >  6:38411/1‑130 (MQ=255)
                                                                                                                                                  cTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCtt                    <  5:38411/130‑1 (MQ=255)
                                                                                                                                                    tGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAAtt                                                       <  2:249506/93‑1 (MQ=255)
                                                                                                                                                    tGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAAtt                                                       >  1:249506/1‑93 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTACACAGAGGTTGCCCGCGCCTGGTCGTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCA  >  NC_002947/4286327‑4286620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑