Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,972,691 A→G intergenic (‑237/‑543) flhD ← / → uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,972,6910AG82.9% 69.2 / 10.7 36intergenic (‑237/‑543)flhD/uspCflagellar transcriptional regulator FlhD/universal stress protein UspC
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (11/18);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTG  >  NZ_CP009273/1972605‑1972776
                                                                                      |                                                                                     
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                    >  2:278704/1‑90 (MQ=255)
     aTAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCt                                                                               >  2:197093/1‑90 (MQ=255)
         aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGa                                                                           <  2:207501/90‑1 (MQ=255)
            tGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGAttt                                                                        <  2:55103/90‑1 (MQ=255)
            tGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGAttt                                                                        >  2:384467/1‑90 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGCCGCAATTGTTCTCGTCACTACCCGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            <  2:638483/90‑1 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            >  1:218618/1‑90 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            <  2:24721/90‑1 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            <  1:387023/90‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                           >  1:49518/1‑90 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            <  1:60904/88‑1 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            >  2:60904/1‑88 (MQ=255)
                                   tttctCCTAAATCGACGCAACCGTACTCGTCACTACCCGCCCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                                 <  2:51475/90‑1 (MQ=255)
                                    ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcac                                                >  2:326038/1‑90 (MQ=255)
                                        ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaa                                            <  2:398951/90‑1 (MQ=255)
                                           aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAac                                         <  1:601389/90‑1 (MQ=255)
                                             aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAacac                                       <  2:704474/90‑1 (MQ=255)
                                                   gCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAtt                                 <  1:197093/90‑1 (MQ=255)
                                                       cTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACtt                             <  1:384467/90‑1 (MQ=255)
                                                             tCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                       <  1:646689‑M1/90‑5 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAt                                                   >  1:603791/1‑61 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAt                                                   <  2:603791/61‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                      >  1:585384‑M1/1‑85 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                    >  2:644360‑M1/1‑83 (MQ=255)
                                                                tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                    >  1:800861‑M1/1‑83 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtg                   >  2:186281‑M1/1‑82 (MQ=255)
                                                                 cACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtg                   <  1:534856‑M1/90‑9 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                 >  1:354620‑M1/1‑80 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                 >  1:516449‑M1/1‑80 (MQ=255)
                                                                      cacGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac              <  2:450530‑M1/90‑14 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  1:507404‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  2:738641‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  2:415453‑M1/90‑21 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata  <  1:278704‑M1/90‑26 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata  >  2:574994‑M1/1‑65 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata  >  2:219849‑M1/1‑65 (MQ=255)
                                                                                  aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggggcgaacaagtccctgata  >  1:625612‑M1/1‑65 (MQ=255)
                                                                                      |                                                                                     
TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTG  >  NZ_CP009273/1972605‑1972776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: