Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,077,067 T→C pseudogene (256/830 nt) BW25113_RS05295 → FTR1 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,077,0670TC100.0% 71.6 / NA 22pseudogene (256/830 nt)BW25113_RS05295FTR1 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/13);  total (9/13)

ACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAG  >  NZ_CP009273/1076980‑1077140
                                                                                       |                                                                         
aCGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACgg                                                                         >  1:855411/1‑90 (MQ=255)
         ccGGCGAATTTCCGCAAAAATAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGGGATCCTTACCTGGACGGTTTTCTGGa                                                                >  1:265752/1‑90 (MQ=255)
                aTTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         <  2:1007677/90‑1 (MQ=255)
                   tCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         <  2:546947/87‑1 (MQ=255)
                   tCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         >  1:546947/1‑87 (MQ=255)
                          aaaGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCgcgc                                               <  1:1019801/90‑1 (MQ=255)
                          aaaGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCgcgc                                               <  2:974670/90‑1 (MQ=255)
                            aGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCaa                                             <  1:764303/90‑1 (MQ=255)
                                     aCTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCACGTCAAAGTg                                   <  1:228249/90‑1 (MQ=255)
                                       tGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         >  1:64014/1‑67 (MQ=255)
                                       tGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         <  2:64014/67‑1 (MQ=255)
                                        gTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCa                                 <  1:737880/90‑1 (MQ=255)
                                                  aTCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         >  1:756548/1‑56 (MQ=255)
                                                  aTCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         >  1:600287/1‑56 (MQ=255)
                                                  aTCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         <  2:600287/56‑1 (MQ=255)
                                                  aTCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCaa                                                         <  2:756548/56‑1 (MQ=255)
                                                  aTCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGACGGGTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAg                       >  2:775926/1‑90 (MQ=255)
                                                  aTCGTGGCGGTGATCGCCGTGGGGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAg                       >  1:1018321/1‑90 (MQ=255)
                                                        gCGGTGATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTc                 <  2:1018321/90‑1 (MQ=255)
                                                          ggTGATCGCCGTGGGGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGa               >  1:731097/1‑90 (MQ=255)
                                                             gATCGCCGTGGTGATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAg            <  1:660466/90‑1 (MQ=255)
                                                                       gtgATCCTTACCTGGACGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAg  <  1:1068289/90‑1 (MQ=255)
                                                                                       |                                                                         
ACGAAACCACCGGCGAATTTCCGCAAAAAGAACAGGAACTGTTTGAAGGTATCGTGGCGGTGATCGCCGTGGTGATCCTTACCTGGATGGTTTTCTGGATGCGCAAAGTGTCGCGCAACGTCAAAGTGCAACTGGAACAGGCAGTCGATAGCGCATTGCAG  >  NZ_CP009273/1076980‑1077140

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