Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,291,109 A→T 11.5% K766* (AAA→TAA) ‡ yraJ → putative fimbrial usher protein YraJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,1090AT11.5% 79.1 / 5.2 35K766* (AAA→TAA) ‡yraJputative fimbrial usher protein YraJ
Reads supporting (aligned to +/- strand):  ref base A (10/21);  new base T (2/2);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-02

GCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCA  >  NC_000913/3290976‑3291251
                                                                                                                                     |                                                                                                                                              
gCAGCTTCACGATCACGGCGATCTTGAACGACCCACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGc                                                                                                                                       <  1:281683/143‑1 (MQ=255)
  aGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGcga                                                                                                                                     <  2:334987/143‑3 (MQ=255)
             aaCCTCGATCTTGACCAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCTATCGCAACGGTAAACCGTTGCCCTTCGGCGctct                                                                                                                          <  2:600436/143‑1 (MQ=255)
                       tttaacaatcaTCTCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAAc                                                                                                                <  1:744747/140‑1 (MQ=255)
                         gAACAAGCAGCACAGTTTGTGGTGCCCGGCAGTGGTGCAATGGTAGCGGCGCGTTTCAGCGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACga                                                                                                              <  2:792511/143‑1 (MQ=255)
                          aacagacaACACACGGTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCACTATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgat                                                                                                             <  1:343698/143‑1 (MQ=255)
                           acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACCGTGTAGGGTCCCAACCGTAACCCGACTCCCCTCGGCGCTCTTGCCAGCAACgatg                                                                                                            >  2:142625/1‑143 (MQ=255)
                             agacaaccTAGTTTGTGGTGCGCAACAGAGGTGCATTGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCATCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTTTTGCCAGCAACgatgat                                                                                                          <  2:24784/141‑1 (MQ=255)
                                   cacaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTGAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                    <  2:182895/143‑1 (MQ=255)
                                         ttGTGGTGCCCAACAGAGGTGTAATGGTAGCGGCGCGTTTCATCTCCAATATTGGTTATCGTGTATTTGTTACAGTCAGCGAGCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaacaa                                                                                              <  1:69636/143‑1 (MQ=255)
                                          tgtgGTGGCCAGTAGGGGTGCAGGGGGAGCGGCGCGGTTCAACGCCAATATCGGTTGTCGCGTACTTGTTACAGTCAGCGATCGCGAGGGTAAACCGTTGCCCTTTGGCGCGCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                             <  2:327691/143‑1 (MQ=255)
                                          tgtgGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                             <  1:113649/143‑1 (MQ=255)
                                               tgCCCAGCAGCGGGGCAACGGTAGCGGCGCGTTTCAACGGCAATATCCGTTATCGCGTACTTGTTACAGTCAGCTAGCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtc                                                                                        <  1:584503/143‑1 (MQ=255)
                                                 cccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTCGCCCTTTGGCGCTCTTGCCAGCAACGATCATACGGGCCCACAAAGTaaccg                                                                                      >  2:134539/1‑139 (MQ=255)
                                                    aaCAGAGGTGCAATGGTAGCGGCGCGTTTAAACGGCAAGATCGGTTATGGCGTACTGGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                   >  2:485749/1‑143 (MQ=255)
                                                     aCAAAGGTGCAATGGTAGCGGGGCGTTTCAACGCCAATATCGGTTATCGCGGACTTGTTACAGCCAGCGAGCGCAACGCTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGAt                                                                                  <  2:246531/143‑1 (MQ=255)
                                                       ggagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGa                                                                                <  1:774809/142‑1 (MQ=255)
                                                              cAATGGTAGCGGCGCGTTTCAACGCCAATATTGGTTATCGCGTACTTGTTACAGTCCGCGATAGCATCGGTAAGCCGTTGCCATTTGGTGCTCTTGCCAGCAAGAACGCTACGGGGCACCACAATCTCGTGCATGAGggcggc                                                                         >  2:818839/1‑143 (MQ=255)
                                                              cAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCCGCAACGATGATACGGGGCAACACAGTACCGTCGCTGAGggcgga                                                                         >  2:641269/1‑142 (MQ=255)
                                                              cAATGGTAGCGGCGCGATTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATGGCAACGGTAAACCGTCGCCGTTTGGCGCTCTTGCCAGCGACGATGAGACGGGGCACCAAGGTCTCGCCGGTGAGgccggg                                                                         >  2:564770/1‑142 (MQ=255)
                                                                      gcggcgCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtata                                                                 <  1:809847/143‑1 (MQ=255)
                                                                               ttCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTgg                                                        <  1:86234/143‑1 (MQ=255)
                                                                               tccAACGCAAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGTCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTgg                                                        <  1:520579/141‑1 (MQ=255)
                                                                                         atatCGGTTATCGCGCACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGCCGCTCTTGCCATCAATGATGCTACGGGGCAACAAAGTATCGTCGATGAGCGCGGCATACTATTTCTCGCAGGGTTATCCAGt                                              >  2:822327/1‑143 (MQ=255)
                                                                                               gTTATCGCGTACCTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATATGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATca                                        <  1:453140/143‑1 (MQ=255)
                                                                                               gTTATCGCGTACCTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGACGATACGGGGCAACACAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGTATATCGAGTAAATca                                        <  1:459961/143‑1 (MQ=255)
                                                                                                         aCTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGa                              <  1:338068/143‑1 (MQ=255)
                                                                                                          ttggtgaCAGTCAGCGATCGCAACGGTAACCCGTTGCGCTTTGGCGCTCTGGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGAc                             <  1:55784/137‑1 (MQ=255)
                                                                                                                       gcgcTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGa                <  1:512680/143‑1 (MQ=255)
                                                                                                                             gCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAGGCTGGACTGTACGCTGGGGaagtagg          >  2:43817/1‑138 (MQ=255)
                                                                                                                             gCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAg          >  2:164886/1‑143 (MQ=255)
                                                                                                                             gCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAg          <  1:477652/143‑1 (MQ=255)
                                                                                                                             gCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGGAAAGCACAAAGCTGGACTGTACGCTGGGGAAAGCAg          >  2:691095/1‑143 (MQ=255)
                                                                                                                              cAACGGTAAACCGTTGCCCTTTGGCGCTCTTGACAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGGGGGCGGCATACTATATCTCTCTGGGATACCGAGTAAATCACAAAGCTGGACTGGACTCGAGGGaaaaaagg         >  2:463654/1‑138 (MQ=255)
                                                                                                                               aaCGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGCCCACAAAGTATCGTCGATGAGGGCGGCATACTATAACTCACTGGGATATCGAGTAAAACACAAAGCGGGACAGTACGCGGGGGGCATCACGc        >  2:270738/1‑143 (MQ=255)
                                                                                                                                     aaaCCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGCTATCGAGTAAATCACAAAGCTGACCTGCACGCTGGCGGAATCAAGCAGATca  >  2:123341/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                              
GCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCA  >  NC_000913/3290976‑3291251

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: