Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,446,677 T→C pseudogene (1804/2646 nt) yfcU ← fimbrial usher protein YfcU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,446,6770TC100.0% 61.8 / NA 19pseudogene (1804/2646 nt)yfcUfimbrial usher protein YfcU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/11);  total (8/11)

GTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCG  >  NZ_CP009273/2446596‑2446760
                                                                                 |                                                                                   
gtTTGTCACTGGTGCCGATGTTCAACTGATAGTTAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGc                                                                             >  2:107156/1‑90 (MQ=255)
     tCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAg                                                                        >  2:232647/1‑90 (MQ=255)
      cACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGt                                                                       >  1:93300/1‑90 (MQ=255)
         tGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGc                                                                    <  1:232647/90‑1 (MQ=255)
          ggTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGcc                                                                   >  2:42355/1‑90 (MQ=255)
                      cAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAAc                                                       <  2:69616/90‑1 (MQ=255)
                           gATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTAc                                                  <  1:34085/90‑1 (MQ=255)
                                  aGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTg                                                                              >  1:358253/1‑55 (MQ=255)
                                  aGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTg                                                                              <  2:358253/55‑1 (MQ=255)
                                   gTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGGACTGTTGTCg                                          >  1:67175/1‑90 (MQ=255)
                                      gcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGccc                                       <  2:73196/90‑1 (MQ=255)
                                      gcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGccc                                       <  1:107156/90‑1 (MQ=255)
                                      ccggcATCGCCACGGCTAAAATACCCGCCCTGGCCGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGccc                                       <  2:449278/86‑1 (MQ=255)
                                                                 aCCTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATg            <  2:338209/90‑1 (MQ=255)
                                                                   cTGGCTGCTGTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGcc          <  1:399023/90‑1 (MQ=255)
                                                                         gctgTCCGCCCCCCTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATc    <  2:236316/90‑1 (MQ=255)
                                                                          ctgTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACt                            <  2:110137/65‑1 (MQ=255)
                                                                          ctgTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACt                            >  1:110137/1‑65 (MQ=255)
                                                                           tgTCCGCCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCg  >  2:293114/1‑90 (MQ=255)
                                                                                 |                                                                                   
GTTTGTCACTGGTGCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCG  >  NZ_CP009273/2446596‑2446760

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