Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,165,181 C→A G221C (GGC→TGC)  yehU ← predicted sensory kinase in two‑component system with YehT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,165,1810CA97.4% 127.7 / NA 38G221C (GGC→TGC) yehUpredicted sensory kinase in two‑component system with YehT
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (21/16);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCAGCTCCTGATACAGCACCTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCCTTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACATCGCGCGTTT  >  REL606/2165109‑2165253
                                                                        |                                                                        
tatCCAGCTCCTGATACAGCACCTGAGCCACTTTCAGGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTc                                                                       <  1:368837/76‑1 (MQ=255)
 atCCAGCTCCTGATACAGCACCTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTcc                                                                      <  1:209454/76‑1 (MQ=255)
        tCCTGATACAGCACCTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAgg                                                               >  2:958688/1‑76 (MQ=255)
           tGATACAGCACCTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTg                                                            >  2:7363/1‑76 (MQ=255)
             aTACAGCACCTGAGCCACCTTCATGCTGTTCACTTCGTTAAACCCCTGTGGCAAAATGCCTTCCGTCGAGGCTGcc                                                          <  1:1347540/76‑1 (MQ=255)
              tACAGCACCTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGccc                                                         >  2:16323/1‑75 (MQ=255)
                      cTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAgc                                                 >  1:985303/1‑76 (MQ=255)
                       tGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAgcg                                                >  2:276342/1‑76 (MQ=255)
                       tGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAgcg                                                >  1:553652/1‑76 (MQ=255)
                            cACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTg                                           <  2:514532/76‑1 (MQ=255)
                            cACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTg                                           <  1:650234/76‑1 (MQ=255)
                                ttCATGCTGTTCACTTCGTTAGCCCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAg                                       <  1:754705/76‑1 (MQ=255)
                                ttCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAGATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAg                                       <  1:1081659/76‑1 (MQ=255)
                                ttCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAg                                       <  2:230916/76‑1 (MQ=255)
                                ttCATGCTGTTCACNTCGTTAACCCCCTCTCGCANAATGCGTTNCGTCGAGGCTGCCACTTTCCGCGCAGTGGCAg                                       <  2:1172749/76‑1 (MQ=255)
                                 tCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGa                                      >  1:964224/1‑76 (MQ=255)
                                  cATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGaa                                     <  2:1419451/76‑1 (MQ=255)
                                    tGCTATTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGaaaa                                   >  2:1404685/1‑76 (MQ=255)
                                     gCTGTTCACTTCGTTAACCCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGaaaaa                                  <  1:591222/76‑1 (MQ=255)
                                     gCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGaaaaa                                  <  1:636598/76‑1 (MQ=255)
                                       tGTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGACNCTNCCACTTTCAGCGCANTGGCAGAAAAAGc                                >  2:1286259/1‑76 (MQ=255)
                                        gTTCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGcc                               <  2:427765/76‑1 (MQ=255)
                                         ttCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCg                              <  1:1213454/76‑1 (MQ=255)
                                          tCACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGa                             >  2:56995/1‑76 (MQ=255)
                                           cACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGaa                            >  1:1088323/1‑76 (MQ=255)
                                           cACTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGaa                            <  1:1212354/76‑1 (MQ=255)
                                             cTTCGTTAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAgg                          >  1:362143/1‑75 (MQ=255)
                                                  ttAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAtt                     >  1:879622/1‑76 (MQ=255)
                                                   tAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAttt                    >  2:1092374/1‑76 (MQ=255)
                                                   tAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAttt                    >  2:483746/1‑76 (MQ=255)
                                                   tAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAttt                    >  2:703791/1‑76 (MQ=255)
                                                   tAAACCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAttt                    >  2:967361/1‑76 (MQ=255)
                                                    aacccccTGTCGCAAAATGCATCCGGCCGGGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAtttt                   <  1:1139203/73‑1 (MQ=255)
                                                     aaCCCCTGTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTAttttt                  >  1:143920/1‑76 (MQ=255)
                                                            gTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACAt           <  1:148902/76‑1 (MQ=255)
                                                            gTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACAt           <  1:1145545/76‑1 (MQ=255)
                                                            gTCGCAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCAGAAGTGTATTTTTCAAACAt           >  2:101699/1‑76 (MQ=255)
                                                                cAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACATcgcg       <  1:784385/76‑1 (MQ=255)
                                                                cAAAATGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACATcgcg       >  2:274145/1‑76 (MQ=255)
                                                                   aaTGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACATCGCGCGt    >  2:880274/1‑76 (MQ=255)
                                                                     tGCATTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGAATTTTTCAAACATCGCGCGttt  <  1:7363/76‑1 (MQ=255)
                                                                        |                                                                        
TATCCAGCTCCTGATACAGCACCTGAGCCACTTTCATGCTGTTCACTTCGTTAAACCCCTGTCGCAAAATGCCTTCCGTCGAGGCTGCCACTTTCAGCGCAGTGGCAGAAAAAGCCGAAGTGTATTTTTCAAACATCGCGCGTTT  >  REL606/2165109‑2165253

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