Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,190 T→C 7.4% Q334Q (CAA→CAG)  flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,1900TC 7.4% 159.3 / 3.0 54Q334Q (CAA→CAG) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base T (30/20);  new base C (2/2);  total (32/22)
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.27e-01

TTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTTTAGCGACCACTTTCGG  >  NC_000913/1965043‑1965336
                                                                                                                                                   |                                                                                                                                                  
ttACTCTTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGtt                                                                                                                                                   <  5:324693/149‑1 (MQ=255)
    tCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTCCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGAcc                                                                                                                                               >  8:408019/1‑149 (MQ=255)
      aTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGtt                                                                                                                                                   >  8:281932/1‑143 (MQ=255)
      aTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGtt                                                                                                                                                   <  7:281932/143‑1 (MQ=255)
      aTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCg                                                                                                                                             <  6:274242/149‑1 (MQ=255)
      aTGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTCCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGAccc                                                                                                                                              <  7:408019/148‑1 (MQ=255)
       tGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCgg                                                                                                                                            >  5:198468/1‑149 (MQ=255)
             ggTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGAttt                                                                                                                                       <  6:198468/148‑1 (MQ=255)
             ggTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttg                                                                                                                                      >  2:40120/1‑149 (MQ=255)
               tttCTCGTTAATAAAATCCAGGGCTTCCAGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgtt                                                                                                                                    >  6:28285/1‑149 (MQ=255)
                ttCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCCGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgtt                                                                                                                                    >  3:88365/1‑148 (MQ=255)
                     gTTAATAAAATCCAGGGCTTCCGGCACCGGAATATGAGTAGGTTGTTCAGGGCGCTGTCCACCAGCCAGACGCGAGCGTTTCAGTTGCCAGACCCAGTCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGTTTTGTTTAc                                                                                                                                 <  4:101613/146‑1 (MQ=255)
                     gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGCCCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCaa                                                                                                                              >  4:227913/1‑149 (MQ=255)
                     gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGAc                                                                                                                                 <  2:287601/146‑1 (MQ=255)
                     gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGAc                                                                                                                                 <  8:339845/146‑1 (MQ=255)
                     gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCaa                                                                                                                              >  4:228584/1‑149 (MQ=255)
                          taaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACCGGGCGCTGTCCACCCGCCAGAAGCCAGCGTTTCAGTTGCCAGCCCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGCTGACCCGGGATATGTTGACCAAtctac                                                                                                                         >  5:350906/1‑147 (MQ=255)
                             aaTCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgt                                                                                                                                     >  6:178990/1‑134 (MQ=255)
                             aaTCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATttgt                                                                                                                                     <  5:178990/134‑1 (MQ=255)
                                         ccGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGCCCAATCTCCGCATGTCGATACAgc                                                                                                          >  3:72183/1‑149 (MQ=255)
                                          cGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGCTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAgc                                                                                                          <  6:350906/148‑1 (MQ=255)
                                            gCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCt                                                                                                       >  3:272854/1‑149 (MQ=255)
                                               ccGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTcgc                                                                                                    <  6:1583/149‑1 (MQ=255)
                                                   aaGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCa                                                                                                <  3:227913/149‑1 (MQ=255)
                                                          gTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTAGAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggcgg                                                                                         >  5:291596/1‑149 (MQ=255)
                                                           tAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGCTGACCCGGGATTTGTTGACCAAtctc                                                                                                                          >  4:362918/1‑115 (MQ=255)
                                                           tAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGCTGACCCGGGATTTGTTGACCAAtctc                                                                                                                          <  3:362918/115‑1 (MQ=255)
                                                              gTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggcggc                                                                                        <  2:51886/146‑1 (MQ=255)
                                                              gTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggcggc                                                                                        >  1:51886/1‑146 (MQ=255)
                                                              gTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCt                                                                                     >  4:159537/1‑149 (MQ=255)
                                                                  tACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCaa                                                                                 <  6:103927/149‑1 (MQ=255)
                                                                          gcTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCgg                                                                          >  1:11246/1‑148 (MQ=255)
                                                                            tGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGAc                                                                       >  3:265759/1‑149 (MQ=255)
                                                                            tGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGAc                                                                       >  7:30375/1‑149 (MQ=255)
                                                                                  ccagccagACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtt                                                                 <  1:40120/149‑1 (MQ=255)
                                                                                   cagccagACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAttt                                                                >  4:101965/1‑149 (MQ=255)
                                                                                    agccagACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGCCGTTAtttt                                                               >  8:166646/1‑149 (MQ=255)
                                                                                         gACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAgcg                                                          >  2:240515/1‑149 (MQ=255)
                                                                                              cAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAAt                                                     >  7:189742/1‑149 (MQ=255)
                                                                                                 cGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTc                                                  <  3:101965/149‑1 (MQ=255)
                                                                                                      cAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGa                                              <  3:228584/148‑1 (MQ=255)
                                                                                                       aGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtt                                                                 >  5:43161/1‑128 (MQ=255)
                                                                                                       aGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTAtt                                                                 <  6:43161/128‑1 (MQ=255)
                                                                                                                aCCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCCGCGCCAATTTCACGAATGCGCCgcgc                                   >  1:295095/1‑149 (MQ=255)
                                                                                                                          gCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTc                                                  <  1:159184/124‑1 (MQ=255)
                                                                                                                          gCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTc                                                  >  2:159184/1‑124 (MQ=255)
                                                                                                                           cACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTg                        >  2:231911/1‑149 (MQ=255)
                                                                                                                             cTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCa                      >  5:164670/1‑149 (MQ=255)
                                                                                                                             cTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCa                      <  4:88365/149‑1 (MQ=255)
                                                                                                                              ttCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCAc                     >  6:300252/1‑149 (MQ=255)
                                                                                                                                  gCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCCCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAAt                                            <  1:323999/122‑1 (MQ=255)
                                                                                                                                  gCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCACCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAAt                                            >  2:323999/1‑122 (MQ=255)
                                                                                                                                     aCCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAgc                                     <  7:224418/126‑1 (MQ=255)
                                                                                                                                     aCCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAgc                                     >  8:224418/1‑126 (MQ=255)
                                                                                                                                       cgcgGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTTTAGCGa            >  8:194629/1‑149 (MQ=255)
                                                                                                                                              tACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTTTAGCGACCACttt     <  1:231911/149‑1 (MQ=255)
                                                                                                                                                 aGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTTTAGCGACCACTTTCgg  >  2:149295/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTCAAGCGTCGGGACGTTATTTTCAGCGCCAATTTCACGAATGCGCAGCGCGACCAGCCCTGCACCTTTAGCGACCACTTTCGG  >  NC_000913/1965043‑1965336

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