Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,191,696 G→A *355* (TAG→TAA yqiI → fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,191,6960GA97.6% 130.2 / NA 42*355* (TAG→TAAyqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/22);  minor base C (1/0);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCATCC  >  NC_000913/3191563‑3191829
                                                                                                                                     |                                                                                                                                     
ttGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCt                                                                                                                                  <  2:513532/139‑1 (MQ=255)
  gCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTaa                                                                                                                                >  1:336726/1‑139 (MQ=255)
      cttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            <  1:575700/139‑1 (MQ=255)
      cttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            <  1:547477/139‑1 (MQ=255)
       tttttCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            >  1:1170101/1‑138 (MQ=255)
       tttttCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            <  2:1170101/138‑1 (MQ=255)
                 tCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAAc                                                                                                                 >  1:1245154/1‑139 (MQ=255)
                      tGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            >  1:121891/1‑123 (MQ=255)
                      tGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            <  2:121891/123‑1 (MQ=255)
                      tGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGcc                                                                                                            >  2:286820/1‑139 (MQ=255)
                            ttCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                      <  2:336726/139‑1 (MQ=255)
                              cGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGcc                                                                                                    <  1:1387696/139‑1 (MQ=255)
                                         ccATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                         <  1:286820/139‑1 (MQ=255)
                                                      aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                            >  1:77067/1‑139 (MQ=255)
                                                      aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                            >  1:900056/1‑139 (MQ=255)
                                                      aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                            >  2:12987/1‑139 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATTCAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                       >  2:779310/1‑139 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            >  2:1444084/1‑86 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                            <  1:1444084/86‑1 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACt                                                                                                                      >  2:1090435/1‑92 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACt                                                                                                                      <  1:1090435/92‑1 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGt                                                                                        <  2:1158132/122‑1 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGt                                                                                        >  1:1158132/1‑122 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                       <  1:489351/139‑1 (MQ=255)
                                                           aGTGCGAACCCGGCAGAAACGGATGTTGAGGGACCGTTTGGTAAAGACGTAATCATGAATATCACCTATAACTACTATCTAATACACACACTAAAGCGGGCCATCAGACCCGTTATATCAGTGCTCTAACTCCAGctct                                                                       >  1:701164/1‑139 (MQ=255)
                                                               cGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTg                                                                   <  2:1245154/139‑1 (MQ=255)
                                                                     cGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACgg                                                                                                               >  1:177779/1‑89 (MQ=255)
                                                                     cGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACgg                                                                                                               <  2:177779/89‑1 (MQ=255)
                                                                           aaaCGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                       <  1:1261569/139‑1 (MQ=255)
                                                                             aCGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGt                                                     >  1:515171/1‑139 (MQ=255)
                                                                                gttgttGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGa                                                  <  1:611421/139‑1 (MQ=255)
                                                                                              gTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCAt                                                                                                          >  2:851480/1‑69 (MQ=255)
                                                                                              gTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCAt                                                                                                          <  1:851480/69‑1 (MQ=255)
                                                                                                         aCGTAATCCTGAAACTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCg                                                                                                  >  2:893980/1‑66 (MQ=255)
                                                                                                         aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCg                                                                                                  <  1:893980/66‑1 (MQ=255)
                                                                                                              aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGc                                       >  2:1066010/1‑120 (MQ=255)
                                                                                                              aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGc                                       <  1:1066010/120‑1 (MQ=255)
                                                                                                              aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAAcc                      <  1:820524/137‑1 (MQ=255)
                                                                                                              aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAAcc                      >  2:820524/1‑137 (MQ=255)
                                                                                                              aTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATACCAGTGCTCTAACTCCAGCGCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAg                    >  1:164808/1‑139 (MQ=255)
                                                                                                                     aaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTAt             <  2:515171/139‑1 (MQ=255)
                                                                                                                                aaCTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCATcc  <  2:900056/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TTGCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCATCC  >  NC_000913/3191563‑3191829

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: