Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 561,551 C→A Q619K (CAA→AAA)  sfmD → putative fimbrial usher protein SfmD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,5510CA95.7% 72.8 / ‑2.4 23Q619K (CAA→AAA) sfmDputative fimbrial usher protein SfmD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (10/12);  total (10/13)
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.39e-01

ATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGG  >  NC_000913/561416‑561684
                                                                                                                                       |                                                                                                                                     
aTCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAaaaa                                                                                                                                   >  1:436802/1‑140 (MQ=255)
  ccccGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATa                                                                                                                                 >  1:348295/1‑140 (MQ=255)
          acgaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAgctggctg                                                                                                                         <  1:307480/140‑1 (MQ=255)
          acgaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAgctggctg                                                                                                                         <  2:241466/140‑1 (MQ=255)
          acgaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAgctggctg                                                                                                                         <  2:222156/140‑1 (MQ=255)
           cgaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGctggctgg                                                                                                                        >  2:50223/1‑140 (MQ=255)
           cgaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGctggctgg                                                                                                                        <  2:47633/140‑1 (MQ=255)
            gaacaTATTGTCGGACTTAATGTTTCAGTACCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTGCTTTGACGCCAACCGTAACAGCAACGGGAAAAATAGCtggctggc                                                                                                                       <  2:356795/140‑1 (MQ=255)
                          aCTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGa                                                                                                         <  2:31579/140‑1 (MQ=255)
                             tAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCt                                                                                                      <  2:436802/140‑1 (MQ=255)
                                                  tGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGa                                                                                 >  2:196205/1‑140 (MQ=255)
                                                  tGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGa                                                                                 >  2:251551/1‑140 (MQ=255)
                                                      ttGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTa                                                                             <  2:94577/140‑1 (MQ=255)
                                                      ttGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTa                                                                             >  1:176032/1‑140 (MQ=255)
                                                           cAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACg                                                                        <  2:198267/140‑1 (MQ=255)
                                                           cAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCACTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACg                                                                        <  2:310462/140‑1 (MQ=255)
                                                                            gggAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATAcc                                                       >  1:398155/1‑140 (MQ=255)
                                                                              gAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTc                                                     <  1:350859/140‑1 (MQ=255)
                                                                                                                aCGCCAACCGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAAc                   >  2:130226/1‑140 (MQ=255)
                                                                                                                        cGTAACAGCAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGc           <  2:348295/140‑1 (MQ=255)
                                                                                                                           aaCAGCAACGGGAAAAATAGTTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGc        <  2:252386/140‑1 (MQ=255)
                                                                                                                                cAACGGGAAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACg   >  2:20119/1‑140 (MQ=255)
                                                                                                                                 aaCGGGAAAAATAGCTGGCTGGCAGGAGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACgg  >  2:257778/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                     
ATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGG  >  NC_000913/561416‑561684

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