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,9650TC100.0% 88.3 / NA 28D67D (GAT→GAChicBtype II toxin‑antitoxin system antitoxin HicB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/13);  total (15/13)

GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAG  >  NZ_CP009273/1503879‑1504047
                                                                                      |                                                                                  
gCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                 <  1:863/90‑1 (MQ=255)
    aaGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTa                                                                             <  1:470666/90‑1 (MQ=255)
          cTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAg                                                                       >  1:174139/1‑90 (MQ=255)
           tttACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGt                                                                      <  1:563918/90‑1 (MQ=255)
              aCTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTAcc                                                                   >  1:189327/1‑90 (MQ=255)
                      cATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCg                                                           >  2:68606/1‑90 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                         <  2:66478/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                         <  2:189327/90‑1 (MQ=255)
                        tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                         <  1:69162/90‑1 (MQ=255)
                         ttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                        <  1:78328/90‑1 (MQ=255)
                             tttttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCctc                                                    >  1:449868/1‑90 (MQ=255)
                              ttttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCctct                                                   <  2:292889/90‑1 (MQ=255)
                                                 aGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATg                                >  1:399095/1‑90 (MQ=255)
                                                  gCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGc                               <  1:573818/90‑1 (MQ=255)
                                                    ttATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCtt                             >  1:330302/1‑90 (MQ=255)
                                                    ttATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCtt                             >  2:189927/1‑90 (MQ=255)
                                                           tttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCCGTTAAATGCTTTTTTACa                      <  1:275846/90‑1 (MQ=255)
                                                               ccTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCa                  <  1:40442/90‑1 (MQ=255)
                                                               ccTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCa                  <  1:162051/90‑1 (MQ=255)
                                                                  tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa               >  1:580985/1‑90 (MQ=255)
                                                                  tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa               >  2:269360/1‑90 (MQ=255)
                                                                  tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa               >  2:425138/1‑90 (MQ=255)
                                                                  tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa               >  2:552479/1‑90 (MQ=255)
                                                                  tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa               >  1:207995/1‑90 (MQ=255)
                                                                   cGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaaa              >  1:116453/1‑90 (MQ=255)
                                                                   cGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaaa              >  2:265353/1‑90 (MQ=255)
                                                                             tAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCa    <  1:549127/90‑1 (MQ=255)
                                                                               gTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAag  >  1:318190/1‑90 (MQ=255)
                                                                                      |                                                                                  
GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAG  >  NZ_CP009273/1503879‑1504047

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