Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,966,120 C→T G24G (GGG→GGA)  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,966,1200CT98.0% 192.1 / ‑4.1 50G24G (GGG→GGA) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/29);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGG  >  NC_000913/1965972‑1966264
                                                                                                                                                    |                                                                                                                                                
ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGtcc                                                                                                                                                <  1:399901/151‑1 (MQ=255)
    aTTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCctct                                                                                                                                            >  2:609276/1‑151 (MQ=255)
          aataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCcg                                                                                                                                       >  1:349147/1‑150 (MQ=255)
            taCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCgcgc                                                                                                                                    >  1:610468/1‑151 (MQ=255)
                gtgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCtt                                                                                                                                  <  1:195780/149‑1 (MQ=255)
                  gATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTc                                                                                                                              <  2:499426/151‑1 (MQ=255)
                     cAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTa                                                                                                                            <  2:343831/150‑1 (MQ=255)
                            cAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGa                                                                                                                    <  2:174192/151‑1 (MQ=255)
                                  aGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGtgtg                                                                                                              >  2:279867/1‑151 (MQ=255)
                                  aGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGtgtg                                                                                                              <  2:159367/151‑1 (MQ=255)
                                     ggAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgggg                                                                                                           >  2:234468/1‑151 (MQ=255)
                                      gagaGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTggggg                                                                                                          >  1:290342/1‑151 (MQ=255)
                                         agCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCtt                                                                                                       >  1:558118/1‑151 (MQ=255)
                                               cccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt                                                                                                 <  2:168417/151‑1 (MQ=255)
                                               cccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGt                                                                                                   <  2:319601/149‑1 (MQ=255)
                                                ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGtttt                                                                                                >  2:104743/1‑151 (MQ=255)
                                                     aaTCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTgtcgt                                                                                           <  1:157222/151‑1 (MQ=255)
                                                             ggCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCg                                                                                   <  1:495110/151‑1 (MQ=255)
                                                                 aGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTc                                                                                 >  1:412728/1‑149 (MQ=255)
                                                                    gACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGaca                                                                            <  1:219671/151‑1 (MQ=255)
                                                                    gACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGac                                                                             >  1:307776/1‑150 (MQ=255)
                                                                     acacacCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGacac                                                                           <  1:471852/151‑1 (MQ=255)
                                                                        cacCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTc                                                                        <  1:279874/151‑1 (MQ=255)
                                                                           cGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGCGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGCGCGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGcc                                                                     >  2:179409/1‑151 (MQ=255)
                                                                               aaaCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATc                                                                 >  1:215219/1‑151 (MQ=255)
                                                                                aaCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATcc                                                                >  2:591309/1‑151 (MQ=255)
                                                                                      ataataCTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGaaa                                                          <  2:198834/151‑1 (MQ=255)
                                                                                        aaCACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAaga                                                        <  2:147156/151‑1 (MQ=255)
                                                                                         acacTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                         >  2:302150/1‑149 (MQ=255)
                                                                                             tAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTAtt                                                   >  1:47242/1‑151 (MQ=255)
                                                                                             tAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCCGGAAAGAGTAtt                                                   >  2:171784/1‑151 (MQ=255)
                                                                                              aaCGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTAt                                                    <  2:91126/149‑1 (MQ=255)
                                                                                               aCGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTaa                                                 <  1:249271/151‑1 (MQ=255)
                                                                                                cGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAg                                                <  2:519031/151‑1 (MQ=255)
                                                                                                cGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAg                                                <  2:388292/151‑1 (MQ=255)
                                                                                                 gCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGc                                               <  1:195308/151‑1 (MQ=255)
                                                                                                       caGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAat                                         <  1:437902/151‑1 (MQ=255)
                                                                                                          cAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAatgat                                      >  2:550229/1‑151 (MQ=255)
                                                                                                              aTCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCa                                  <  1:595655/151‑1 (MQ=255)
                                                                                                                cagcagTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCaga                                <  2:19614/151‑1 (MQ=255)
                                                                                                                 agcagTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCagag                               <  2:270498/151‑1 (MQ=255)
                                                                                                                  gcagTGAGGTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAgc                              <  1:212582/151‑1 (MQ=255)
                                                                                                                         ggTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGc                       <  2:266328/151‑1 (MQ=255)
                                                                                                                          gTCAGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCg                      >  2:42575/1‑151 (MQ=255)
                                                                                                                             aGTTCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTc                   >  1:342237/1‑151 (MQ=255)
                                                                                                                                tCACGGGAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCaa                 <  1:586520/150‑1 (MQ=255)
                                                                                                                                     ggAACGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGttt           <  2:477010/151‑1 (MQ=255)
                                                                                                                                       aaCGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTg          >  1:8047/1‑150 (MQ=255)
                                                                                                                                        aCGCGGGATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAg        <  1:233588/151‑1 (MQ=255)
                                                                                                                                             ggATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAgg   <  2:307769/151‑1 (MQ=255)
                                                                                                                                              gATTTGTCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAggg  >  1:455333/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGG  >  NC_000913/1965972‑1966264

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑