Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,978,083 +A coding (115/351 nt) flhD ← DNA‑binding transcriptional dual regulator FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,978,0831.A95.0% 61.4 / ‑5.1 20N39Y (AAT→TAT) flhDDNA‑binding transcriptional dual regulator FlhD
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base A (7/12);  minor base C (1/0);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACC  >  NC_000913/1977838‑1978324
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
aaCTCATGATCAAGCGCCTATCTTGCGCCGCGCATCTTCCGGCTGATCAACATAATTCACCCAGCGCGCTGAGAGCCAGATGCCGCTATGAAATTGCTGGAGATCGTAAACGCGGCAATCTTGCGTCAACTGAGTATTCGTCTGGGGGCTGTTTACACGGCAGCGATAAACCAGTTGATTGGTTTCTGCCTGCTTAACCATTTGCGGAAGAGTCAGTGCCGATAACGTTGTCGCCATTTCTTCAttacatc                                                                                                                                                                                                                                               >  2:153748/2‑246 (MQ=255)
       gATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCGTCATTATATGCCGAGAc                                                                                                                                                                                                                                        <  1:91452/251‑1 (MQ=255)
                   ttCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCg                                                                                                                                                                                                                             <  1:248580/250‑1 (MQ=255)
                          gcAGCGCTTCTTCAGGCTGATTAACATCCTTCAGCACGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCttt                                                                                                                                                                                                                     >  2:272004/1‑251 (MQ=255)
                              cgcTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTcc                                                                                                                                                                                                                 >  1:88365/1‑251 (MQ=255)
                                     tcAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTTTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGaa                                                                                                                                                                                                             <  1:214683/248‑1 (MQ=255)
                                                     catTCACCAAGCGTGTTGAGAGCATGATACCGGTATGAATTTGCTGGAGATCGTCAACGCGGGACTCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGCAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCa                                                                                                                                                                                            >  2:201289/1‑249 (MQ=255)
                                                                cGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATg                                                                                                                                                                               >  2:121531/1‑251 (MQ=255)
                                                                 gtgtTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATg                                                                                                                                                                               >  1:106986/1‑250 (MQ=255)
                                                                            atgatgCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGGTGATGt                                                                                                                                                                   <  2:117327/251‑1 (MQ=255)
                                                                                                                        ttGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTcccac                                                                                                                        <  2:245799/250‑1 (MQ=255)
                                                                                                                         tGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAGCGCTGTGCAAGTAGTAAAGATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTccc                                                                                                                          <  2:268319/247‑1 (MQ=255)
                                                                                                                                                                                                         ttGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGc                                       <  2:232714/250‑1 (MQ=255)
                                                                                                                                                                                                               aaGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTGTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATATCCAACTTTATTTTAATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGaa                                <  1:272004/251‑1 (MQ=255)
                                                                                                                                                                                                                  agTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAg                             <  1:121531/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccGCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGAt                      >  1:237162/1‑250 (MQ=255)
                                                                                                                                                                                                                            gCTAACGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCACCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAg                   <  2:106986/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCGTTGTCGCCATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATGAGAATAGTTGCGATAAGCTGc               <  2:22522/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cctttattcATTCTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCACTAAGc        >  2:308586/7‑248 (MQ=255)
                                                                                                                                                                                                                                          cATTTCTTCATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGc        >  1:98040/1‑248 (MQ=255)
                                                                                                                                                                                                                                              tcttcATTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAacc  <  2:79248/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
CACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACC  >  NC_000913/1977838‑1978324

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