Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,897,276 C→T Q21* (CAA→TAA)  yqeA → predicted amino acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,897,2760CT96.7% 213.1 / NA 62Q21* (CAA→TAA) yqeApredicted amino acid kinase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (38/21);  total (41/21)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGAGTGAGGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGT  >  REL606/2897203‑2897346
                                                                         |                                                                      
tgagtgagGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTaa                                                                      <  2:1286035/76‑1 (MQ=255)
tgagtgagGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTaa                                                                      >  1:999984/1‑76 (MQ=255)
   gtgagGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATg                                                                   <  1:1559528/76‑1 (MQ=255)
    tgagGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGa                                                                  >  1:1035851/1‑76 (MQ=255)
     gagGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGaa                                                                 <  2:380777/76‑1 (MQ=255)
     gagGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGaa                                                                 >  1:68260/1‑76 (MQ=255)
      agGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGaaa                                                                <  2:561163/76‑1 (MQ=255)
         aaaaaTAAAAAAAAAAATTCTGCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCg                                                             <  1:1218855/76‑1 (MQ=255)
           aaaTGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGt                                                           <  2:323292/76‑1 (MQ=255)
           aaaTGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGt                                                           >  2:1201159/1‑76 (MQ=255)
            aaTGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTa                                                          >  2:944835/1‑76 (MQ=255)
             aTGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaa                                                         >  2:330010/1‑76 (MQ=255)
              tGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaaa                                                        <  2:375619/76‑1 (MQ=255)
              tGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaaa                                                        <  2:1499671/76‑1 (MQ=255)
              tGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaaa                                                        <  1:336660/76‑1 (MQ=255)
              tGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaaa                                                        >  2:637826/1‑76 (MQ=255)
              tGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGGaaa                                                        >  2:693983/1‑76 (MQ=255)
               gAGTAAGAAAGTTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaaaa                                                       <  1:1272328/76‑1 (MQ=255)
                aGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTaaaaa                                                      <  2:729083/76‑1 (MQ=255)
                 gTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAAt                                                     >  2:119585/1‑76 (MQ=255)
                   aaGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTa                                                   >  2:268689/1‑76 (MQ=255)
                       aaaTTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTtc                                               >  2:1702181/1‑76 (MQ=255)
                          ttGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCa                                            >  2:635295/1‑76 (MQ=255)
                          ttGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCa                                            >  2:1441864/1‑76 (MQ=255)
                           tGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAg                                           >  2:892582/1‑76 (MQ=255)
                            gTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAgg                                          >  1:1694754/1‑76 (MQ=255)
                            gTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAgg                                          >  1:1448728/1‑76 (MQ=255)
                             ttCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGc                                         >  1:1246570/1‑76 (MQ=255)
                              tctcGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCa                                        >  2:1436733/1‑76 (MQ=255)
                              tctcGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCa                                        >  1:723018/1‑76 (MQ=255)
                               ctcGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCaa                                       >  2:166562/1‑76 (MQ=255)
                               ctcGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCaa                                       >  2:715723/1‑76 (MQ=255)
                               ctcGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGGAAATGGAAGCGGTAAAAATTACTTCTCAGGCaa                                       >  2:283092/1‑76 (MQ=255)
                                tcGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAAt                                      >  2:1235611/1‑76 (MQ=255)
                                 cGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAAtt                                     >  2:1141707/1‑76 (MQ=255)
                                    ccTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAAttgtt                                  >  2:1323388/1‑76 (MQ=255)
                                      tGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGa                                <  1:56590/76‑1 (MQ=255)
                                         gCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGAttt                             >  1:822996/1‑76 (MQ=255)
                                         gCGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGAttt                             >  1:768110/1‑76 (MQ=255)
                                          cGGGACTGCGCTGGGCGCTGACCTGGCCGGGTAAATGAAAGGGGTAAAACTTACTTCTCCGGAAATTGTTGATTTa                            >  2:129406/1‑76 (MQ=255)
                                          cGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTa                            <  2:454637/76‑1 (MQ=255)
                                          cGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTa                            >  2:129397/1‑76 (MQ=255)
                                          cGGGAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTa                            <  2:663765/76‑1 (MQ=255)
                                           gggAATGCGCTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTaa                           <  2:534044/76‑1 (MQ=255)
                                                   gcTGGGCGATGACCTGGCCGGGTAAATGACAGCGGTAAAAATTACTTCTAAGGCACTTGTTGATTCAATTGCCAAg                   >  2:583645/1‑76 (MQ=255)
                                                   gcTGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCCCAGGCAATTGTTGATTTAATTGCCCAg                   >  2:664623/1‑76 (MQ=255)
                                                     tGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAggg                 <  2:295188/76‑1 (MQ=255)
                                                     tGGGCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAggg                 >  2:1281729/1‑76 (MQ=255)
                                                      gggCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGa                >  1:767405/1‑76 (MQ=255)
                                                      gggCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGa                >  1:834379/1‑76 (MQ=255)
                                                      gggCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGa                >  2:2682/1‑76 (MQ=255)
                                                        gCGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACa              <  2:957042/76‑1 (MQ=255)
                                                         cGATGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACAt             >  2:1316887/1‑76 (MQ=255)
                                                           aTGACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGa           >  2:889416/1‑76 (MQ=255)
                                                             gACCTGGCCGGGTAAATGAAAGCGGTACAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAg         >  2:634527/1‑76 (MQ=255)
                                                             gACCTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAg         >  2:204904/1‑76 (MQ=255)
                                                                cTGGCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTa      <  2:1448727/76‑1 (MQ=255)
                                                                  ggCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTTCCCAGGGACATGAAGTTATc    <  2:1363149/76‑1 (MQ=255)
                                                                  ggCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATc    >  1:1545729/1‑76 (MQ=255)
                                                                   gCCGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCg   <  2:370821/76‑1 (MQ=255)
                                                                    ccGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGt  >  1:727864/1‑76 (MQ=255)
                                                                    ccGGGTAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGt  <  2:1135395/76‑1 (MQ=255)
                                                                         |                                                                      
TGAGTGAGGAAAAATGAGTAAGAAAATTGTTCTCGCCCTGGGCGGGAATGCGCTGGGCGATGACCTGGCCGGGCAAATGAAAGCGGTAAAAATTACTTCTCAGGCAATTGTTGATTTAATTGCCCAGGGACATGAAGTTATCGT  >  REL606/2897203‑2897346

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