Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,019,460 G→A V359I (GTC→ATC)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,019,4600GA97.9% 129.5 / ‑5.7 47V359I (GTC→ATC) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/21);  total (26/21)
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

CATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACA  >  NC_000913/2019321‑2019602
                                                                                                                                           |                                                                                                                                              
cATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCGCAGCAACTGCATAGCCAACGCTCAGCATTCCTTGAACCTCTGGCGGGGGAAGACGACGATTCGCTTCCGGTTCCCATCTCTTTAc                                                                                                                                       <  2:193482/149‑1 (MQ=255)
 aTTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAGCCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACa                                                                                                                                      >  2:440156/1‑149 (MQ=255)
 aTTCAGTTAGGGCAAAGCAACAGCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCCACAGCAAAGCCACCGCACAGCAAACCATGAACCGCTGGCGGGGGAAGCCGACGCTCCGCTTCCGGGTCCCATCTCTGTAcc                                                                                                                                      >  2:1138014/1‑148 (MQ=255)
    cAGTTAGGGCAAAGCAACATCAGTGGAGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAgg                                                                                                                                   <  1:205199/149‑1 (MQ=255)
      tgtatgtcAAAGAAACATCAGTGGCGAAAGCTTTAGTGGTCAGACGCAGGCCGCTTCCCAGCAAAAGCAAAGCCAACGCACAGCAAACAATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGc                                                                                                                                 <  2:771867/142‑1 (MQ=255)
       ttAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCg                                                                                                                                >  1:162452/1‑149 (MQ=255)
       ttAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCAAAGAAAACCATTAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTACCATCTCTTTACAAGGGCg                                                                                                                                >  1:561697/1‑149 (MQ=255)
            gCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTaa                                                                                                                           >  1:66650/1‑149 (MQ=255)
            gCAAAGCAACATCAGTGGCGAAAGCGTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTaa                                                                                                                           <  2:900796/149‑1 (MQ=255)
             cAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCACCGCACAGCAATCCATGACCCTCTGGCGGGGTCAGACGACGATCCGCTTCCGGTTCCCATCTCTTTACAACGGCGTGGAAc                                                                                                                          >  2:87490/1‑149 (MQ=255)
              aaaGCAACATCAGTGGCGAACGCTTTAGTGGTCAGCAGCAGGCCGCTTCACAGCAACAGCAAAGCCCACGCTCCGCAAACCATGAACCGCTGGCGGGGGAAGTCGACGATACGCTTCCTGTTCCCAACTCTTTACAACGGCGTGAAAcc                                                                                                                         >  2:804282/1‑148 (MQ=255)
              aaaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACa                                                                                                                         <  2:66650/149‑1 (MQ=255)
               aaGCAACATCAGTGGGGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAg                                                                                                                        <  2:1031605/149‑1 (MQ=255)
                             gCGAAAGCTTTAGTGGTCAGCAGCAGTCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                          <  1:1224564/149‑1 (MQ=255)
                             gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGCTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                          >  1:173983/1‑149 (MQ=255)
                             gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATGTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                          >  1:843939/1‑149 (MQ=255)
                               gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                                                                                        >  1:424008/1‑149 (MQ=255)
                                aaaGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAATGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGGGTAACCGGCAACAGCGGCGTTGat                                                                                                       >  2:396404/1‑149 (MQ=255)
                                aaaGCTTTAGGGGGCAGCAGCAGGCCGCTGCCCAGCAACAGCAAAGCCACCGCACCGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                                                                                       <  1:722959/149‑1 (MQ=255)
                                 aaGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATCCGCTTCCGGTTCCCGTCTCTTTACCAGGGCGTGTAACTGGCAACAGCGGCGTGGatt                                                                                                      >  1:425720/1‑148 (MQ=255)
                                     tttAGGGGTCCGCAGCAGGCCGCTTCCCAGCTACAGCATAGCCAACGCAGAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCATCCGGTTCCCATCTCTTTACAAGCGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                                                                                  <  2:441658/149‑1 (MQ=255)
                                      ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                 >  2:625418/1‑149 (MQ=255)
                                       tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCg                                                                                                <  1:45078/149‑1 (MQ=255)
                                       tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCGTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCACCAGCGGCGATGATATTTTCg                                                                                                >  2:311401/1‑149 (MQ=255)
                                        aGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                                                                                               <  1:595973/149‑1 (MQ=255)
                                           ggTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa                                                                                            >  1:864075/1‑149 (MQ=255)
                                                                   caacagcaaAGCCAACGCACCGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACCAGGGCGTGTAACAGGCAACAGCGGCGTTGAAATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc                                                                    >  2:99221/1‑149 (MQ=255)
                                                                     acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGCTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCcgc                                                                  >  2:510166/1‑149 (MQ=255)
                                                                     acagacaaGCGAACGCACAGCAAACCCTGCACCTCTGGCGGGGGACGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTATCAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCcgc                                                                  <  1:1029884/143‑1 (MQ=255)
                                                                          aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCtt                                                             >  2:783199/1‑149 (MQ=255)
                                                                            aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCGTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCtttt                                                           >  1:70992/1‑149 (MQ=255)
                                                                                aaCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCATAATCCGCGTCTTTTcccc                                                       <  2:424008/149‑1 (MQ=255)
                                                                                     acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCttt                                                  >  1:728515/1‑149 (MQ=255)
                                                                                        gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                               <  2:791877/149‑1 (MQ=255)
                                                                                          aaaCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTACCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACAGTAATCCGCGTCTTTTACCCGCTTTGTTgc                                             >  1:709593/1‑149 (MQ=255)
                                                                                             ccATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCt                                          >  1:1175278/1‑149 (MQ=255)
                                                                                                tGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCaa                                       <  1:604164/149‑1 (MQ=255)
                                                                                                tGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCaa                                       <  2:725402/149‑1 (MQ=255)
                                                                                                      tctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCa                                 <  1:783199/149‑1 (MQ=255)
                                                                                                       ctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTATGTTGCGCTCAAGACGCAg                                >  2:1061812/1‑149 (MQ=255)
                                                                                                             ggggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAt                          <  1:954087/149‑1 (MQ=255)
                                                                                                               ggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAttt                        >  1:1113234/1‑148 (MQ=255)
                                                                                                                     cgacgaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGAt                  <  2:728515/149‑1 (MQ=255)
                                                                                                                         gaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGc              <  1:440156/149‑1 (MQ=255)
                                                                                                                              gCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCCTCTTTACCCCGCTTTGTTTCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTac         >  1:618816/1‑148 (MQ=255)
                                                                                                                               cTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGc        >  1:36163/1‑149 (MQ=255)
                                                                                                                                     gTTCCCATCTCTATACAATGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGAcaca  <  1:86605/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACA  >  NC_000913/2019321‑2019602

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