Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,022,885 G→A *90* (TAG→TAA fliQ → flagellar biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,022,8850GA96.9% 195.6 / ‑1.9 65*90* (TAG→TAAfliQflagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (32/31);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

CGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGGT  >  NC_000913/2022750‑2023017
                                                                                                                                       |                                                                                                                                    
cGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCg                                                                                                                                   <  1:618531/139‑1 (MQ=255)
   aaTGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTAc                                                                                                                                <  2:417294/139‑1 (MQ=255)
    aTGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACt                                                                                                                               <  2:600403/139‑1 (MQ=255)
     tGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCg                                                                                                                                   <  1:879605/134‑1 (MQ=255)
     tGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCg                                                                                                                                   >  2:879605/1‑134 (MQ=255)
     tGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTa                                                                                                                              >  2:327950/1‑139 (MQ=255)
     tGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTa                                                                                                                              >  1:442977/1‑139 (MQ=255)
          cTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTg                                                                                                                         <  2:1036481/139‑1 (MQ=255)
                     tCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAg                                                                                                              <  1:920293/139‑1 (MQ=255)
                       cGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCg                                                                                                            >  1:588598/1‑139 (MQ=255)
                         aaaaTCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGaa                                                                                                          <  2:175658/139‑1 (MQ=255)
                                    cGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               <  1:257833/139‑1 (MQ=255)
                                          tATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGtt                                                                                         <  2:738508/139‑1 (MQ=255)
                                                 attattattGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTg                                                                                  <  2:588598/139‑1 (MQ=255)
                                                  ttattattGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGt                                                                                 >  1:1064033/1‑139 (MQ=255)
                                                          gCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGt                                                                                 <  2:926824/131‑1 (MQ=255)
                                                          gCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGt                                                                                 >  1:926824/1‑131 (MQ=255)
                                                          gCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTgg                                                                         >  1:238333/1‑139 (MQ=255)
                                                              gACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               >  2:759561/1‑113 (MQ=255)
                                                              gACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               <  1:759561/113‑1 (MQ=255)
                                                              gACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGt                                                                     >  2:1215274/1‑139 (MQ=255)
                                                                  gTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACt                                                                 >  2:133144/1‑139 (MQ=255)
                                                                  gTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACt                                                                 >  1:276235/1‑139 (MQ=255)
                                                                         cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTa                                                                                                                                 >  2:199173/1‑68 (MQ=255)
                                                                         cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTa                                                                                                                                 <  1:199173/68‑1 (MQ=255)
                                                                         cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGt                                                                                                                    <  1:1008759/81‑1 (MQ=255)
                                                                         cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGt                                                                                                                    >  2:1008759/1‑81 (MQ=255)
                                                                               cTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                                    <  1:1215274/139‑1 (MQ=255)
                                                                               cTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                                    <  1:133144/139‑1 (MQ=255)
                                                                               cTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                                    <  1:327950/139‑1 (MQ=255)
                                                                                tGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcgc                                                   <  2:442977/139‑1 (MQ=255)
                                                                                  ttGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTg                                                 >  1:144825/1‑139 (MQ=255)
                                                                                   tGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               <  1:1028630/92‑1 (MQ=255)
                                                                                   tGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               >  2:1028630/1‑92 (MQ=255)
                                                                                        gATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTgcg                                                              <  2:215344/120‑1 (MQ=255)
                                                                                        gATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTgcg                                                              >  1:215344/1‑120 (MQ=255)
                                                                                            aCGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGCCAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACcg                                       >  2:448109/1‑139 (MQ=255)
                                                                                              gTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACcgcg                                     <  2:238333/139‑1 (MQ=255)
                                                                                                   cACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACaa                                                                                                       >  2:698409/1‑68 (MQ=255)
                                                                                                   cACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACaa                                                                                                       <  1:698409/68‑1 (MQ=255)
                                                                                                        tGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGt                                                                                                                    <  1:295050/50‑1 (MQ=255)
                                                                                                        tGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGt                                                                                                                    >  2:295050/1‑50 (MQ=255)
                                                                                                        tGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGa                           <  1:347935/139‑1 (MQ=255)
                                                                                                        tGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGa                           >  1:916770/1‑139 (MQ=255)
                                                                                                           tCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCg                        >  1:151190/1‑139 (MQ=255)
                                                                                                               tAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACg                    <  2:1082977/139‑1 (MQ=255)
                                                                                                                  ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTg                                                                                  <  1:102254/74‑1 (MQ=255)
                                                                                                                  ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTg                                                                                  >  2:102254/1‑74 (MQ=255)
                                                                                                                  ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTg                            >  1:545494/1‑128 (MQ=255)
                                                                                                                  ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTg                            <  2:545494/128‑1 (MQ=255)
                                                                                                                  ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAg                 >  2:1136685/1‑139 (MQ=255)
                                                                                                                    tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTa                                                                                        <  1:158978/66‑1 (MQ=255)
                                                                                                                    tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTa                                                                                        >  2:158978/1‑66 (MQ=255)
                                                                                                                    tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGTGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCa                  >  2:284837/1‑136 (MQ=255)
                                                                                                                    tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGTGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCa                  <  1:284837/136‑1 (MQ=255)
                                                                                                                    tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCg               >  2:166300/1‑139 (MQ=255)
                                                                                                                       cGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               >  1:18448/1‑56 (MQ=255)
                                                                                                                       cGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                               <  2:18448/56‑1 (MQ=255)
                                                                                                                         tataTCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCAc                                         <  2:69308/108‑1 (MQ=255)
                                                                                                                         tataTCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCAc                                         >  1:69308/1‑108 (MQ=255)
                                                                                                                             tcatcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATgg                                                                                                    >  2:450369/1‑45 (MQ=255)
                                                                                                                             tcatcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATgg                                                                                                    <  1:450369/45‑1 (MQ=255)
                                                                                                                                tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg   >  2:137335/1‑139 (MQ=255)
                                                                                                                                tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg   >  1:413820/1‑139 (MQ=255)
                                                                                                                                 cGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGGt  >  2:48852/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                    
CGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGGT  >  NC_000913/2022750‑2023017

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