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,9650TC81.0% 36.7 / 6.2 21D67D (GAT→GAChicBtype II toxin‑antitoxin system antitoxin HicB
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (5/12);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAGT  >  NZ_CP009273/1503879‑1504050
                                                                                      |                                                                                     
gCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                    <  1:520788/90‑1 (MQ=255)
     aGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTAt                                                                               <  2:101372/90‑1 (MQ=255)
        tGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGa                                                                            <  1:268699/90‑1 (MQ=255)
           tttACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGt                                                                         <  1:97672/90‑1 (MQ=255)
                 gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                   >  1:546124/1‑90 (MQ=255)
                 gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                   >  1:295544/1‑90 (MQ=255)
                       aTTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGt                                                             <  1:48716/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                            <  2:809609/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                            <  2:36951/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                            <  1:310308/90‑1 (MQ=255)
                         ttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                           <  1:110073/90‑1 (MQ=255)
                                ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTaa                                                    >  1:250751/1‑90 (MQ=255)
                                ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTaa                                                    >  2:167102/1‑90 (MQ=255)
                                                gAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAAt                                    <  2:251347/90‑1 (MQ=255)
                                                       tCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCttttt                             >  2:264527/1‑90 (MQ=255)
                                                               ccTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCa                     <  1:264527/90‑1 (MQ=255)
                                                                  tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                  >  2:170188/1‑90 (MQ=255)
                                                                             tAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCa       <  1:225707/90‑1 (MQ=255)
                                                                               gTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCa                     <  2:12376/74‑1 (MQ=255)
                                                                               gTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCa                     >  1:12376/1‑74 (MQ=255)
                                                                                  aCGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAGt  <  1:227444/90‑1 (MQ=255)
                                                                                      |                                                                                     
GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAGT  >  NZ_CP009273/1503879‑1504050

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: