Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,353,081 G→A *58* (TAG→TAA ghoT → toxin of GhoTS toxin‑antitoxin pair, membrane‑lytic protein, stimulator of persister cell formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,353,0810GA97.7% 135.5 / ‑3.5 43coding (174/174 nt)ghoTtoxin of GhoTS toxin‑antitoxin pair, membrane‑lytic protein, stimulator of persister cell formation
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/20);  minor base . (0/1);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

ATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATT  >  NC_000913/4352958‑4353200
                                                                                                                           |                                                                                                                       
atCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                      <  1:106481/127‑1 (MQ=255)
atCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                      >  2:106481/1‑127 (MQ=255)
atCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGt                                                                                                          >  1:849050/1‑139 (MQ=255)
   cTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGa                                                                                                       <  1:527020/139‑1 (MQ=255)
   cTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGa                                                                                                       <  2:569377/139‑1 (MQ=255)
     ttGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGc                                                                                                                     >  2:21017/1‑123 (MQ=255)
     ttGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGc                                                                                                                     <  1:21017/123‑1 (MQ=255)
     ttGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGAcc                                                                                                     >  2:599591/1‑139 (MQ=255)
           ttATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTgg                                                                                               <  2:647160/139‑1 (MQ=255)
                  gTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGt                                                                                        >  2:455350/1‑139 (MQ=255)
                      aTCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGaaa                                                                                    <  2:849050/139‑1 (MQ=255)
                       tctcACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                      >  2:630576/1‑104 (MQ=255)
                       tctcACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                      <  1:630576/104‑1 (MQ=255)
                       tctcACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAAc                                                                                   >  2:638328/1‑139 (MQ=255)
                        ctcACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACg                                                                                  >  2:300948/1‑139 (MQ=255)
                            caTGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTA‑‑AGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                            <  2:997625/139‑1 (MQ=255)
                                agaAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtct                                                                          <  1:455350/139‑1 (MQ=255)
                                 gaAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtctc                                                                         <  2:1009826/139‑1 (MQ=255)
                                  aaaaCACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAg                                                                                                                <  1:82761/99‑1 (MQ=255)
                                  aaaaCACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAg                                                                                                                >  2:82761/1‑99 (MQ=255)
                                     acacATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTctcct                                                                     <  1:599591/139‑1 (MQ=255)
                                                  tACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                            <  2:521177/119‑1 (MQ=255)
                                                  tACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                            >  1:521177/1‑119 (MQ=255)
                                                   aCTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTg                                                       >  2:479277/1‑139 (MQ=255)
                                                       aGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGa                                                   >  2:873972/1‑139 (MQ=255)
                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGc                                                                                                  <  2:19041/90‑1 (MQ=255)
                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGc                                                                                                  >  1:19041/1‑90 (MQ=255)
                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAAc                                                 >  2:639844/1‑139 (MQ=255)
                                                          gCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTAACGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACg                                                <  1:873972/139‑1 (MQ=255)
                                                             ttCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATg                                             >  2:86044/1‑139 (MQ=255)
                                                              tCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCg                                                                                         <  2:387880/94‑1 (MQ=255)
                                                              tCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCg                                                                                         >  1:387880/1‑94 (MQ=255)
                                                                        aaTAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGt                                                                                              >  1:580376/1‑79 (MQ=255)
                                                                        aaTAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGt                                                                                              <  2:580376/79‑1 (MQ=255)
                                                                        aaTAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGt                                  <  1:86044/139‑1 (MQ=255)
                                                                        aaTAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGt                                  <  1:300948/139‑1 (MQ=255)
                                                                          tAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTgg                                                                                               >  2:547648/1‑76 (MQ=255)
                                                                          tAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTgg                                                                                               <  1:547648/76‑1 (MQ=255)
                                                                                           cTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTgg                                                                                                        <  1:182489/50‑1 (MQ=255)
                                                                                           cTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTgg                                                                                                        >  2:182489/1‑50 (MQ=255)
                                                                                                  tgGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGt        >  1:855271/1‑139 (MQ=255)
                                                                                                    gCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGa      >  1:384470/1‑139 (MQ=255)
                                                                                                        tACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGAc                                                                                                      >  1:96091/1‑39 (MQ=255)
                                                                                                        tACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGAc                                                                                                      <  2:96091/39‑1 (MQ=255)
                                                                                                        tACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatt  >  2:1105237/1‑139 (MQ=255)
                                                                                                                     ctTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGaa     <  1:866380/123‑1 (MQ=255)
                                                                                                                     ctTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGaa     >  2:866380/1‑123 (MQ=255)
                                                                                                                           |                                                                                                                       
ATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATT  >  NC_000913/4352958‑4353200

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