Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,571 Δ1 bp 12.2% coding (451/2061 nt) pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,5710C.12.2% 115.9 / 10.3 41coding (451/2061 nt)pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base . (3/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAC  >  NC_002947/5688425‑5688714
                                                                                                                                                  |                                                                                                                                                
gACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCAc                                                                                                                                               <  1:438794/150‑1 (MQ=255)
          aGGATGCGTTCGGCCAGCAACAGCTGGCGCTCGGACGCCGGCTCTTCCGTTGTCGTTGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGc                                                                                                                                    <  1:355323/150‑1 (MQ=255)
           ggATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCg                                                                                                                                   >  2:303180/1‑150 (MQ=255)
             aTGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCAGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCACCAGGTGGCAGCGAc                                                                                                                                 <  1:273637/150‑1 (MQ=255)
                 gTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCa                                                                                                                             >  1:344410/1‑150 (MQ=255)
                        cAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                      >  1:712294/1‑150 (MQ=255)
                         aGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGc                                                                                                                     >  1:725692/1‑150 (MQ=255)
                          gCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCt                                                                                                                    >  2:285940/1‑150 (MQ=255)
                                            ccACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGt                                                                                                  >  1:166808/1‑150 (MQ=255)
                                              aCCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCg                                                                                                <  2:725692/150‑1 (MQ=255)
                                              aCCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCg                                                                                                <  1:409765/150‑1 (MQ=255)
                                                    aCCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTc                                                                                          <  2:344410/150‑1 (MQ=255)
                                                       tGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTTCACCTGCAGCTGCGGCACGGTCTCAGCCAGTGTTGGTGCCACCTGGTACAGCAACAGCACGGTCTGCTCGCTGGCGAGGATGATGTCGGTGGTCtgt                                                                                       >  2:36308/1‑148 (MQ=255)
                                                           tGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCa                                                                                   <  2:166808/150‑1 (MQ=255)
                                                             gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGACAGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                 <  1:623823/150‑1 (MQ=255)
                                                                 gggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGtttt                                                                             >  2:103306/1‑150 (MQ=255)
                                                                  ggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTc                                                                            <  1:303180/150‑1 (MQ=255)
                                                                          tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                    >  1:156556/1‑150 (MQ=255)
                                                                                    gATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGt                                                                                                                         >  2:582584/1‑87 (MQ=255)
                                                                                    gATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGt                                                                                                                         <  1:582584/87‑1 (MQ=255)
                                                                                     aTCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                         <  2:156556/150‑1 (MQ=255)
                                                                                       ctcGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTc                                                       <  2:15945/150‑1 (MQ=255)
                                                                                          gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTgg                                                                                                                                           >  1:713064/1‑64 (MQ=255)
                                                                                          gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTgg                                                                                                                                           <  2:713064/64‑1 (MQ=255)
                                                                                          gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                    <  1:646839/150‑1 (MQ=255)
                                                                                          gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                    >  1:299405/1‑150 (MQ=255)
                                                                                          gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                    >  2:48934/1‑150 (MQ=255)
                                                                                          gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGAGCGGGATCGCCCAGTCGCGGCGTACCGATTCCAt                                                    >  2:645954/1‑150 (MQ=255)
                                                                                                               cTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGTGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTgg                               <  1:645954/150‑1 (MQ=255)
                                                                                                                tGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGgc                              <  2:242089/150‑1 (MQ=255)
                                                                                                                                tctcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAg              <  1:48934/150‑1 (MQ=255)
                                                                                                                                  tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                 <  1:567606/129‑1 (MQ=255)
                                                                                                                                  tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                 >  2:567606/1‑129 (MQ=255)
                                                                                                                                     gCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCa                                                                                   >  2:435955/1‑76 (MQ=255)
                                                                                                                                     gCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTATGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCa                                                                                   <  1:435955/76‑1 (MQ=255)
                                                                                                                                          tgTCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCtt                                                                                        >  2:512345/2‑66 (MQ=25)
                                                                                                                                          tgTCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCtt                                                                                        <  1:512345/65‑1 (MQ=25)
                                                                                                                                           gTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                            <  2:551285/61‑1 (MQ=255)
                                                                                                                                           gTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                            >  1:551285/1‑61 (MQ=255)
                                                                                                                                            tCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCTATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAc  >  1:606372/1‑150 (MQ=255)
                                                                                                                                            tCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAc  >  1:87242/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GACCGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAC  >  NC_002947/5688425‑5688714

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: