Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 154,182 C→T L415L (CTG→CTA htrE ← putative fimbrial usher protein HtrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913154,1820CT92.6% 71.3 / ‑3.2 27L415L (CTG→CTAhtrEputative fimbrial usher protein HtrE
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (19/6);  total (20/7)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCAT  >  NC_000913/153967‑154414
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
ggATGTTTAACTTCATCCATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGGACTTGTTTCTTCGAATCACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGCTCCGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCACGACCTAACAAACCAGCGGAATAGttatt                                                                                                                                                                                                       >  2:429737/1‑251 (MQ=255)
       tCACGTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGt                                                                                                                                                                                                 <  2:239467/250‑1 (MQ=255)
              atcaATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGa                                                                                                                                                                                          >  1:298813/1‑250 (MQ=255)
                 aaTTAGAGTTAGTGCATCATTAAGACCACGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATAc                                                                                                                                                                                      >  1:103221/1‑251 (MQ=255)
                                       aGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGtt                                                                                                                                                                <  1:498199/251‑1 (MQ=255)
                                       aGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGtt                                                                                                                                                                >  1:449017/1‑251 (MQ=255)
                                               gTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCGTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGc                                                                                                                                                        <  2:89254/251‑1 (MQ=255)
                                                          tCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCgtagtagta                                                                                                                                              <  1:455637/250‑1 (MQ=255)
                                                                aGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGGTAACTTTGCCCATGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGAAGTAGTAgctttct                                                                                                                                       >  1:111519/1‑248 (MQ=255)
                                                                             gATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAAATTTGCCCCTGGTATGTTATATCATCTGGGATACGAACATTGGAATGAGTCCCATCGAAAGAAAATGCACCAACTGAAGCATTTAGACCAAAACCTAACAAACCAGCGGTATAGATATTATCGGTTCTCTGAATACCGGTATACACCGTCAGATAGTTATTCAGGCCTAATTAGTAGCTTGATTAaaatcactt                                                                                                                            >  1:241554/1‑244 (MQ=255)
                                                                              atatTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCCGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAgg                                                                                                                         >  2:262442/1‑251 (MQ=255)
                                                                                           tACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAat                                                                                                             >  1:400583/1‑250 (MQ=255)
                                                                                               tGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAAcc                                                                                                                                                                            >  1:309639/1‑183 (MQ=255)
                                                                                               tGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAAcc                                                                                                                                                                            <  2:309639/183‑1 (MQ=255)
                                                                                                             gcTTGTTCCAGGAAATACGCTAAGTTTGACCCTGGTATGTTTTATCATCCGAGATACGAACATTGGAATGCGTCACATCGTAAGAAAATGCACCAACTGAAGTCTTTAGACGAAGACCTAACAAACCAGAGGTATAGTTATTATCGGTTATCTGAATACCGGAATAACCCGTCAGATAGTTATTCAGGCCGGAGAAGTAGCTAGCTTGAAATAAATTAGGATCATCCTGAATATCCTATTTAAAGAcccg                                                                                           >  2:436860/2‑248 (MQ=255)
                                                                                                                                aTAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGTATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACAGCCGCTAATATCCCAAc                                                                       >  2:26574/1‑251 (MQ=255)
                                                                                                                                  aaCTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGt                                                                     >  1:306064/1‑251 (MQ=255)
                                                                                                                                    ttctaattttGGTATTTTTTATCATCGGGAATACAAACATTGGAATTAGTCACATAAAAAAAAAATGCACCAACGGAAGTATTCAGACCAAGACCAAACAACCCAGCGGTATATTTATTATCGGTTATCTGAATACCGGTATAACCTGTAAGATAGTTATTCAGGCCTTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATTTTTTAAGACCGGCCCGCCGCTAATATCCCAACGTc                                                                    <  2:273992/241‑1 (MQ=255)
                                                                                                                                    cTTTGACCCTGGTATGTTTTCTCATCCGGGATAAGACCATTGGAATGAGTCACATCGAAAGAACATGCACCAACTGAAGCATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTATCTTGCTTGAAATAAATTAGGCTCATCATGACTATCATCTTTTAAGACCTGACAGCCGCTAATATCCCAACGTcc                                                                   >  2:33976/1‑251 (MQ=255)
                                                                                                                                           tcTGGTATGTTTTATCATTCAGGATCCGTACATTGGAATGAGTCTCATCGAAAGCACATGCACCAACTGACATATTTAGACAACGACCTAACAAACCAGCGGTATAGTTATAATCGGTTATCTGAATACCGGTATAACACGTAAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTAATAATGAATATCATCTTTTAAGACCTGACCGCCGCATATATCCCAACGTCCAACgacc                                                            >  2:232460/2‑248 (MQ=255)
                                                                                                                                             tGGTATGTTTTATCATTCGGGATACGAACATTGGAATTAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTc                                                                    >  1:273992/1‑241 (MQ=255)
                                                                                                                                                tATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCg                                                       >  1:229255/1‑251 (MQ=255)
                                                                                                                                                       tATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAAAAtt                                                >  2:397783/1‑251 (MQ=255)
                                                                                                                                                       tATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGCCCAACGCCAGGGCGAAACAtt                                                >  2:427372/1‑251 (MQ=255)
                                                                                                                                                                 gATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACgg                                      >  1:373801/1‑251 (MQ=255)
                                                                                                                                                                                               agaaaATGCACCAACTGAAGTATTTAGACCAAGACCCAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAgg                            >  1:460720/1‑231 (MQ=255)
                                                                                                                                                                                                      gCACAAACTGAAGTTTTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGTTTATCTGAATACCGGTATAACCCGTCAGATATTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTAATCCTAAATATCATCTTTTAAGCCCTTCCCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCAt  <  1:422787/250‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
GGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCAT  >  NC_000913/153967‑154414

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