Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,710 C→G 15.5% E48Q (GAA→CAA)  yddH ← flavin reductase like‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,7100CG15.5% 191.8 / 24.1 69E48Q (GAA→CAA) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (5/6);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAT  >  NC_000913/1533611‑1533802
                                                                                                   |                                                                                            
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                                                                               <  2:1896261/100‑1 (MQ=255)
               tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAt                                                                               >  1:2019194/1‑100 (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                              <  2:2430047/100‑1 (MQ=255)
                                                                cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaa                              <  1:392154/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:802777/5‑100 (MQ=255)
                                                                                      gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc        >  1:289511/5‑100 (MQ=255)
                                                                                      gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc        >  1:2121875/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)
                                                                                                   |                                                                                            
CTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAT  >  NC_000913/1533611‑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: