Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 916,170 T→A 11.7% N240K (AAT→AAA ycjF → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported916,1700TA11.7% 161.8 / 10.9 59N240K (AAT→AAAycjFconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/29);  new base A (4/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATTGCGGTCAGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTATTG  >  Exported/916100‑916239
                                                                      |                                                                     
tgatTGCGGTCAGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTtaattaat                                                                       >  1:774813/1‑71 (MQ=255)
 gatTGCGGTCAGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatc                                                                      >  1:2702434/1‑71 (MQ=255)
 gatTGCGGTCAGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatc                                                                      >  1:1787848/1‑71 (MQ=255)
        gTCAGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcatcgc                                                               >  1:2074957/1‑71 (MQ=255)
           aGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTtaattaat                                                                       >  1:453069/1‑60 (MQ=255)
           aGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCAc                                                            >  1:1256250/1‑71 (MQ=255)
           aGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCAc                                                            >  1:1314388/1‑71 (MQ=255)
             cccGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatc                                                                      <  1:1053024/59‑1 (MQ=255)
             cccGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGc                                                          >  1:1992889/1‑71 (MQ=255)
              ccGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatctcatc                                                                 >  1:2727834/1‑63 (MQ=255)
                gCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGt                                                       <  1:2654532/71‑1 (MQ=255)
                gCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGt                                                       <  1:3356074/71‑1 (MQ=255)
                  tGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGt                                                       >  1:3516058/1‑69 (MQ=255)
                  tGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTAt                                                     >  1:457189/1‑71 (MQ=255)
                  tGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTAt                                                     >  1:668624/1‑71 (MQ=255)
                    gCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATgg                                                   >  1:1883525/1‑71 (MQ=255)
                    gCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATg                                                    <  1:2740921/70‑1 (MQ=255)
                      gTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  <  1:3122261/54‑1 (MQ=255)
                      gTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  <  1:3387416/54‑1 (MQ=255)
                      gTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgc                                                                    <  1:1409352/52‑1 (MQ=255)
                      gTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgc                                                                    <  1:1980594/52‑1 (MQ=255)
                       ttGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAt                                                <  1:1857903/71‑1 (MQ=255)
                       ttGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAt                                                <  1:3386323/71‑1 (MQ=255)
                       ttGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAt                                                <  1:1645219/71‑1 (MQ=255)
                       ttGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAt                                                <  1:96240/71‑1 (MQ=255)
                       ttGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAt                                                <  1:2469179/71‑1 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:3487100/1‑52 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:790745/1‑52 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:3517685/1‑52 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:1762464/1‑52 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:1280216/1‑52 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:476868/1‑52 (MQ=255)
                        tGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAtt                                               >  1:2765795/1‑71 (MQ=255)
                           tCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGaa                                            <  1:3577546/71‑1 (MQ=255)
                           tCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGaa                                            <  1:3027172/71‑1 (MQ=255)
                           tCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGaa                                            <  1:675051/71‑1 (MQ=255)
                            cGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAAc                                           >  1:2723359/1‑71 (MQ=255)
                            cGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAAc                                           >  1:3070127/1‑71 (MQ=255)
                              atatGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTg                                         <  1:661366/71‑1 (MQ=255)
                              atatGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTg                                         <  1:265104/71‑1 (MQ=255)
                              atatGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTg                                         <  1:3551916/71‑1 (MQ=255)
                                 tGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTgggg                                      <  1:1837987/71‑1 (MQ=255)
                                 tGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTggg                                       <  1:12827/70‑1 (MQ=255)
                                   gagTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgca                                                                   >  1:3237655/3‑40 (MQ=255)
                                   gCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAatcgcat                                                                  >  1:3385077/1‑41 (MQ=255)
                                   gCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAtt                                               >  1:573515/1‑60 (MQ=255)
                                   gCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAtt                                               >  1:2982772/1‑60 (MQ=255)
                                      tttATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTa                                 >  1:1149616/1‑71 (MQ=255)
                                        tATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTa                                 >  1:1066468/1‑69 (MQ=255)
                                        tATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTa                                 >  1:1922347/1‑69 (MQ=255)
                                        tATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTAc                                <  1:2982962/70‑1 (MQ=255)
                                        tATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTAc                                <  1:2036312/70‑1 (MQ=255)
                                          tCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAtt                                               >  1:686054/1‑53 (MQ=255)
                                          tCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCAtt                                               >  1:635520/1‑53 (MQ=255)
                                            gCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCg                           <  1:3081853/71‑1 (MQ=255)
                                            gCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCg                           >  1:521986/1‑71 (MQ=255)
                                            gCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCg                           <  1:2854784/71‑1 (MQ=255)
                                            gCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCg                           >  1:2475080/1‑71 (MQ=255)
                                              cTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGtt                         >  1:1681985/1‑71 (MQ=255)
                                                ggCGCTATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTAt                                   <  1:1917688/59‑1 (MQ=255)
                                                    cAATCTGCGTTTAATTAATGGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTc                   >  1:3111699/1‑71 (MQ=255)
                                                    cAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTc                   >  1:772695/1‑71 (MQ=255)
                                                      cgatGCGATTAATTAAACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGt                    <  1:2816070/65‑1 (MQ=38)
                                                      aTCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTg                 <  1:1478637/71‑1 (MQ=255)
                                                       tCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTg                                              <  1:2531598/41‑1 (MQ=255)
                                                       tCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTg                                              <  1:632629/41‑1 (MQ=255)
                                                       tCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTg                                              <  1:1856678/41‑1 (MQ=255)
                                                       tCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGa                                             <  1:2264182/42‑1 (MQ=255)
                                                        atGCGATTAATTAAACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCt                  >  1:2693309/2‑68 (MQ=39)
                                                        atGCGATTAATTAAACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCt                  >  1:767954/2‑68 (MQ=39)
                                                          gCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGttt              >  1:3374279/1‑70 (MQ=255)
                                                          gCGATTAATTAAACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCg                           <  1:774043/57‑1 (MQ=39)
                                                          gCGATTAATTAAACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCg                           <  1:291124/57‑1 (MQ=39)
                                                           cGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGt                <  1:2246418/67‑1 (MQ=255)
                                                           cGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGt                <  1:3018314/67‑1 (MQ=255)
                                                            gattaattaaACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGtt                         >  1:2567437/3‑57 (MQ=255)
                                                            gattaattaaACGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGtt                         >  1:575385/3‑57 (MQ=255)
                                                             tttAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGc          <  1:751530/71‑1 (MQ=255)
                                                               taattaatCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTg        >  1:1312478/1‑71 (MQ=255)
                                                                 attaatCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGt      <  1:148309/71‑1 (MQ=255)
                                                                   taatCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTAt    >  1:2198134/1‑71 (MQ=255)
                                                                    aatCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTAtt   >  1:2807761/1‑71 (MQ=255)
                                                                    aatCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTAtt   >  1:3128752/1‑71 (MQ=255)
                                                                    aatCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTAt    >  1:3132972/1‑70 (MQ=255)
                                                                     atcgcatcgcCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTATTg  <  1:3368053/71‑1 (MQ=255)
                                                                      |                                                                     
TGATTGCGGTCAGCCCGCTGGCGTTGGTCGATATGGCGTTTATCGCCTGGCGCAATCTGCGTTTAATTAATCGCATCGCCACGCTGTATGGCATTGAACTGGGGTATTACAGCCGTTTGCGTCTGTTTAAGCTGGTATTG  >  Exported/916100‑916239

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