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,6910AG98.0% 160.9 / ‑3.7 52intergenic (‑237/‑543)flhD/uspCflagellar transcriptional regulator FlhD/universal stress protein UspC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (27/23);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACT  >  NZ_CP009273/1972605‑1972772
                                                                                       |                                                                                 
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                >  2:1327177/1‑90 (MQ=255)
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                >  1:122503/1‑90 (MQ=255)
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                >  2:1000456/1‑90 (MQ=255)
      tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTg                                                                          >  1:650257/1‑90 (MQ=255)
       aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGc                                                                         <  1:901577/90‑1 (MQ=255)
       aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGc                                                                         <  1:1560786/90‑1 (MQ=255)
              ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACTCGCACATACAACGGGGGGGGGCTGCGATTTTc                                                                  <  2:445342/90‑1 (MQ=255)
              ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTc                                                                  >  2:1038045/1‑90 (MQ=255)
              ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTc                                                                  <  1:228505/90‑1 (MQ=255)
                  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaata                                                              >  1:587368/1‑90 (MQ=255)
                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                             >  1:110897/1‑90 (MQ=255)
                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                             >  1:125028/1‑90 (MQ=255)
                     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATg                                                           <  1:180216/90‑1 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        >  2:528673/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCGCACATACAACGGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        <  1:1099852/90‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        <  2:1551161/89‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        >  1:1551161/1‑89 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                       >  2:172178/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                       >  2:1051709/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                       >  1:1376083/1‑90 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGa                                                      >  2:777503/1‑90 (MQ=255)
                           aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGAt                                                     >  1:1509431/1‑90 (MQ=255)
                                   ttttttCTAAATCGACGCAACTGTTCTCGTCACCACCCCGCCATACCACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                             <  1:1449381/90‑1 (MQ=255)
                                     tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcaca                                           <  2:1023057/90‑1 (MQ=255)
                                             aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAacac                                   <  1:68579/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                              <  1:159688/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                              <  2:1239950/90‑1 (MQ=255)
                                                   gCAACTGTTCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAAtt                             <  2:1373297/90‑1 (MQ=255)
                                                         gTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTaa                       >  2:76730/1‑90 (MQ=255)
                                                           aCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                     >  1:1499088‑M1/1‑88 (MQ=255)
                                                           aCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                     >  1:839213‑M1/1‑88 (MQ=255)
                                                            cTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                    <  2:1456717‑M1/90‑4 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaa                                        >  1:1356699/1‑68 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaa                                        <  2:1356699/68‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                  >  2:875710‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                  >  1:722095‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                 >  2:506573‑M1/1‑84 (MQ=255)
                                                               gTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAGTGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                 <  2:371067‑M1/90‑7 (MQ=255)
                                                               gTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                 <  1:510021‑M1/90‑7 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                             <  1:1116453/61‑1 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                             >  2:1116453/1‑61 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                >  1:430875‑M1/1‑83 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg             >  2:840751‑M1/1‑80 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        >  1:983682‑M1/1‑75 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        >  2:1571058‑M1/1‑75 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   >  2:1264749‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   >  2:1099858‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   <  2:722095‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   <  1:1038045‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   >  1:557014‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   <  2:86041‑M1/90‑21 (MQ=255)
                                                                              tACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct  <  1:342933‑M1/90‑22 (MQ=255)
                                                                              tACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct  <  2:994423‑M1/90‑22 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        >  2:311011‑M1/1‑65 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        <  1:311011‑M1/80‑16 (MQ=255)
                                                                                       |                                                                                 
TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACT  >  NZ_CP009273/1972605‑1972772

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