Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,041,028 G→A R15C (CGT→TGT)  sdhE ← flavinator of succinate dehydrogenase; antitoxin of CptAB toxin‑antitoxin pair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,0280GA95.0% 187.9 / 0.8 60R15C (CGT→TGT) sdhEflavinator of succinate dehydrogenase; antitoxin of CptAB toxin‑antitoxin pair
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (25/32);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

CAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGC  >  NC_000913/3040933‑3041124
                                                                                               |                                                                                                
cagACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcagcgg                                                                                              <  2:1576124/100‑1 (MQ=255)
cagACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcagcgg                                                                                              <  2:275787/100‑1 (MQ=255)
cagACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcagcgg                                                                                              >  1:254911/1‑100 (MQ=255)
    cGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATg                                                                                          <  2:227581/100‑1 (MQ=255)
      aaTAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGcc                                                                                        >  1:20595/1‑100 (MQ=255)
       aTAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGccc                                                                                       <  2:767954/100‑1 (MQ=255)
        tAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCa                                                                                      <  2:1497862/100‑1 (MQ=255)
        tAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCa                                                                                      >  2:773005/1‑100 (MQ=255)
        tAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCa                                                                                      <  2:231608/100‑1 (MQ=255)
          aaGATGCGTTTTTTTTTTTTTCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaat                                                                                    >  2:1099907/1‑100 (MQ=255)
          aaGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaat                                                                                    <  2:189600/100‑1 (MQ=255)
           aGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatg                                                                                   <  1:1245590/100‑1 (MQ=255)
              tGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatgaat                                                                                >  1:67938/1‑100 (MQ=255)
               gCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCaatgaatg                                                                               <  1:1438954/100‑1 (MQ=255)
                  tttttCGTCATCGCTTAAGCTGTCGTACAAAAGTACGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCgg                                                                            <  2:540699/100‑1 (MQ=255)
                       cGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCttt                                                                       >  2:1425454/1‑100 (MQ=255)
                        gTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTttg                                                                      <  1:686798/100‑1 (MQ=255)
                        gTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttg                                                                      <  2:470101/100‑1 (MQ=255)
                          cATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTttgtt                                                                    >  2:797608/1‑100 (MQ=255)
                            tCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTtgttgt                                                                  <  1:1355643/100‑1 (MQ=255)
                             cGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTgttgtt                                                                 >  1:335896/1‑100 (MQ=255)
                             cGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTgttgtt                                                                 <  2:191714/100‑1 (MQ=255)
                             cGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTgttgtt                                                                 <  1:1099912/100‑1 (MQ=255)
                             cGCTTAAGCTGTCGTACTCATGCTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTgttgtt                                                                 >  2:606486/1‑100 (MQ=255)
                              gCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTa                                                                <  2:31808/100‑1 (MQ=255)
                              gCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTa                                                                <  2:185313/100‑1 (MQ=255)
                                    gCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTcc                                                          <  1:122187/100‑1 (MQ=255)
                                     cTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCa                                                         >  1:357479/1‑100 (MQ=255)
                                      tGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCAt                                                        >  2:1293194/1‑100 (MQ=255)
                                         cGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCAtct                                                      >  1:65930/1‑99 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCTTGTTCGCTCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTcc                                                <  2:1160668/100‑1 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTcc                                                <  1:273697/100‑1 (MQ=255)
                                              tCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTcc                                                <  2:335896/100‑1 (MQ=255)
                                               cATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCt                                               >  1:761501/1‑100 (MQ=255)
                                                aTGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTg                                              >  1:608927/1‑100 (MQ=255)
                                                aTGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTg                                              >  1:1255088/1‑100 (MQ=255)
                                                  gTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGGCCATCTTCTTCCTGtc                                            >  1:866441/1‑100 (MQ=255)
                                                   ttCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGtct                                           <  1:912241/100‑1 (MQ=255)
                                                     cGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCa                                         >  1:799330/1‑100 (MQ=255)
                                                      gAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCAc                                        <  1:1127908/100‑1 (MQ=255)
                                                      gAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCAc                                        <  1:100841/100‑1 (MQ=255)
                                                      gAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCAc                                        >  1:162056/1‑100 (MQ=255)
                                                      gAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGCCTCAc                                        >  1:329934/1‑100 (MQ=255)
                                                       aaaaaCGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACg                                       <  2:1358304/100‑1 (MQ=255)
                                                            cGGCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAAt                                  >  2:1456964/1‑100 (MQ=255)
                                                              gCATGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATcc                                <  2:612542/100‑1 (MQ=255)
                                                                aTGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAg                              <  1:1172208/100‑1 (MQ=255)
                                                                aTGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAg                              <  2:167570/100‑1 (MQ=255)
                                                                aTGATTGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAg                              >  1:483657/1‑100 (MQ=255)
                                                                     tGAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTAccc                         >  2:1393900/1‑100 (MQ=255)
                                                                      gAAATATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCg                        >  1:356797/1‑100 (MQ=255)
                                                                        aaTATCGAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGc                      <  2:1097011/100‑1 (MQ=255)
                                                                              gAGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTTCCCGGCTATTGt                <  1:1036432/100‑1 (MQ=255)
                                                                               aGTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTa               <  1:1514966/100‑1 (MQ=255)
                                                                                gTTCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTaa              >  1:96803/1‑100 (MQ=255)
                                                                                 ttCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCTGCTATTGTAAc             <  2:548060/100‑1 (MQ=255)
                                                                                 ttCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAAc             <  2:1560741/100‑1 (MQ=255)
                                                                                  tCGCGCATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACg            <  1:223648/100‑1 (MQ=255)
                                                                                     cgcATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACgtgt         <  1:1425462/100‑1 (MQ=255)
                                                                                       cATACCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGttt       >  2:375472/1‑100 (MQ=255)
                                                                                          aCCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTc    <  2:1307891/100‑1 (MQ=255)
                                                                                          aCCACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTc    >  1:643242/1‑100 (MQ=255)
                                                                                           ccACAGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCg   <  1:1062194/100‑1 (MQ=255)
                                                                                            cacaGCGGCATGTCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGTAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGc  <  2:371664/100‑1 (MQ=255)
                                                                                            cacaGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGc  >  1:396091/1‑100 (MQ=255)
                                                                                               |                                                                                                
CAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATGCGGGCTTTGTTGTTAATGTCCATCTTCTTCCTGTCTCACGAAAATCCAGTACCCGGCTATTGTAACGTGTTTTTCGC  >  NC_000913/3040933‑3041124

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