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,6910AG93.8% 45.4 / ‑2.2 16intergenic (‑237/‑543)flhD/uspCflagellar transcriptional regulator FlhD/universal stress protein UspC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (7/8);  minor base . (1/0);  total (8/8)
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.65e-01

TTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTG  >  NZ_CP009273/1972604‑1972776
                                                                                        |                                                                                     
ttACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGgggg                                                                                     >  2:305161/1‑90 (MQ=255)
         gAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCg                                                                            >  2:397925/1‑90 (MQ=255)
                     tACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataat                                                                >  2:44242/1‑90 (MQ=255)
                      acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATg                                                               <  1:10314/90‑1 (MQ=255)
                                      tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcaca                                               >  1:312104/1‑90 (MQ=255)
                                                  aCGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCaa                                   <  1:349146/90‑1 (MQ=255)
                                                    gCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAtt                                 <  1:151586/90‑1 (MQ=255)
                                                        cTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACtt                             >  2:108889/1‑90 (MQ=255)
                                                        cTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAAAAATGCGTGATGCAGATCACACAAAACACTCAATTACtt                             <  2:453590/90‑1 (MQ=255)
                                                               cGTCACTACACGCACATACAACGGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                       <  1:359114‑M1/90‑5 (MQ=255)
                                                               cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            >  2:319540/1‑52 (MQ=255)
                                                               cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                            <  1:319540/52‑1 (MQ=255)
                                                                 tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                    >  1:339251‑M1/1‑83 (MQ=255)
                                                                    cTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcga                >  1:353952‑M1/1‑79 (MQ=255)
                                                                              aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  1:51163‑M1/90‑21 (MQ=255)
                                                                              aTACAACGGGGGGGGGCTGCGATTTTCAATAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  2:359105‑M1/90‑21 (MQ=255)
                                                                                   aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata  >  1:45928‑M1/1‑65 (MQ=255)
                                                                                        |                                                                                     
TTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTG  >  NZ_CP009273/1972604‑1972776

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