Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,816,917 C→T 20.0% intergenic (+103/+26) PP_5092 → / ← pilT NLP/P60 protein/twitching motility protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,816,9170CT20.0% 68.4 / 11.0 35intergenic (+103/+26)PP_5092/pilTNLP/P60 protein/twitching motility protein
Reads supporting (aligned to +/- strand):  ref base C (8/20);  new base T (3/4);  total (11/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCC  >  NC_002947/5816769‑5817060
                                                                                                                                                    |                                                                                                                                               
gCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc                                                                                                                                                 <  8:329105/149‑1 (MQ=255)
gCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc                                                                                                                                                 <  7:295000/149‑1 (MQ=255)
     tgctgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCGCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCGCCCATACCAGTGGGCGCGGGCGCgcccc                                                                                                                                            <  6:308061‑M2/149‑6 (MQ=255)
            gAAAGAGCAGCAGCTGAGTCAGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg                                                                                                                                     >  5:147032/1‑149 (MQ=255)
                gagCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCt                                                                                                                                 <  7:219263/149‑1 (MQ=255)
                    agcagcTGAGGCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACCGGTAGCGCGGCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGc                                                                                                                             <  5:234231/149‑1 (MQ=255)
                                  gCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGc                                                                                                                <  3:303599/148‑1 (MQ=255)
                                     aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACCCCCCCGAACACACGTATCCCGTGGCTCCACCCCGCAGCGCACAAAACGGTGGGACCGGGCCCGCCCCCGCAACGGGCGGTGCAGGTACGACaccgcgtt                                                                                                            >  4:324673/1‑142 (MQ=255)
                                      tccGTTACGCTGATCGGCGCTAGCCGGTCCCGGCCCCTTCGGGGCCACGCCCGCTCCCACAGGTCGCGCGTCGCTCTCGCGCGCAGGTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTg                                                                                                           <  3:44113/147‑1 (MQ=38)
                                       gCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGa                                                                                                          <  6:147032/149‑1 (MQ=255)
                                        cGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGAt                                                                                                         <  5:308061/149‑1 (MQ=255)
                                                  aTCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGGGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCg                                                                                                <  3:98624/148‑1 (MQ=255)
                                                                   ccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCaa                                                                              <  4:116327/149‑1 (MQ=255)
                                                                    cGGCCCCTGCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGtac                                                                                                                       <  1:283504/107‑1 (MQ=255)
                                                                      gCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCa                                                                                        <  7:361776/136‑1 (MQ=255)
                                                                      gCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGc                                                                           <  7:216542/149‑1 (MQ=255)
                                                                                gggCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTg                                                                                                           >  7:144898/1‑107 (MQ=255)
                                                                                gggCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCAAATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTg                                                                                                           <  8:144898/107‑1 (MQ=255)
                                                                                  gCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaagggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcaccc                                                               <  6:38614‑M2/149‑83 (MQ=255)
                                                                                      gCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACCTACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGCTCCACCCTCGCTCCGTCAGCCTCCGCAACGCCCGCGCCAGCGGCAGc                                                           >  3:133402/1‑149 (MQ=38)
                                                                                       cccGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttg                                                          <  2:358654‑M2/149‑88 (MQ=255)
                                                                                            tCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccac                                                     <  2:298385‑M2/149‑93 (MQ=255)
                                                                                              ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccacca                                                   <  3:336307‑M2/149‑95 (MQ=255)
                                                                                              ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctg                                                                                   >  5:99054‑M2/1‑55 (MQ=255)
                                                                                              ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctg                                                                                   <  6:99054‑M2/117‑63 (MQ=255)
                                                                                              ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcag                                                                    <  7:209433‑M2/132‑78 (MQ=255)
                                                                                              ccACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcag                                                                    >  8:209433‑M2/1‑55 (MQ=255)
                                                                                               cacaGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgctt                                                                                                  >  1:219933‑M2/1‑54 (MQ=255)
                                                                                               cacaGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgctt                                                                                                  <  2:219933‑M2/101‑48 (MQ=255)
                                                                                                  aGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccaccacgta                                               >  4:336307‑M2/1‑51 (MQ=255)
                                                                                                    gTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcacccgcttggccaccacgtagt                                             <  2:6302‑M2/149‑101 (MQ=255)
                                                                                                              gcgctcGCGGGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCTGCCCCGCCCTGc                                   <  6:11608/146‑1 (MQ=255)
                                                                                                              gCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaactcttggccagactcagctgctgctctttcggcagcacccgcttggccaccacgtagtgcttctgcca                                   >  1:358654‑M2/1‑39 (MQ=255)
                                                                                                                       cgcAGCTGACATAGCTGGGGGAGCGGGCGCGGCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGAGCAACCCTCGCTCCCTCTGCCTCCGCAACGCCCCCTCCTTCGTCAGCATCCCCTCCTCCGCCCCGCCTTGCATCACTGaa                          <  6:276210/149‑1 (MQ=255)
                                                                                                                               acaTACCTGTGGTAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACaa                     <  3:324673/146‑1 (MQ=255)
                                                                                                                                           gAGCGGGCGCgcccgcgaaggggccggtaccggatagcgccgatcagcttaacgcttggccagactcagctgctgctctttcggcagcactccgcttggccaccacgtagtgcttctgccagtaaggcttgctcaaggtgtcgatactc      >  1:298385‑M2/1‑10 (MQ=255)
                                                                                                                                               gggCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCAGCGTCAGCATCCCCTCCGCCGCCGCGGCCGGGATCACTGAAGACCACTGCGCCCACCGCCCCTcc  >  4:165812/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCC  >  NC_002947/5816769‑5817060

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