Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,576 T→A 11.8% Q149L (CAG→CTG)  pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,5760TA11.8% 115.7 / 5.5 41Q149L (CAG→CTG) pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base A (3/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCA  >  NC_002947/5688435‑5688723
                                                                                                                                             |                                                                                                                                                    
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                                                                                                                                  <  1:582584/87‑1 (MQ=255)
                                                                          gATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGt                                                                                                                                  >  2:582584/1‑87 (MQ=255)
                                                                           aTCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                                  <  2:156556/150‑1 (MQ=255)
                                                                             ctcGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTc                                                                <  2:15945/150‑1 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTgg                                                                                                                                                    <  2:713064/64‑1 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTgg                                                                                                                                                    >  1:713064/1‑64 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                             <  1:646839/150‑1 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                             >  2:48934/1‑150 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                             >  1:299405/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                                          >  2:567606/1‑129 (MQ=255)
                                                                                                                        tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                          <  1:567606/129‑1 (MQ=255)
                                                                                                                           gCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCa                                                                                            >  2:435955/1‑76 (MQ=255)
                                                                                                                           gCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTATGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCa                                                                                            <  1:435955/76‑1 (MQ=255)
                                                                                                                                tgTCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCtt                                                                                                 <  1:512345/65‑1 (MQ=25)
                                                                                                                                tgTCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCtt                                                                                                 >  2:512345/2‑66 (MQ=25)
                                                                                                                                 gTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                     >  1:551285/1‑61 (MQ=255)
                                                                                                                                 gTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                     <  2:551285/61‑1 (MQ=255)
                                                                                                                                  tCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCTATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAc           >  1:606372/1‑150 (MQ=255)
                                                                                                                                  tCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAc           >  1:87242/1‑150 (MQ=255)
                                                                                                                                           ccTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCa  >  1:121323/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCA  >  NC_002947/5688435‑5688723

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: