Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,133,668 C→T V365V (GTC→GTT)  flgE → flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,133,6680CT94.7% 49.4 / ‑3.3 20V365V (GTC→GTT) flgEflagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (9/9);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

CAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAG  >  NC_000913/1133422‑1133872
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                            
cAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCACTCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGAtc                                                                                                                                                                                                          >  1:168092/1‑251 (MQ=255)
                            atcaaTGATGACGGTACGGTGGTCGGCACCTAGTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAATGCGACAACGTCTGGTCTGCGACGCAAGCTTCGGGCGTGGCGCTGTTTGTGACAGCCTGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGTTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGt                                                                                                                                                                              <  1:337951/251‑1 (MQ=255)
                               aaTGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCGGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAATCGTTGGCTCCCTGACCAACGGGGCGCTGGAAGCGTCGAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGc                                                                                                                                                                           <  1:426647/251‑1 (MQ=255)
                                atgatgACGGTACGGTTGTCGGCAACGATTCCAACGAACAAACCCAACTGCTGGGGCAGAGTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAAGCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGcc                                                                                                                                                                          <  1:377157/251‑1 (MQ=255)
                                       nggtacggtTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTATGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGt                                                                                                                                                                    >  1:89351/2‑250 (MQ=255)
                                                     gCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGAGTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTTGCTCCCTGACCAACGGTGCGCTGGAATCGTCCATGGTGGATCTGATTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAAc                                                                                                                                                     <  1:254745/251‑1 (MQ=255)
                                                                            cAACTGCTGGGGCAGATTGTACTGGAGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGAt                                                                                                                                                                                 >  1:346668/1‑200 (MQ=255)
                                                                                 gctgGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACACCGTATGGCCTGCGACGCAATCTTCTGGCGTGGAGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGAcc                                                                                                                            >  1:178773/1‑248 (MQ=255)
                                                                                      ggCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACACCGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGATGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATAAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTc                                                                                                                    >  1:438857/1‑251 (MQ=255)
                                                                                                        ncTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTc                                                                                                                                                         >  1:266709/2‑196 (MQ=255)
                                                                                                                  caacGAAGGTCTGGCATCCGAAGGCGCCAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGATGGGAAACTTTGGCACCCTGACCAGCGGTGCGCTGGAAGCGTCCAACGTTGATCTGAGTAAAGAACTGGTCGGTATGATCGTGGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGc                                                                                        <  1:290583/251‑1 (MQ=255)
                                                                                                                   aacGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTATGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTg                                                                                                                                                                            >  1:2538/1‑166 (MQ=255)
                                                                                                                                                           cGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGc                                                                            <  1:61523/222‑1 (MQ=255)
                                                                                                                                                            gANNNNNTCTNCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTc                                                                                                                                                         <  2:266709/144‑1 (MQ=255)
                                                                                                                                                                                      tNNGGACAGNNGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTg                                                                                                                                                                            <  2:2538/99‑1 (MQ=255)
                                                                                                                                                                                      tGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAAGTTATACCGCGATggg                                            <  1:374026/227‑1 (MQ=255)
                                                                                                                                                                                                            cTTTGNCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGTAAAGAACNGGTCAATATNNTCNTTGCCNANNNTNNNTANCAGTCTAACGCCCAGACCANCANNACCCAGNACCAGATCCTCAACACGCTGGTTAAc                                                                                                     <  2:97890/148‑1 (MQ=255)
                                                                                                                                                                                                            cTTTGGAACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTTGATCTCAGGAAAGAACTGGTCAATAAGATCGATGCCCCGCGAAACTATCAGTCTAACGCACAGACCATCAAAAACCAGGACCAGATCCTCAAAACGCTGGTTAAc                                                                                                     >  1:97890/1‑148 (MQ=255)
                                                                                                                                                                                                                                       ggAAGCGTCTAACGTTGATCTCAGTAAAGAACTTGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACg                                                                                                                                                    <  1:303439/74‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaCGTTGATCTCAGTAAAGAACTGGTCAATATGATCGCTGCCCAGCGTAACTATAAGTCTACCGCCCAGCCCATCAAACCCCAGGACCAGATCCTCAACAAGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAg  >  1:228360/1‑210 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                            
CAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATTGTACTGGCGAACTTTGCCAACAACGAAGGTCTGGCATCCGAAGGCGACAACGTCTGGTCTGCGACGCAATCTTCTGGCGTGGCGCTGTTGGGGACAGCCGGGACGGGAAACTTTGGCACCCTGACCAACGGTGCGCTGGAAGCGTCCAACGTCGATCTCAGTAAAGAACTGGTCAATATGATCGTTGCCCAGCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAG  >  NC_000913/1133422‑1133872

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