Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 318,568 T→C Q140Q (CAA→CAG yahN ← neutral amino‑acid efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606318,5680TC97.5% 159.4 / ‑3.7 40Q140Q (CAA→CAGyahNneutral amino‑acid efflux system
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (15/24);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CCATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGT  >  REL606/318495‑318640
                                                                         |                                                                        
ccATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTg                                                                        <  2:39274/76‑1 (MQ=255)
  aTTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCg                                                                      >  2:137956/1‑76 (MQ=255)
      aaCGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGtt                                                                  <  2:1028512/76‑1 (MQ=255)
            cccATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTagagag                                                            <  2:489678/76‑1 (MQ=255)
             ccATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAgagaga                                                           >  2:86341/1‑76 (MQ=255)
             ccATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAgagaga                                                           <  2:1123636/76‑1 (MQ=255)
              cATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGAt                                                          >  2:1341677/1‑76 (MQ=255)
                   tGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGta                                                     <  1:124572/76‑1 (MQ=255)
                   tGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGta                                                     <  2:1583575/76‑1 (MQ=255)
                    ggTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaa                                                    <  2:1477863/76‑1 (MQ=255)
                    ggTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaa                                                    >  2:654785/1‑76 (MQ=255)
                    ggTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaa                                                    <  1:674923/76‑1 (MQ=255)
                    ggTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaa                                                    <  2:250305/76‑1 (MQ=255)
                     gtgtTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaat                                                   <  2:29783/76‑1 (MQ=255)
                     gtgtTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaat                                                   >  1:1193120/1‑76 (MQ=255)
                      tgtTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaatt                                                  <  1:83873/76‑1 (MQ=255)
                      tgtTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaatt                                                  >  1:1245873/1‑76 (MQ=255)
                       gtTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaatta                                                 <  2:542992/76‑1 (MQ=255)
                        tttCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGtaattaa                                                <  2:168562/76‑1 (MQ=255)
                           cGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATcc                                             >  2:299388/1‑76 (MQ=255)
                           cGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATcc                                             <  2:1634380/76‑1 (MQ=255)
                            ggCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCcg                                            <  1:308823/76‑1 (MQ=255)
                             gCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCcgc                                           <  1:339107/76‑1 (MQ=255)
                               aTTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCgcgg                                         >  2:1144741/1‑76 (MQ=255)
                                 ttAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCgcggcg                                       >  1:1184386/1‑76 (MQ=255)
                                  tAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGa                                      <  2:1537650/76‑1 (MQ=255)
                                   aaTGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGaa                                     >  1:554644/1‑76 (MQ=255)
                                    aTGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGaaa                                    >  1:704921/1‑76 (MQ=255)
                                     tGTTACTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGaaaa                                   >  1:1171712/1‑76 (MQ=255)
                                        tACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAg                                >  2:1440105/1‑76 (MQ=255)
                                          cTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGAc                              >  2:833144/1‑76 (MQ=255)
                                          cTGAGAAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGAc                              <  2:1659162/76‑1 (MQ=255)
                                          cTGAAAAAATACTGATAAAAAATAAAACGGTATGAGGGTTAGAGAGATCGGTAATTAATCTGCGGCGAAAAAAGaa                              >  1:481692/1‑75 (MQ=255)
                                             agaAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATa                           <  2:1376151/76‑1 (MQ=255)
                                             agaAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATa                           <  2:1172693/76‑1 (MQ=255)
                                             agaAAATACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATa                           <  2:619295/76‑1 (MQ=255)
                                               aaaaTACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACGTAcc                         <  1:823528/76‑1 (MQ=255)
                                                aaaTACTGATAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCa                        <  1:133178/76‑1 (MQ=255)
                                                        aTAAAAAATAAAACGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGgcgc                >  2:1110811/1‑76 (MQ=255)
                                                                 aaaaCGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGgtt       >  1:1239402/1‑76 (MQ=255)
                                                                     cGGTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGttgttg   >  1:1339165/1‑76 (MQ=255)
                                                                      ggTCTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTtgttgt  <  2:1563081/76‑1 (MQ=255)
                                                                         |                                                                        
CCATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGT  >  REL606/318495‑318640

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