Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 578,971 (A)6→5 intergenic (+380/‑166) BW25113_RS26085 → / → appY tail fiber assembly protein/DNA‑binding transcriptional activator AppY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273578,9660A.87.5% 18.2 / ‑0.5 8intergenic (+375/‑171)BW25113_RS26085/appYtail fiber assembly protein/DNA‑binding transcriptional activator AppY
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (4/3);  total (4/4)
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.65e-01

GCAAGAGGTTTCAGGTGCGTTGTAGTGAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAAACAACATACATATAATAATTTAATCTTAAATGAAATTTATTAAAATTTGCAAACTATAATTTTGTGTATAAAAATAT  >  NZ_CP009273/578877‑579047
                                                                                         |                                                                                 
gCAAGAGGTTTCAGGTGCGTTGTAGTGAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATa                                                                                   >  1:332677/1‑90 (MQ=255)
     agGTTTCAGGTGCGTTGTAGTGAGTTTATGTTAATGAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATaaaaaa                                                                              >  1:172894/1‑90 (MQ=255)
            aGGTGCGTTGTAGTGAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAACAacatac                                                                      <  1:216234/90‑1 (MQ=255)
                  gTTGTAGTGAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAACAACATACATataa                                                                >  1:99570/1‑90 (MQ=255)
                          gAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAacaaca                                                                         <  1:221681/73‑1 (MQ=255)
                          gAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAacaaca                                                                         >  2:221681/1‑73 (MQ=255)
                          gAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAACAACATACATATAATAATTTaa                                                        >  1:71742/1‑90 (MQ=255)
                                                       cGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAACAACATACATATAATAATTTAATCTTAAATGAAATTTATTAAAATTTGCaa                           <  2:71742/90‑1 (MQ=255)
                                                        gTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAAACAACATACATATAATAATTTAATCTTAAATGAAATTTATTAAAATTTGCaa                           <  2:172894/90‑1 (MQ=255)
                                                                                aGTTAGAATAAAAACAACATACATATAATAATTTAATCTTAAATGAAATTTATTAAAATTTGCAAACTATAATTTTGTGTATAAAAatat  >  2:171416/1‑90 (MQ=255)
                                                                                         |                                                                                 
GCAAGAGGTTTCAGGTGCGTTGTAGTGAGTTTATGTTAATAAAAAGCATAGTAAGCGTTGAAAAATGTAACTTTGAAATAAGTTAGAATAAAAAACAACATACATATAATAATTTAATCTTAAATGAAATTTATTAAAATTTGCAAACTATAATTTTGTGTATAAAAATAT  >  NZ_CP009273/578877‑579047

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