Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 560,244 A→G N183S (AAT→AGT)  sfmD → putative fimbrial usher protein SfmD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913560,2440AG100.0% 60.1 / NA 18N183S (AAT→AGT) sfmDputative fimbrial usher protein SfmD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/8);  total (10/8)

TTGAGTAGCGCTCACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGC  >  NC_000913/560177‑560320
                                                                   |                                                                            
ttgAGTAGCGCTCACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCAcc                                                                 <  2:111919/81‑1 (MQ=255)
  gAGTAGCGCTCACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCgg                                                               >  2:111940/1‑81 (MQ=255)
  gAGTAGCGCTCACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCgg                                                               >  1:111950/1‑81 (MQ=255)
            cACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAgg                                                      >  1:111942/1‑80 (MQ=255)
            cACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAgg                                                      >  1:111948/1‑80 (MQ=255)
                gTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATg                                                 <  1:111920/81‑1 (MQ=255)
                    cATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTaa                                             <  1:111921/81‑1 (MQ=255)
                        ccGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGa                                          >  1:111944/1‑80 (MQ=255)
                         cGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATa                                        <  1:111922/81‑1 (MQ=255)
                           ccTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGc                                      >  1:111943/1‑81 (MQ=255)
                           ccTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGc                                      <  1:111923/81‑1 (MQ=255)
                            cTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGc                                      >  1:111951/1‑80 (MQ=255)
                                   gTGGGATGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTAtttttt                              >  2:111941/1‑81 (MQ=255)
                                        aTGAAGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGt                          >  1:111947/1‑80 (MQ=255)
                                           aaGGTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGCAATGGTAACGATAGCTATTTTTTTAGTGAGc                      >  2:111953/1‑81 (MQ=255)
                                              gTATTCCTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGTGAGCTCa                   <  2:111924/81‑1 (MQ=255)
                                                    cTGCTTTACTCCTGAGTTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGa             <  2:111925/81‑1 (MQ=255)
                                                               cTGAGTTATAATTTCACCGGTAACAGAGGTAGTGGTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGc  <  2:111926/81‑1 (MQ=255)
                                                                   |                                                                            
TTGAGTAGCGCTCACGGTTACATTCCGCCTGAAGAGTGGGATGAAGGTATTCCTGCTTTACTCCTGAATTATAATTTCACCGGTAACAGAGGTAATGGTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGC  >  NC_000913/560177‑560320

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