Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,524,293 A→G intergenic (+98/‑105) BW25113_RS07645 → / → BW25113_RS26145 RHS repeat protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,524,2930AG75.0% 44.1 / 10.9 24intergenic (+98/‑105)BW25113_RS07645/BW25113_RS26145RHS repeat protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (8/10);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

AGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACT  >  NZ_CP009273/1524164‑1524422
                                                                                                                                 |                                                                                                                                 
agcaaGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCa                                                                                                                           >  1:748559/1‑138 (MQ=255)
  caaGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAag                                                                                                                        <  1:761271/139‑1 (MQ=255)
      ggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCg                                                                                                                    <  2:759100/139‑1 (MQ=255)
      ggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCg                                                                                                                    <  1:110436/139‑1 (MQ=255)
                             taGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGt                                                                                                                   <  2:7075/117‑1 (MQ=255)
                             taGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGt                                                                                                                   >  1:7075/1‑117 (MQ=255)
                             taGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaa                                                                                              >  2:75194/1‑138 (MQ=255)
                                 ggCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAaga                                                                                                                       >  2:775562/1‑109 (MQ=255)
                                 ggCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAaga                                                                                                                       <  1:775562/109‑1 (MQ=255)
                                          gTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATAttatttat                                                                                 >  1:811761/1‑138 (MQ=255)
                                          gTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATAttattta                                                                                  <  2:811761/137‑1 (MQ=255)
                                                  tttAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGaa                                                                                                >  1:638464/1‑115 (MQ=255)
                                                  tttAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGaa                                                                                                <  2:638464/115‑1 (MQ=255)
                                                                    aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                       >  2:564486/1‑138 (MQ=255)
                                                                                    gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                      >  1:696646/1‑139 (MQ=255)
                                                                                         tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                                  <  2:284817/138‑1 (MQ=255)
                                                                                              aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                             <  2:412044/138‑1 (MQ=255)
                                                                                                     aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGt                      >  1:655522/1‑138 (MQ=255)
                                                                                                                     ttttATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg      <  1:564486/138‑1 (MQ=255)
                                                                                                                     ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg      <  2:449679/138‑1 (MQ=255)
                                                                                                                       ttATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCGTTTTATAAGTATAATGGTTAATAATTACAGGCCGCa    <  2:219732/138‑1 (MQ=255)
                                                                                                                        tATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACt  >  2:308024/1‑139 (MQ=255)
                                                                                                                             aCATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACt  >  2:216480/1‑134 (MQ=255)
                                                                                                                             aCATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACt  <  1:216480/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                 
AGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACT  >  NZ_CP009273/1524164‑1524422

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