Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,503,965 T→C D67D (GAT→GAC hicB → type II toxin‑antitoxin system antitoxin HicB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,503,9650TC93.8% 42.2 / ‑3.4 16D67D (GAT→GAChicBtype II toxin‑antitoxin system antitoxin HicB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (7/8);  total (7/9)
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.12e-01

GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAG  >  NZ_CP009273/1503879‑1504049
                                                                                      |                                                                                    
gCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                   <  2:185902/90‑1 (MQ=255)
          cTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAg                                                                         <  1:560832/90‑1 (MQ=255)
                 gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                   >  1:367542/1‑73 (MQ=255)
                 gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                   <  2:367542/73‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:30448/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  2:426977/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGACAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:408998/90‑1 (MQ=255)
                                ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTaa                                                   <  1:176791/90‑1 (MQ=255)
                                ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTaa                                                   >  1:529348/1‑90 (MQ=255)
                                        aaGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGc                                           <  2:674679/90‑1 (MQ=255)
                                                    ttATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCtt                               >  1:479111/1‑90 (MQ=255)
                                                          cTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTAc                         <  2:546411/90‑1 (MQ=255)
                                                            ttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAg                       >  1:429688/1‑90 (MQ=255)
                                                                               gTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGAATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAag    >  2:498458/1‑90 (MQ=255)
                                                                                 cACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAagag  >  1:195664/1‑90 (MQ=255)
                                                                                 cACGACCACCTTTTTGCAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAATGGCTTTTTTACAGTCATAAAACACTCAGCAagag  >  1:414710/1‑90 (MQ=255)
                                                                                      |                                                                                    
GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAG  >  NZ_CP009273/1503879‑1504049

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