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

GGCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAG  >  NZ_CP009273/1503878‑1504049
                                                                                       |                                                                                    
ggCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCa                                                                                    <  1:2225635/90‑1 (MQ=255)
         tGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGa                                                                           <  1:1571173/90‑1 (MQ=255)
          gCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGaa                                                                          >  1:2622381/1‑90 (MQ=255)
           cTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAg                                                                         <  2:1410059/90‑1 (MQ=255)
           cTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAg                                                                         <  1:91179/90‑1 (MQ=255)
           cTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAg                                                                         <  2:1746919/90‑1 (MQ=255)
             ttACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTa                                                                       <  2:2806182/90‑1 (MQ=255)
                  gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                  >  2:1997639/1‑90 (MQ=255)
                  gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                  >  1:2125265/1‑90 (MQ=255)
                  gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                  >  2:908024/1‑90 (MQ=255)
                    ccGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGa                                                                >  1:1576619/1‑90 (MQ=255)
                    ccGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGa                                                                >  1:29206/1‑90 (MQ=255)
                    ccGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGa                                                                >  1:1847825/1‑90 (MQ=255)
                       cATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCg                                                             >  1:1755520/1‑90 (MQ=255)
                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  2:1708014/90‑1 (MQ=255)
                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:377916/90‑1 (MQ=255)
                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:2894698/90‑1 (MQ=255)
                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:1323851/90‑1 (MQ=255)
                          ttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                          <  1:126746/90‑1 (MQ=255)
                            gATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGcc                                                        <  1:2215945/90‑1 (MQ=255)
                            gATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGcc                                                        <  2:884992/90‑1 (MQ=255)
                            gATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGcc                                                        <  1:59070/90‑1 (MQ=255)
                                 ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                          >  2:1414477/1‑83 (MQ=255)
                                 ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                          <  1:1414477/83‑1 (MQ=255)
                                         aaGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGc                                           <  2:1311475/90‑1 (MQ=255)
                                                  aGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATg                                  >  2:248443/1‑90 (MQ=255)
                                                   gCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGc                                 <  1:908024/90‑1 (MQ=255)
                                                      tATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCttt                              >  1:1820316/1‑90 (MQ=255)
                                                           cTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTAc                         <  2:1685362/90‑1 (MQ=255)
                                                           cTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTAc                         <  1:560551/90‑1 (MQ=255)
                                                           cTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTAc                         <  2:29206/90‑1 (MQ=255)
                                                            tttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGcc                                                        <  2:2469998/58‑1 (MQ=255)
                                                            tttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGcc                                                        >  1:2469998/1‑58 (MQ=255)
                                                            tttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACa                        >  1:1258382/1‑90 (MQ=255)
                                                            tttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACa                        >  1:287295/1‑90 (MQ=255)
                                                             ttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                          <  1:1856999/55‑1 (MQ=255)
                                                             ttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                          >  2:1856999/1‑55 (MQ=255)
                                                              tACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGt                      <  1:172474/90‑1 (MQ=255)
                                                                  ttCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGa                  >  1:2180214/1‑90 (MQ=255)
                                                                  ttCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGa                  >  1:2640280/1‑90 (MQ=255)
                                                                  ttCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGa                  >  2:2357990/1‑90 (MQ=255)
                                                                  ttCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGa                  >  1:695352/1‑90 (MQ=255)
                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  1:3012243/1‑90 (MQ=255)
                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  1:2107251/1‑90 (MQ=255)
                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  2:316903/1‑90 (MQ=255)
                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  2:889802/1‑90 (MQ=255)
                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  1:1802105/1‑90 (MQ=255)
                                                                        aTTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCAc            <  1:1254500/90‑1 (MQ=255)
                                                                        aTTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCAc            <  2:526800/90‑1 (MQ=255)
                                                                         ttAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACt           >  1:372233/1‑90 (MQ=255)
                                                                         ttAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACt           <  2:572051/90‑1 (MQ=255)
                                                                              tAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCa      <  2:287295/90‑1 (MQ=255)
                                                                              tAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCa      <  1:2360447/90‑1 (MQ=255)
                                                                                gTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGAATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAag    >  1:1353101/1‑90 (MQ=255)
                                                                                  cACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAagag  >  2:2714931/1‑90 (MQ=255)
                                                                                  cACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAagag  >  2:2955920/1‑90 (MQ=255)
                                                                                       |                                                                                    
GGCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAG  >  NZ_CP009273/1503878‑1504049

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