Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,283,874 Δ1 bp intergenic (‑350/‑105) BW25113_RS11435 ← / → narP IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,283,8710T.80.0% 27.2 / 10.6 15intergenic (‑347/‑108)BW25113_RS11435/narPIS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base . (4/8);  total (4/11)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

TTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAA  >  NZ_CP009273/2283785‑2283943
                                                                                      |                                                                        
ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                       <  1:318116/90‑2 (MQ=255)
       tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACAt                                                               <  1:172201/90‑1 (MQ=255)
       tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACAt                                                               <  1:750193/90‑1 (MQ=255)
             tCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCca                                                          <  1:732238/90‑1 (MQ=255)
             tCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCcat                                                         >  1:255719/1‑90 (MQ=255)
                       aaGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAAt                                                <  1:161441/90‑1 (MQ=255)
                           tGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGt                                           >  1:19644/1‑90 (MQ=255)
                           tGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGt                                           <  1:344244/90‑1 (MQ=255)
                            gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTa                                          >  2:651332/1‑90 (MQ=255)
                                        cGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACgcgc                              <  1:318355/90‑1 (MQ=255)
                                         gTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGcc                             >  1:497931/1‑90 (MQ=255)
                                                 aaTGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCg                     <  2:132547/90‑1 (MQ=255)
                                                       aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa               <  2:497931/90‑1 (MQ=255)
                                                              ccAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgctac        <  1:495140/90‑1 (MQ=255)
                                                                  aGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTc    <  1:690526/90‑1 (MQ=255)
                                                                     gaccGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCaa  <  2:484283/87‑1 (MQ=255)
                                                                                      |                                                                        
TTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAA  >  NZ_CP009273/2283785‑2283943

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