Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,224,879 Δ1 bp coding (1005/1434 nt) ygjI → putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,224,8780G.99.2% 557.0 / ‑2.3 123coding (1004/1434 nt)ygjIputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (71/51);  minor base A (0/1);  total (71/52)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 139 total aligned reads displayed.

TTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCGATGCTCGGTTCCAATACCGTG  >  NC_000913/3224803‑3224948
                                                                           |                                                                       
ttttCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCcggcgcg                                                                         <  1:3467408/76‑1 (MQ=255)
 tttCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCcggcgcgc                                                                        >  1:1392381/1‑75 (MQ=255)
 tttCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCcggcgcgc                                                                        <  1:5351641/76‑2 (MQ=255)
  ttcttcTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCcggcgcgcg                                                                       <  2:5262396/76‑3 (MQ=255)
   tcttctCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCcggcgcgcgc                                                                      >  2:600525/1‑73 (MQ=255)
   tcttctCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCcggcgcgcgc                                                                      <  2:4338392/76‑4 (MQ=255)
      tctcCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGc                                                                  >  2:2537600/1‑76 (MQ=255)
        tcCGAGATGCCGGAAGGCATCTTTGGTAATAAAACCGACTAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGt                                                                <  2:1276398/76‑1 (MQ=255)
         ccGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTg                                                               >  2:788059/1‑76 (MQ=255)
         ccGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTg                                                               >  2:2966342/1‑76 (MQ=255)
         ccGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTg                                                               >  1:524764/1‑76 (MQ=255)
           gAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGa                                                             >  2:3214934/1‑76 (MQ=255)
           gAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGa                                                             >  2:1159713/1‑76 (MQ=255)
            aaaTCCCGGACGGCGTCTTGGGTAAGAAAGCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGAt                                                            <  2:5796977/76‑1 (MQ=255)
            aaaTCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGAt                                                            >  1:3453028/1‑76 (MQ=255)
             aaTCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATc                                                           >  2:1667365/1‑76 (MQ=255)
              aTCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATcc                                                          <  1:1073077/76‑1 (MQ=255)
                cccGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAg                                                        <  2:5358791/76‑1 (MQ=255)
                cccGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAg                                                        <  1:4688322/76‑1 (MQ=255)
                 ccGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGt                                                       >  2:1648380/1‑76 (MQ=255)
                 ccGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGt                                                       >  2:2057678/1‑76 (MQ=255)
                  cGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGtt                                                      >  1:3883404/1‑76 (MQ=255)
                    gAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTcc                                                    >  2:2576465/1‑76 (MQ=255)
                    gAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGAGCGCAGCGTGGATCCAGTTcc                                                    <  1:899252/76‑1 (MQ=255)
                     aaGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCt                                                   <  2:680947/76‑1 (MQ=255)
                     aaGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCt                                                   <  2:105661/76‑1 (MQ=255)
                        gCATCTTTGGTAAGAAAACCGTCGAACTGAACGACAACGGCGTTCCGGCGCGCGCAGAGTGGCTCCAGTTCCTGac                                                >  1:269716/1‑75 (MQ=255)
                        gCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGAt                                                <  1:1204843/76‑1 (MQ=255)
                         cATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATc                                               >  1:5323099/1‑76 (MQ=255)
                         cATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATc                                               >  1:4566359/1‑76 (MQ=255)
                         cATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATc                                               >  2:379625/1‑76 (MQ=255)
                            cTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTc                                            >  2:2449881/1‑76 (MQ=255)
                              ttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCAt                                          >  2:4251190/1‑76 (MQ=255)
                               tGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATc                                         >  1:160707/1‑76 (MQ=255)
                                 gTAAGAAACCCGACGACCTGAACGAAACCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATccc                                       <  2:3023600/76‑1 (MQ=255)
                                 gTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATccc                                       <  2:4210717/76‑1 (MQ=255)
                                 gTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATccc                                       <  2:338880/76‑1 (MQ=255)
                                 gTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATccc                                       <  1:2507881/76‑1 (MQ=255)
                                  tAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCg                                      <  2:3910957/76‑1 (MQ=255)
                                   aaGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGc                                     >  1:5419958/1‑76 (MQ=255)
                                   aaGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGc                                     <  2:2129784/76‑1 (MQ=255)
                                    aGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCt                                    >  1:4144355/1‑76 (MQ=255)
                                    aGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCt                                    >  2:2445590/1‑76 (MQ=255)
                                     gAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCtg                                   <  2:3661837/76‑1 (MQ=255)
                                      aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCCTCCCGCtgt                                  >  2:1450117/1‑75 (MQ=255)
                                      aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCtga                                  <  1:3114362/76‑1 (MQ=255)
                                      aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCtga                                  >  2:4332798/1‑76 (MQ=255)
                                      aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCtga                                  <  1:2475330/76‑1 (MQ=255)
                                      aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCtga                                  >  1:2329485/1‑76 (MQ=255)
                                       aaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCtgat                                 <  1:5133986/76‑1 (MQ=255)
                                       aaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCACTTCCCGATCGTCCTCCCGCtgac                                 >  2:1420684/1‑75 (MQ=255)
                                         aCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGa                                                             >  1:3952520/1‑46 (MQ=39)
                                           cGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGAtt                             <  1:677367/76‑1 (MQ=255)
                                            gTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGACGATTa                            >  2:1808393/1‑76 (MQ=255)
                                             tCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTAt                           >  1:1522168/1‑76 (MQ=255)
                                              ccaaCTCAACGAAAACGGCGGTCCGACGCGCGCAGCGCGGATCCAGTTCCTGATCGTCACCCCGCTGACGCTTATc                         >  2:3225786/3‑76 (MQ=255)
                                               gAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcc                        <  2:5621229/76‑1 (MQ=255)
                                               gAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcc                        <  1:1667365/76‑1 (MQ=255)
                                               gAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcc                        >  2:4496305/1‑76 (MQ=255)
                                               gAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcc                        >  1:3837101/1‑76 (MQ=255)
                                               gAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcc                        <  1:3065/76‑1 (MQ=255)
                                                aaCTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATccc                       <  1:1902039/76‑2 (MQ=255)
                                                aaCTGAACGAAAACGGCGTGCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATccc                       <  2:3696411/76‑2 (MQ=255)
                                                  cTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcccga                     >  1:3647079/1‑73 (MQ=255)
                                                  cTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcccga                     <  2:2184929/76‑4 (MQ=255)
                                                     aaCGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGc                   >  1:4889204/1‑76 (MQ=255)
                                                      gcGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCt                  <  2:1971012/75‑1 (MQ=255)
                                                      aCGAAAACGGCTTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCt                  <  1:2802110/76‑1 (MQ=255)
                                                       cGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTc                 <  1:1681190/76‑1 (MQ=255)
                                                       cGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTc                 <  1:556122/76‑1 (MQ=255)
                                                       cGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTc                 >  2:1619509/1‑76 (MQ=255)
                                                       cGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTc                 <  2:3234411/76‑1 (MQ=255)
                                                       cGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTc                 <  2:730944/76‑1 (MQ=255)
                                                        gAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCg                >  1:1375132/1‑76 (MQ=255)
                                                        gAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCg                >  2:1498926/1‑76 (MQ=255)
                                                         aaaaCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATTGTCATCCCGCTGATGATTATCCCGATGCTCgg               >  2:4482110/1‑76 (MQ=255)
                                                         aaaaCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCgg               <  1:2958236/76‑1 (MQ=255)
                                                         aaaaCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCgg               >  1:5422994/1‑76 (MQ=255)
                                                          aaaCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATcccg                      >  1:2522965/1‑65 (MQ=255)
                                                           aaCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGtt             <  2:1989991/76‑1 (MQ=255)
                                                           aaCGGCGTGCCGGTGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGtt             <  1:3337836/76‑1 (MQ=255)
                                                            aCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTc            >  2:4448770/1‑76 (MQ=255)
                                                            aCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTc            >  1:5109949/1‑76 (MQ=255)
                                                             cGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTcc           >  2:3828108/1‑76 (MQ=255)
                                                              ggCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCa          >  2:2529325/1‑76 (MQ=255)
                                                              ggCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCa          >  2:3385526/1‑76 (MQ=255)
                                                              ggCGTTCCGGCGCGCGCAGCGTGGATCCAGCTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCa          >  2:1755949/1‑76 (MQ=255)
                                                                cGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAAt        <  1:4609159/76‑1 (MQ=255)
                                                                 ggtCCGGCGAGCGCAGCGCGGATCCGGTTCCTGATCGTCTTCCCGCTGATGATTATACCGATGCTGGGATCCGATa       >  1:4783817/3‑76 (MQ=255)
                                                                 gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATa       >  1:4497690/1‑76 (MQ=255)
                                                                 gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATa       >  2:3471406/1‑76 (MQ=255)
                                                                 gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATa       >  2:2150428/1‑76 (MQ=255)
                                                                 gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATa       <  1:3585528/76‑1 (MQ=255)
                                                                 gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATa       >  1:2315928/1‑76 (MQ=255)
                                                                 gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGACGATTATCCCGATGCTCGGTTCCAATa       >  1:1094085/1‑76 (MQ=255)
                                                                   tCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATAcc     >  2:5693238/1‑76 (MQ=255)
                                                                      ggcgcgcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTg  >  1:3846351/1‑76 (MQ=255)
                                                                      ggcgcgcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTg  >  1:1368315/1‑76 (MQ=255)
                                                                      ggcgcgcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTg  <  2:5067130/76‑1 (MQ=255)
                                                                      ggcgcgcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTg  >  2:1262703/1‑76 (MQ=255)
                                                                           |                                                                       
TTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCGATGCTCGGTTCCAATACCGTG  >  NC_000913/3224803‑3224948

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