Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,976,529 T→C 9.1% H204R (CAT→CGT) ‡ motA ← proton conductor component of flagella motor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,976,5290TC 9.1% 130.2 / 4.0 44H204R (CAT→CGT) ‡motAproton conductor component of flagella motor
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base C (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCG  >  NC_000913/1976384‑1976672
                                                                                                                                                 |                                                                                                                                               
gTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCAtgtg                                                                                                                                              >  2:552319/1‑149 (MQ=255)
       ttaaGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCACATGCGAtggt                                                                                                                                       <  3:256183/146‑4 (MQ=255)
       ttCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTAGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCACATGCGAtggt                                                                                                                                       <  3:289840/149‑4 (MQ=255)
            aTTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCAc                                                                                                                                  <  2:411249/149‑1 (MQ=255)
             ttAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCAcc                                                                                                                                 >  4:443780/1‑149 (MQ=255)
                gAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAg                                                                                                                              >  2:413973/1‑149 (MQ=255)
                              tttGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGAc                                                                                                                <  4:495359/149‑1 (MQ=255)
                                  aCGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCa                                                                                                                    <  2:6150/141‑1 (MQ=255)
                                  aCGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCa                                                                                                                    >  1:6150/1‑141 (MQ=255)
                                    gCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCgg                                                                                                          >  4:307471/1‑149 (MQ=255)
                                              atTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTaa                                                                                                >  3:483653/1‑149 (MQ=255)
                                               ttttGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAAc                                                                                               >  3:408124/1‑149 (MQ=255)
                                                     tGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGa                                                                                         >  3:19860/1‑149 (MQ=255)
                                                      ggTTTCGGCGCTTTTCTGACATAAAACAGTCGCTAATGGGGAAATAAATCCGTGAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGaa                                                                                        <  1:413973/149‑1 (MQ=255)
                                                            ggCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGAccc                                                                                  >  4:201697/1‑149 (MQ=255)
                                                                   ttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTAc                                                                           <  1:206122/149‑1 (MQ=255)
                                                                      tGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGa                                                                                                                                           >  3:387873/1‑82 (MQ=255)
                                                                      tGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGa                                                                                                                                           <  4:387873/82‑1 (MQ=255)
                                                                                     gCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaaa                                                         <  1:380552/149‑1 (MQ=255)
                                                                                        aaTGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGc                                                      >  2:139297/1‑149 (MQ=255)
                                                                                          tGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGGACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCg                                                    >  2:491164/1‑149 (MQ=255)
                                                                                             ggAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGaa                                                 >  1:223653/1‑149 (MQ=255)
                                                                                                 ataaatCCGTAAGCCGATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTTACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGa                                             <  2:173319/149‑1 (MQ=255)
                                                                                                 ataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGa                                             <  3:307471/149‑1 (MQ=255)
                                                                                                 ataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGa                                             <  1:139297/149‑1 (MQ=255)
                                                                                                  taaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAg                                            >  2:221252/1‑149 (MQ=255)
                                                                                                  taaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAg                                            <  3:201697/149‑1 (MQ=255)
                                                                                                              ccAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGt                                               >  1:151619/1‑134 (MQ=255)
                                                                                                              ccAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGt                                               <  2:151619/134‑1 (MQ=255)
                                                                                                               cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcg                               >  1:379022/1‑149 (MQ=255)
                                                                                                                      aTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                        <  4:408124/149‑1 (MQ=255)
                                                                                                                         ccGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGtt                     <  1:221252/149‑1 (MQ=255)
                                                                                                                          cGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTc                    >  2:419994/1‑149 (MQ=255)
                                                                                                                          cGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTc                    <  4:245098/149‑1 (MQ=255)
                                                                                                                           gAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCg                   >  1:20948/1‑149 (MQ=255)
                                                                                                                             ggAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGGCCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGcc                 <  2:486909/149‑1 (MQ=255)
                                                                                                                                 aGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGa             <  1:419994/149‑1 (MQ=255)
                                                                                                                                     cccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCTGGGACTTc         >  2:487234/1‑149 (MQ=255)
                                                                                                                                      ccaccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATAc                                                             >  4:375072/1‑96 (MQ=255)
                                                                                                                                      ccaccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATAc                                                             <  3:375072/96‑1 (MQ=255)
                                                                                                                                       gcttaTCGCACATGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTg       >  1:57767/5‑149 (MQ=255)
                                                                                                                                       gcttaTCGCACATGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTg       >  1:281839/5‑149 (MQ=255)
                                                                                                                                        accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGc      <  1:552319/149‑1 (MQ=255)
                                                                                                                                            tCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCg  >  1:276196/1‑149 (MQ=255)
                                                                                                                                            tCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCg  >  4:465260/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTAACCGTTCAGATTAGAAAGCAGAGTGACTTTGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCG  >  NC_000913/1976384‑1976672

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