Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,778,934 (T)6→5 intergenic (‑564/+187) BW25113_RS25705 ← / ← BW25113_RS13825 DUF5507 domain‑containing protein/tRNA‑Ile

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,778,9290T.94.7% 66.6 / ‑0.8 19intergenic (‑559/+192)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (7/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAA  >  NZ_CP009273/2778841‑2778996
                                                                                        |                                                                   
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                    >  1:249977/1‑90 (MQ=255)
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                    >  2:152495/1‑90 (MQ=255)
            aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                       >  1:284811/1‑90 (MQ=255)
            aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                       <  1:32265/90‑1 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                           >  1:332539/1‑90 (MQ=255)
                           gCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTg                                        <  1:431020/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                      <  2:284811/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                      <  1:33338/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                      <  1:418039/90‑1 (MQ=255)
                              ttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATg                                     >  2:373364/1‑90 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                      >  2:31111/1‑88 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                      <  1:31111/88‑1 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATgg                                    >  1:131775/1‑90 (MQ=255)
                                       aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                            <  1:117279/90‑1 (MQ=255)
                                       aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                            <  2:131775/90‑1 (MQ=255)
                                             tAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAATTACATTGCa                      >  2:198257/1‑90 (MQ=255)
                                                         tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc          <  2:307923/90‑1 (MQ=255)
                                                          aaGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc          >  2:214487/1‑90 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                               >  2:99200/1‑61 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa                               <  1:99200/61‑1 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa  <  1:88849/90‑1 (MQ=255)
                                                                                        |                                                                   
AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAA  >  NZ_CP009273/2778841‑2778996

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