Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,921 T→C 7.5% Q75R (CAA→CGA)  higA ← antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,9210TC 7.5% 162.3 / 3.7 53Q75R (CAA→CGA) higAantitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTT  >  NC_000913/3233782‑3234062
                                                                                                                                           |                                                                                                                                             
gtgtTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgc                                                                                                                                      <  4:78067/149‑2 (MQ=255)
   ttCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt                                                                                                                                   >  1:233889/1‑149 (MQ=255)
    tCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAtt                                                                                                                                  <  4:468238/149‑1 (MQ=255)
     ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTa                                                                                                                                 >  3:461041/1‑149 (MQ=255)
          gTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                                                                                                                            >  2:318496/1‑149 (MQ=255)
          gTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgc                                                                                                                                      <  4:234439/139‑2 (MQ=255)
          gTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgc                                                                                                                                      >  3:234439/1‑138 (MQ=255)
           ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc                                                                                                                           >  1:426731/1‑149 (MQ=255)
                        cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcg                                                                                                              >  4:124511/1‑149 (MQ=255)
                        cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcg                                                                                                              >  2:305402/1‑149 (MQ=255)
                         ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                                                             >  4:70072/1‑149 (MQ=255)
                           gCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg                                                                                                           <  3:124511/149‑1 (MQ=255)
                                 aaCGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttctt                                                                                                     >  3:41764/1‑149 (MQ=255)
                                            ccATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGtt                                                                                          >  4:45672/1‑149 (MQ=255)
                                              aTAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAt                                                                                        >  2:132661/1‑149 (MQ=255)
                                              aTAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAt                                                                                        >  2:375054/1‑149 (MQ=255)
                                                aGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAttt                                                                                      <  1:453535/149‑1 (MQ=255)
                                                           cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACAcc                                                                           >  1:409094/1‑149 (MQ=255)
                                                                  cGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCa                                                                    <  3:253050/149‑1 (MQ=255)
                                                                    gCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAttt                                                                                      >  4:22144/1‑129 (MQ=255)
                                                                    gCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAttt                                                                                      <  3:22144/129‑1 (MQ=255)
                                                                      aGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCaa                                                                >  4:101578/1‑149 (MQ=255)
                                                                       gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAg                                                               <  2:82863/149‑1 (MQ=255)
                                                                       gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAg                                                               >  4:212627/1‑149 (MQ=255)
                                                                          cGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAgggg                                                            <  2:233889/149‑1 (MQ=255)
                                                                            gAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGtt                                                          <  3:70072/149‑1 (MQ=255)
                                                                               aGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTc                                                       >  4:55933/1‑149 (MQ=255)
                                                                               aGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTc                                                       <  1:318496/149‑1 (MQ=255)
                                                                                gggTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACCCACCAGATTCAGCAAGGGGTTTTCa                                                      <  4:26305/149‑1 (MQ=255)
                                                                                   ttAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGa                                                   >  2:381704/1‑149 (MQ=255)
                                                                                    tAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGAt                                                  <  2:426731/149‑1 (MQ=255)
                                                                                    tAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGAt                                                  <  1:155532/149‑1 (MQ=255)
                                                                                        ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGtt                                              <  4:284217/149‑1 (MQ=255)
                                                                                         cATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTg                                             >  4:380887/1‑149 (MQ=255)
                                                                                          atatTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGa                                            <  3:21145/149‑1 (MQ=255)
                                                                                            atTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGagc                                          <  4:169425/149‑1 (MQ=255)
                                                                                                 tCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGcagcag                                     >  4:87470/1‑149 (MQ=255)
                                                                                                    aTAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGagc                                          <  3:384872/141‑1 (MQ=255)
                                                                                                    aTAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGagc                                          >  4:384872/1‑141 (MQ=255)
                                                                                                     tAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAgatga                                 <  1:305402/149‑1 (MQ=255)
                                                                                                        gggTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCt                              <  3:212627/149‑1 (MQ=255)
                                                                                                         ggTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTa                             <  3:45672/149‑1 (MQ=255)
                                                                                                          gTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACTCCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTAc                            <  4:410929/149‑1 (MQ=255)
                                                                                                           tACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTAcc                           >  4:371812/1‑149 (MQ=255)
                                                                                                           tACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTAcc                           >  1:33317/1‑149 (MQ=255)
                                                                                                           tACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCAGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTAcc                           >  3:182763/1‑149 (MQ=255)
                                                                                                                aTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTc                      >  3:176570/1‑149 (MQ=255)
                                                                                                                  cACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCa                    <  3:55933/149‑1 (MQ=255)
                                                                                                                   aCGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAg                   >  3:216182/1‑149 (MQ=255)
                                                                                                                       cTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGcc               <  3:371812/149‑1 (MQ=255)
                                                                                                                       cTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGcc               <  3:314054/149‑1 (MQ=255)
                                                                                                                                 aGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTAtt     >  3:3992/1‑149 (MQ=255)
                                                                                                                                    cATGGCTCAAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGtt  >  4:320046/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTT  >  NC_000913/3233782‑3234062

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑