Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,707 G→T 15.3% P49T (CCT→ACT)  yddH ← flavin reductase like‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,7070GT15.3% 201.8 / 22.8 72P49T (CCT→ACT) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base T (5/6);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

TTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAT  >  NC_000913/1533608‑1533802
                                                                                                   |                                                                                               
ttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAgg                                                                                                 <  2:2179594/100‑1 (MQ=255)
 tgctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGt                                                                                                <  1:77850/100‑1 (MQ=255)
 tgctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACAACAATCGCCCCGCGTGGAgg                                                                                                 >  2:1657163/1‑99 (MQ=255)
  gctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGtt                                                                                               >  1:369962/1‑100 (MQ=255)
   ctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTc                                                                                              >  2:449832/1‑100 (MQ=255)
       aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTTAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                          <  1:878624/100‑1 (MQ=255)
         cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTc                                                                                        >  1:1324453/1‑100 (MQ=255)
          gCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTcc                                                                                       <  2:1328699/100‑1 (MQ=255)
             cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcg                                                                                    <  1:477285/100‑2 (MQ=255)
              gggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCg                                                                                   <  1:2551064/100‑1 (MQ=255)
               ggATAACGATGCCAAATTTACCGTTGTGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                  <  2:1881141/100‑4 (MQ=255)
               ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                  <  1:94072/100‑4 (MQ=255)
               ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                  <  1:255842/100‑4 (MQ=255)
                gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGc                                                                                 >  1:1222232/1‑100 (MQ=255)
                 aTAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtg                                                                                <  1:2528072/100‑6 (MQ=255)
                  tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAt                                                                               >  1:2019194/1‑100 (MQ=255)
                  tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAt                                                                               <  2:1896261/100‑1 (MQ=255)
                   aaCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtggc                                                                              <  1:1833868/100‑8 (MQ=38)
                    aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                             <  1:458756/100‑1 (MQ=255)
                     cGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                            >  1:1979073/1‑100 (MQ=255)
                      gATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAc                                                                           <  2:421252/100‑1 (MQ=255)
                         gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAg                                                                        <  2:1215178/100‑1 (MQ=255)
                          ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                       >  1:111478/1‑100 (MQ=255)
                          ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                       >  1:1739461/1‑100 (MQ=255)
                           cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                      <  2:1294151/100‑1 (MQ=255)
                             aaTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTg                                                                    <  2:1380019/100‑1 (MQ=255)
                              aTTTACCGTTATGTTCAATCCACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGc                                                                   <  2:1062002/100‑1 (MQ=255)
                                 tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAgt                                                                <  2:2153340/100‑2 (MQ=255)
                                   ccGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCa                                                              >  1:1061207/1‑100 (MQ=255)
                                           ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                      <  1:1478615/100‑1 (MQ=255)
                                            tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                     <  1:319514/100‑1 (MQ=255)
                                               atcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgcc                                                  >  1:2570940/1‑100 (MQ=255)
                                                  aaCTGCCTGGTCCATGTCGATTTAGCTACTACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCGTAATGTTAcgac                                                  >  2:1771570/1‑95 (MQ=255)
                                                    cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGATGCAGCCATAATGTTACGCCGCTgg                                             >  1:629883/1‑100 (MQ=255)
                                                    cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAATGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                             >  2:1375303/1‑100 (MQ=255)
                                                     tCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTggg                                            >  2:948461/1‑100 (MQ=255)
                                                      cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCCCCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                           >  2:1495038/1‑100 (MQ=255)
                                                      cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                           >  1:482432/1‑100 (MQ=255)
                                                          ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                       >  2:1971225/1‑100 (MQ=255)
                                                          ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                       <  1:2196096/100‑1 (MQ=255)
                                                           gTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCCAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCt                                      >  1:1897580/1‑100 (MQ=255)
                                                           gTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCt                                      <  2:2155701/100‑1 (MQ=255)
                                                             ccATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCa                                    <  2:1691747/100‑1 (MQ=255)
                                                               aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                  <  2:103251/100‑1 (MQ=255)
                                                               aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                  <  2:103298/100‑1 (MQ=255)
                                                               aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                  >  1:2560277/1‑100 (MQ=255)
                                                                tGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCg                                 <  2:331906/100‑1 (MQ=255)
                                                                tGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCg                                 >  1:408107/1‑100 (MQ=255)
                                                                 gTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                <  1:948466/100‑1 (MQ=255)
                                                                 gTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                >  2:286881/1‑100 (MQ=255)
                                                                 gTCGATTTATCTACCACAATAGCCAAGCGTGGAGGTTCAAACTCCACCGGCAATGACCCGGCTGCAGCCATAATGTTACGCCGCTGGGATCGCTCATCGa                                >  2:496186/1‑100 (MQ=255)
                                                                   cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaa                              <  1:392154/100‑1 (MQ=255)
                                                                   cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaa                              <  2:2430047/100‑1 (MQ=255)
                                                                           cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                      >  1:1128915/1‑100 (MQ=255)
                                                                           cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                      <  1:1461181/100‑1 (MQ=255)
                                                                            tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATc                     >  1:202210/1‑100 (MQ=255)
                                                                              ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                   >  1:1463207/1‑100 (MQ=255)
                                                                               cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                  <  2:573103/100‑1 (MQ=255)
                                                                                acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGa                 >  1:1707243/1‑100 (MQ=255)
                                                                                   aTCGCCACGCGTGGAGGTTCAGACTCGCCCGGCATTGACCAGGCTGCAGCCATAATGTTCCGCCGCTGGGATTGCTCATCGAATCTGGTGATCAAGACCg              <  2:2369522/100‑1 (MQ=255)
                                                                                     cGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGtt            >  2:1498199/1‑100 (MQ=255)
                                                                                      gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg           >  1:162021/1‑100 (MQ=255)
                                                                                      gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg           >  1:1980086/1‑100 (MQ=255)
                                                                                       ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg          <  1:2304521/100‑1 (MQ=255)
                                                                                         gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc        >  1:289511/5‑100 (MQ=255)
                                                                                         gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc        >  1:2121875/5‑100 (MQ=255)
                                                                                         gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc        >  1:802777/5‑100 (MQ=255)
                                                                                         aCGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc        <  1:701127/100‑1 (MQ=255)
                                                                                          ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc       >  2:2167822/4‑100 (MQ=255)
                                                                                          ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc       >  2:2089819/4‑100 (MQ=255)
                                                                                             gTGGAGGTTCAAACTCCACCGGCAGTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATg    <  2:1635600/100‑1 (MQ=255)
                                                                                              tGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa   >  1:1132998/1‑100 (MQ=255)
                                                                                              tGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa   <  1:1657172/100‑1 (MQ=255)
                                                                                               ggAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAt  >  2:1900464/1‑100 (MQ=255)
                                                                                                   |                                                                                               
TTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAT  >  NC_000913/1533608‑1533802

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