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.94.7% 138.7 / 3.9 38intergenic (‑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/2);  new base . (8/28);  total (8/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

TTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTC  >  NZ_CP009273/2283783‑2283955
                                                                                        |                                                                                    
ttttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtt                                                                                     <  2:1484874/90‑1 (MQ=255)
 tttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAGCGAACGACAAGGTTGTAGtttc                                                                                   <  2:1504988/90‑2 (MQ=255)
  ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                                   >  2:537199/1‑89 (MQ=255)
  ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                                   <  1:1520300/90‑2 (MQ=255)
  ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                                   <  2:1096663/90‑2 (MQ=255)
       gttgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTAc                                                                             <  1:763357/90‑1 (MQ=255)
         tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACAt                                                                           >  1:33138/1‑90 (MQ=255)
              ttCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCca                                                                      <  1:661962/90‑1 (MQ=255)
              ttCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCca                                                                      <  2:247722/90‑1 (MQ=255)
                        gaAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAAt                                                            >  2:766572/1‑90 (MQ=255)
                        gaAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAAt                                                            <  1:755340/90‑1 (MQ=255)
                              gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtt                                                                                     >  1:387876/1‑60 (MQ=255)
                              gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtt                                                                                     <  2:387876/60‑1 (MQ=255)
                              gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTa                                                      <  1:755784/90‑1 (MQ=255)
                                    aaaTTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTc                                                >  1:174006/1‑90 (MQ=255)
                                       ttCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCAcg                                             <  1:381789/90‑1 (MQ=255)
                                         ccGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCAcgcg                                           <  1:818052/90‑1 (MQ=255)
                                          cGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCAcgcg                                           <  1:1356724/90‑1 (MQ=255)
                                          cGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACgcgc                                          <  1:631098/90‑1 (MQ=255)
                                              tttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTgg                                      >  2:862449/1‑90 (MQ=255)
                                               ttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGt                                     <  2:1261936/90‑1 (MQ=255)
                                               ttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGt                                     <  2:165628/90‑1 (MQ=255)
                                                tataATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTa                                    <  2:655304/90‑1 (MQ=255)
                                                   aaTGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCg                                 <  2:170766/90‑1 (MQ=255)
                                                     tGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGt                                <  2:528458/89‑1 (MQ=255)
                                                     tGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGt                                >  1:528458/1‑89 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTAcc                            <  2:598579/90‑1 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa                           <  1:914679/90‑1 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa                           <  1:1539507/90‑1 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa                           <  2:154034/90‑1 (MQ=255)
                                                          atgatgCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCaa                          <  2:174006/90‑1 (MQ=255)
                                                              tgCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgct                      >  2:130426/1‑90 (MQ=255)
                                                                ccAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAccgctac                    <  2:1201268/90‑1 (MQ=255)
                                                                ccAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgctac                    <  2:1271035/90‑1 (MQ=255)
                                                                  aaaGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgctacgc                  <  2:648419/90‑1 (MQ=255)
                                                                    aGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTc                <  1:1400911/90‑1 (MQ=255)
                                                                      cGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCaa              <  2:1581297/90‑1 (MQ=255)
                                                                            aCAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAAGATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATaa        <  1:245379/90‑1 (MQ=255)
                                                                            aCAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATaa        <  1:521192/90‑1 (MQ=255)
                                                                               aGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGa     >  1:894418/1‑90 (MQ=255)
                                                                                 gTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGAtt   <  1:581493/90‑1 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  <  2:1231874/90‑1 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  <  2:41921/90‑1 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  >  1:1260552/1‑90 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  <  1:1256853/90‑1 (MQ=255)
                                                                                        |                                                                                    
TTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTC  >  NZ_CP009273/2283783‑2283955

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