Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,672,787 G→T 21.8% M1M (ATG→ATT) † U712_08615 → Flagellar biosynthesis protein flhF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,672,7870GT21.8% 85.5 / 7.7 56M1M (ATG→ATT) †U712_08615Flagellar biosynthesis protein flhF
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base T (10/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 5.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCGCCGCCTGTCAGAATGTATGTGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTGGCCTCGTCAACAAGCC  >  CP006881/1672657‑1672921
                                                                                                                                  |                                                                                                                                       
tGCTCGCCGCCTGTCAGAATGTATGTGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATAGGAGTGGTGGATATTTAATTAAAATaa                                                                                                                                 >  4:361712/1‑138 (MQ=255)
  cTCGCCGCCTGTCAGAATGTATGTGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATaaaa                                                                                                                               >  5:340198/1‑138 (MQ=255)
                gAATGTATGTGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTgctgct                                                                                                                <  2:244912/139‑1 (MQ=255)
                     tatgtGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAAt                                                                                                           >  6:339133/1‑139 (MQ=255)
                     tatgtGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTAATTCTTCTTCAAt                                                                                                           >  6:172686/1‑139 (MQ=255)
                       tgtgAAGCAGCTACTGGACCGGTACTTCCCTGAGGTTCCGGTAGTCTCTTATAACGACCTTGAGGCGAATGTTGAGGTACAGAGCATTGGAGTGGGGGAGTTTTTATGTAATTTAAAACAGTTTCTGGTGGTGCAATg                                                                                                          >  2:53104/1‑138 (MQ=255)
                           aaGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGTTATTTAATTAAAATAAAAAAATTTAATGGTGCTTCaaagaaag                                                                                                      >  3:157308/1‑138 (MQ=255)
                           aaGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGGTTTTTTATGAAAATAAAAAAATTTAATGATGCTTGAATGc                                                                                                         >  5:222515/1‑135 (MQ=255)
                           aaGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCaaga                                                                                                     <  2:52050/139‑1 (MQ=255)
                           aaGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCa                                                                                                        <  4:260271/136‑1 (MQ=255)
                           aaGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCa                                                                                                        >  3:260271/1‑136 (MQ=255)
                                 cTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATTAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGc                                                                                               >  1:286709/1‑139 (MQ=255)
                                 cTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCgg                                                                                                <  5:262111/138‑1 (MQ=255)
                                 cTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCgg                                                                                                >  6:262111/1‑138 (MQ=255)
                                          cGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCtcat                                                                                      <  5:339133/139‑1 (MQ=255)
                                          cGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCtcat                                                                                      <  5:172686/139‑1 (MQ=255)
                                              aCTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGGTATTTTATTAAAATAAAAAAATTTTATGCTGCGTGAATGCAAGGAGCGGCCCGTCtcatcct                                                                                   >  6:309440/1‑136 (MQ=255)
                                                    cTGATCTTACGGTACTCTCTTATAACGAACTGGAGGCGAGTGTTTAGGGACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAagag                                                                            <  6:133861/139‑1 (MQ=255)
                                                      agacagccGTTACTCTATTATAACGAACTTGAGGCTAATGTTGAGGTAAAAAGAATTGGAGTGGTGTATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGtt                                                                          <  4:286984/133‑1 (MQ=255)
                                                       atattgcGGTATTGTCTTATAACGAACTTGAGGCGATTGTTGAGGTACAAAGCATTGTAGTTGTGTATATTTAATTAAATTAAAATAATTTACTGCTTCTTCAATGCAAGTAGCGGCACTTCTCATCAGAAAAGTGtt                                                                          <  6:270984/132‑1 (MQ=255)
                                                       agcTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGTATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGtt                                                                          <  6:41535/136‑1 (MQ=255)
                                                       accTTCCGGTACACTCTTATAACGAACTTCCGGCTAATGTTGCGGTACAAAGCATTGGAGTGGGGGATATTTAATTAAAATAAAAAAATTAACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGtt                                                                          <  6:36637/136‑1 (MQ=255)
                                                       aTCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATTAAAATAAAAAAATTTAATGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTa                                                                         >  3:26578/1‑139 (MQ=255)
                                                       aTCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGtt                                                                          <  3:147567/138‑1 (MQ=255)
                                                       aTCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGtt                                                                          <  3:310883/137‑1 (MQ=255)
                                                       aTCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATATAAAAAAAAAAAAAAAAATTACAGCTGCATCAATGATACAAGCCCCACAACTCGACAGAAAAAAGTTa                                                                         >  5:62809/1‑139 (MQ=255)
                                                            ccGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTa                                                                     <  8:6862/138‑1 (MQ=255)
                                                             cGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTAAAGCTGCTTCAATGCAAGAAGCGGCACGTCTCATCAGAAAAGAGTTAGGTAAt                                                                   >  3:288468/1‑139 (MQ=255)
                                                              ggTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAAt                                                                   >  3:328586/1‑138 (MQ=255)
                                                               gTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGa                                                                 <  6:340198/139‑1 (MQ=255)
                                                                  ctctctTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCTTTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAg                                                               <  3:361712/138‑1 (MQ=255)
                                                                  ctctctTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAg                                                               <  5:359714/138‑1 (MQ=255)
                                                                       ttATATCGAACTTGAGGTGAATGTTGATGTACAAAGCATTGGAGTGGTGGTTTTTTAATGAATATAAAAAAATTTACTGCTGCTTCAATGAAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGt                                                           <  4:157308/137‑1 (MQ=255)
                                                                       ttATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGTTATTTTTTTAAAATAAAAAAATTTTAAGCTGCTTCAATGCAATAACCGGCGCTTCTCATCAGAAAATAGTTAGGTTATGAAGCgggga                                                         >  5:316453/1‑136 (MQ=255)
                                                                       ttATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGTTATTTAAATAAAATAAAAAAAATTTAATCTGCTTAAATGCAAGAAGCGGCACTTCTCAAAATAAAAGAGTTATGTTATGAAGctgg                                                           >  3:286984/1‑134 (MQ=255)
                                                                         ataACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGt                                                           <  4:328586/135‑1 (MQ=255)
                                                                          taACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTAATGCTGCTTCAATGCAAGAAGCGGCACTTGTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCatat                                                      >  2:252679/1‑139 (MQ=255)
                                                                               aaCTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGTTATTTTATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAATCGGGCCTTCTCATCAGATAAGATTTATGTAATTAAGCGGTCAAATTAAAt                                                 >  4:45127/1‑139 (MQ=255)
                                                                               aaCTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTATCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAAt                                                 >  7:296828/1‑139 (MQ=255)
                                                                                 cTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATTAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCAAATCAGAAAATAGTTTGGTAATGAAGCGGTCATATTAAAt                                                 >  1:118995/1‑137 (MQ=255)
                                                                                 cTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATTAAAAAATTTTATTCTGCTTCTATGCAAGAAGCGGCACTTGTGATCTGAAAAGAGTTATGTTATGTAGCGGTCATATTAAATTc                                               >  6:130052/1‑139 (MQ=255)
                                                                                 cTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCTTCAGAAAAGAGTTAGGTAATGAAGCGGGCATATTAAATTc                                               >  7:155423/1‑139 (MQ=255)
                                                                                          aaTGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCAAATTAAATTCGAAAAAAAt                                      >  2:73531/1‑139 (MQ=255)
                                                                                          aaTGTTGAGGTACAAAGCATTGGAGTGGTGGATAATTTATTATTATAAAAAAATTTTATGCAGCTTCAATTCAAGAAGCGGCAAATATCAACATAAAAGACATAGGGAATTAAGCGGTGATATTAAATTAGAAAAAAAt                                      >  1:243898/1‑139 (MQ=255)
                                                                                            tGTTGAGGTACAAAGCATTGGAGTGGTGGGTATTTAATTAAAATAAAAAAATTTTCTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTTATGAAGCGGTCATATTAAATTCGAAAAAAAtt                                     >  5:27762/1‑138 (MQ=255)
                                                                                                      aaaaGCATTGGAGTGGTTGATATTTAATCAAAATAAAAAAATTTAATGATGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGaaa                          <  5:315998/138‑1 (MQ=255)
                                                                                                           cATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGt                     >  5:180279/1‑139 (MQ=255)
                                                                                                              tGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGATTTAGGTAATGAAGCGGGCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTg                  >  2:12080/1‑139 (MQ=255)
                                                                                                              tGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTg                  >  7:179286/1‑139 (MQ=255)
                                                                                                              tGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACAGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTTGGGAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAAAGGTTTg                  >  5:340356/1‑139 (MQ=255)
                                                                                                                caggggtgcATATTTAATGAAAAAAAAAAAATTTACGGCTCCTTCAAGGCAAGCCGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAACGCGAAAATGGTTTGGc                <  6:316453/135‑1 (MQ=255)
                                                                                                                   tggtggATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTGGCCTc             <  1:252679/139‑1 (MQ=255)
                                                                                                                     gtggATATTTTATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGGCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTGGCCGCGt           >  4:20964/1‑139 (MQ=255)
                                                                                                                               aaTGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAagag                                                                            <  1:136636/64‑1 (MQ=255)
                                                                                                                               aaTGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAagag                                                                            >  2:136636/1‑64 (MQ=255)
                                                                                                                               aaTGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTGGCCTCGTCAACAAGcc  <  1:12080/138‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
TGCTCGCCGCCTGTCAGAATGTATGTGAAGCAGCTACTGGAACGGTACTTCCCTGATCTTCCGGTACTCTCTTATAACGAACTTGAGGCGAATGTTGAGGTACAAAGCATTGGAGTGGTGGATATTTAATGAAAATAAAAAAATTTACTGCTGCTTCAATGCAAGAAGCGGCACTTCTCATCAGAAAAGAGTTAGGTAATGAAGCGGTCATATTAAATTCGAAAAAAATTAAAAAGCGAAAATGGTTTGGCCTCGTCAACAAGCC  >  CP006881/1672657‑1672921

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