Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,061,342 T→A intergenic (‑576/‑774) BW25113_RS10415 ← / → BW25113_RS24005 IS5‑like element IS5 family transposase/TonB‑dependent receptor plug domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,061,3420TA83.9% 82.3 / 10.2 31intergenic (‑576/‑774)BW25113_RS10415/BW25113_RS24005IS5‑like element IS5 family transposase/TonB‑dependent receptor plug domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base A (10/16);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-02

AAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  >  NZ_CP009273/2061258‑2061426
                                                                                    |                                                                                    
aaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTa                                                                                 >  1:567580/1‑90 (MQ=255)
  gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAgt                                                                               <  1:739646/90‑1 (MQ=255)
             tctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGaa                                                                    <  2:398146/90‑1 (MQ=255)
             tctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGaa                                                                    <  1:153145/90‑1 (MQ=255)
                            tctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                     <  1:542165/90‑1 (MQ=255)
                                 gACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGt                                                <  1:153974/90‑1 (MQ=255)
                                 gACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGt                                                <  2:35052/90‑1 (MQ=255)
                                   cacaGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAc                                              <  2:81311/90‑1 (MQ=255)
                                       gCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACt                                          >  1:557267/1‑90 (MQ=255)
                                         tttCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCg                                        >  2:332863/1‑90 (MQ=255)
                                            cTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtca                                     >  1:605619/1‑90 (MQ=255)
                                               cTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAc                                           >  2:530524/1‑81 (MQ=255)
                                               cTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAc                                           <  1:530524/81‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  1:677326/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:173292/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:35089/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:539150/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:557267/90‑1 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcatc                                   >  1:604665/1‑85 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcatc                                   <  2:604665/85‑1 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTa                              >  2:242832/1‑90 (MQ=255)
                                                         aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtata                        >  2:607842/1‑90 (MQ=255)
                                                         aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtata                        >  2:655424/1‑90 (MQ=255)
                                                           ggCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata                      >  1:2095/1‑90 (MQ=255)
                                                           ggCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata                      >  2:759523/1‑90 (MQ=255)
                                                                    aCGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGc             <  2:459384/90‑1 (MQ=255)
                                                                       aaGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAt          <  2:2653/90‑1 (MQ=255)
                                                                            gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg     >  2:767286/1‑90 (MQ=255)
                                                                            gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg     <  1:463919/90‑1 (MQ=255)
                                                                            gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg     <  2:836742/90‑1 (MQ=255)
                                                                               cAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTgtgt  <  2:673315/90‑1 (MQ=255)
                                                                                    |                                                                                    
AAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  >  NZ_CP009273/2061258‑2061426

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