Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,192,183 A→G 9.1% F14L (TTT→CTT)  yehD ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,192,1830AG 9.1% 118.6 / 3.7 44F14L (TTT→CTT) yehDputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base G (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

ACACTGGTTGGTACAGTCGGCATATTAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAATGTAAACACATTGTTTTATT  >  NC_000913/2192041‑2192330
                                                                                                                                              |                                                                                                                                                   
acacTGGTTGGTACAGTCGGCATATTAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAaagaag                                                                                                                                                >  2:429073/1‑148 (MQ=255)
          gtACAGTCGGCATATTAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAATGAAGAAAAGACAGc                                                                                                                                      <  5:204032/148‑1 (MQ=255)
                     atatTAATACTTAATGAATAACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCaataat                                                                                                                          <  3:17637/149‑1 (MQ=255)
                      tatTAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAAtt                                                                                                                         <  2:179356/149‑1 (MQ=255)
                       atTAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAAtt                                                                                                                         >  7:265293/1‑148 (MQ=255)
                        ttAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTg                                                                                                                        >  3:175141/1‑148 (MQ=255)
                        ttAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGa                                                                                                                       <  7:354257/149‑1 (MQ=255)
                                   tGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCataa                                                                                                            >  6:40149/1‑149 (MQ=255)
                                     aaTCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCataa                                                                                                            <  6:36262/147‑1 (MQ=255)
                                      aTCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCataatac                                                                                                         <  3:271351/149‑2 (MQ=255)
                                           caccCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTcc                                                                                                    >  8:96590/1‑149 (MQ=255)
                                              ccGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCataata                                                                                                          <  1:63279/140‑1 (MQ=255)
                                              ccGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCataata                                                                                                          >  2:63279/1‑140 (MQ=255)
                                                gCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTtaa                                                                                                >  4:9641/1‑148 (MQ=255)
                                                             acGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTAttttt                                                                                  >  1:181908/1‑149 (MQ=255)
                                                             acGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTAttttt                                                                                  >  4:119603/1‑149 (MQ=255)
                                                                gagaTTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGgt                                                                               >  3:342239/1‑149 (MQ=255)
                                                                   aTTCTTCAATATTCCCCTTAATAGTTAATTTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGgttgt                                                                            <  3:119603/149‑1 (MQ=255)
                                                                       tGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGgttgtt                                                                           >  6:350320/1‑146 (MQ=255)
                                                                       tGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGgttgtt                                                                           <  5:350320/146‑1 (MQ=255)
                                                                        gCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTaa                                                                         >  1:121855/1‑147 (MQ=255)
                                                                         cAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAaga                                                                       <  7:96590/148‑1 (MQ=255)
                                                                            tatTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCt                                                                   <  3:9641/149‑1 (MQ=255)
                                                                              ttCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTg                                                                  <  4:175141/148‑1 (MQ=255)
                                                                                 cccTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTAtc                                                              >  8:275221/1‑149 (MQ=255)
                                                                                       aTAGTTAATGTTCCGGTATCAACGTCTTCAGCAAATACTCCAGCGCTCATAAATAAAGAATAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGCTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTg                                                         <  1:429073/148‑1 (MQ=255)
                                                                                       aTAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTg                                                         <  8:265293/148‑1 (MQ=255)
                                                                                                gTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAAtt                                               >  4:40838/1‑149 (MQ=255)
                                                                                                   ccGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGt                                                                     <  4:58985/124‑1 (MQ=255)
                                                                                                   ccGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGt                                                                     <  6:413177/124‑1 (MQ=255)
                                                                                                   ccGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGt                                                                     >  3:58985/1‑124 (MQ=255)
                                                                                                   ccGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGt                                                                     >  5:413177/1‑124 (MQ=255)
                                                                                                     ggTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTca                                          >  6:164452/1‑149 (MQ=255)
                                                                                                         tCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGaa                                                                                                                      >  8:408997/1‑69 (MQ=255)
                                                                                                         tCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGaa                                                                                                                      <  7:408997/69‑1 (MQ=255)
                                                                                                            aCGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCaa                                    >  4:386330/1‑148 (MQ=255)
                                                                                                                  gcagcaAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGtt                             >  5:192545/1‑149 (MQ=255)
                                                                                                                   cagcaAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGtt                             <  2:181908/148‑1 (MQ=255)
                                                                                                                    agcaAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGc                           <  5:164452/149‑1 (MQ=255)
                                                                                                                      caAATACTCCAGCGCTCATAAAAAGAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCg                         >  8:280630/1‑149 (MQ=255)
                                                                                                                      caAATACTCCAGCGCTCATAAAAAGAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCg                         >  8:273896/1‑149 (MQ=255)
                                                                                                                          tACTCCAGCGCTCATAAAAAGAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAAt                     <  7:273896/149‑1 (MQ=255)
                                                                                                                          tACTCCAGCGCTCATAAAAAGAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAAt                     <  7:280630/149‑1 (MQ=255)
                                                                                                                          tACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAAt                     <  3:401311/149‑1 (MQ=255)
                                                                                                                           aCTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAATg                    >  7:328806/1‑149 (MQ=255)
                                                                                                                                            aaaaGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAATGTAAACACATTGTTTTAt   <  5:40149/149‑1 (MQ=255)
                                                                                                                                             aaaGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAATGTAAACACATTGTTTTAtt  <  2:121855/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ACACTGGTTGGTACAGTCGGCATATTAATACTTACTGAATCACCACCCGCTTCGAATTTACACGGAGATTCTGCAATATTCCCCTTAATAGTTAATGTTCCGGTATCAACGTCTGCAGCAAATACTCCAGCGCTCATAAAAAAAGAAGAAAAGACAGCGGCAGCAATAATTGAACGTTTCATAATATGTTCCTTAATAAAATTTATTTTTGGTTGTTAAGAGTCTGTATCTCTTGGAAACAAATTTCTCACATCAAAAAAGTTGCCGCAATGTAAACACATTGTTTTATT  >  NC_000913/2192041‑2192330

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