Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,425,126 G→C W607C (TGG→TGC stfR → prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,425,1260GC100.0% 58.0 / NA 20W607C (TGG→TGCstfRprophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (10/10);  total (10/10)

CGATCCGATGAATAACTGCGAGTTTAACGGATTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGGCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTG  >  NZ_CP009273/1425039‑1425208
                                                                                       |                                                                                  
tgctCCGCTGAATAACTGCGCTTCTACCGGATTTTTTATGCCTGGTGGCTCGAATGCCAGGGTTCCTTATTCTTACGGCATCTTCTGCCa                                                                                  <  1:364089/87‑1 (MQ=255)
             aaCTGCGAGTTTAACGGATTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAAt                                                                     >  1:389710/1‑90 (MQ=255)
                 gCGAGTTTAACGGATTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAAc                                                                 >  2:349287/1‑90 (MQ=255)
                            ggATTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCAc                                                      <  2:220334/90‑1 (MQ=255)
                              aTTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTc                                                    <  1:330536/90‑1 (MQ=255)
                              aTTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTc                                                    <  2:545541/90‑1 (MQ=255)
                                       gCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAatgat                                           >  1:166769/1‑90 (MQ=255)
                                        ccTGGTGGCTGGACTGACAGGGGGCGTTATGTTTATGGCATGTTCTGCCAATATCAAACAATGACCGAGCCATTCACTCAATAAtgatga                                         <  1:283275/90‑1 (MQ=255)
                                          tggtggCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAg                                        >  2:169461/1‑90 (MQ=255)
                                            gtggCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAAtgatga                                         >  2:570064/1‑87 (MQ=255)
                                            gtggCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAAtgatga                                         <  1:570064/87‑1 (MQ=255)
                                            gtggCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGta                                      >  2:289490/1‑90 (MQ=255)
                                                cTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGtaataa                                  >  1:348056/1‑90 (MQ=255)
                                                        aCAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCgatg                          <  1:87019/90‑1 (MQ=255)
                                                        aCAGGGGGCGTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCgatg                          <  1:349287/90‑1 (MQ=255)
                                                                 gTTATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGct                 >  2:581282/1‑90 (MQ=255)
                                                                  ttATGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGctc                <  1:540016/90‑1 (MQ=255)
                                                                    aTGCTTATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCtg              <  1:374624/90‑1 (MQ=255)
                                                                         tATGGCATGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTc         >  2:104286/1‑90 (MQ=255)
                                                                                tGTTCTGCCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTg  >  2:47131/1‑90 (MQ=255)
                                                                                       |                                                                                  
CGATCCGATGAATAACTGCGAGTTTAACGGATTTGTTATGCCTGGTGGCTGGACTGACAGGGGGCGTTATGCTTATGGCATGTTCTGGCAATATCAAAACAATGAACGAGCCATTCACTCAATAATGATGAGTAATAAGGGCGATGATTTGCGCTCTGTGTTCTATGTTG  >  NZ_CP009273/1425039‑1425208

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