Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,881,536 T→G D7A (GAT→GCT)  atpE ← F0F1 ATP synthase subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,881,5360TG94.9% 224.4 / 0.9 59D7A (GAT→GCT) atpEF0F1 ATP synthase subunit C
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (26/30);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GAGGATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGATCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAAT  >  REL606/3881463‑3881602
                                                                         |                                                                  
gAGGATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGcc                                                                  >  1:96831/1‑76 (MQ=255)
 aGGATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCa                                                                 >  1:1849434/1‑76 (MQ=255)
 aGGATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCa                                                                 <  1:850287/76‑1 (MQ=255)
 aGGATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCa                                                                 <  1:554595/76‑1 (MQ=255)
  ggATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCat                                                                <  1:1681209/76‑1 (MQ=255)
    aTGCCGATCCCGATCGCAGCCCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCatat                                                              <  2:879887/76‑1 (MQ=255)
    aTGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCatat                                                              <  2:563263/76‑1 (MQ=255)
     tGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATAtt                                                             <  2:1502345/76‑1 (MQ=255)
     tGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATAtt                                                             >  2:674034/1‑76 (MQ=255)
       ccGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCa                                                           >  2:1342449/1‑76 (MQ=255)
          aTACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGt                                                        >  2:1727231/1‑76 (MQ=255)
          aTACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGt                                                        >  2:543373/1‑76 (MQ=255)
           tACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGtt                                                       <  1:1075733/76‑1 (MQ=255)
            aCCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGttt                                                      >  2:1662094/1‑76 (MQ=255)
            aCCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGttt                                                      >  2:218404/1‑76 (MQ=255)
              cGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTc                                                    <  2:181549/76‑1 (MQ=255)
               gATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTcc                                                   <  1:1542313/76‑1 (MQ=255)
               gATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTcc                                                   >  1:698669/1‑76 (MQ=255)
                aTCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCa                                                  <  1:250867/76‑1 (MQ=255)
                  cGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATg                                                >  2:1907446/1‑76 (MQ=255)
                    cagcaCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGAc                                              <  1:1612545/76‑1 (MQ=255)
                    cagcaCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGAc                                              >  2:1191093/1‑76 (MQ=255)
                     agcaCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACa                                             >  2:1560897/1‑76 (MQ=255)
                     agcaCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACa                                             >  2:555216/1‑76 (MQ=255)
                      gcaCCGATTGCCGCCAGGCCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAg                                            <  2:759293/76‑1 (MQ=255)
                      gcaCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAg                                            <  2:414845/76‑1 (MQ=255)
                      gcaCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAg                                            <  2:759279/76‑1 (MQ=255)
                        aCCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGtc                                          <  2:620383/76‑1 (MQ=255)
                          cGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGtctc                                        >  2:78228/1‑76 (MQ=255)
                           gATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTcc                                       <  1:283319/76‑1 (MQ=255)
                            aTTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCa                                      <  2:491999/76‑1 (MQ=255)
                                ccgccAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGTCAGTCTCCCgttt                                  >  2:612479/1‑76 (MQ=255)
                                ccgccAGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAgttt                                  <  2:723293/76‑1 (MQ=255)
                                     aGACCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTc                             >  1:1936343/1‑76 (MQ=255)
                                       aCCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAg                           >  2:108321/1‑76 (MQ=255)
                                       aCCCATCATCACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAg                           <  1:13739/76‑1 (MQ=255)
                                           atcatcACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTaa                       >  1:1721380/1‑76 (MQ=255)
                                           atcatcACAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTaa                       >  1:562144/1‑76 (MQ=255)
                                            tcatcaCAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTaaa                      <  1:489541/76‑1 (MQ=255)
                                             catcaCAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTaaaa                     <  2:1210526/76‑1 (MQ=255)
                                              atcaCAGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAAc                    <  2:146434/76‑1 (MQ=255)
                                                cacaGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACgt                  <  2:1176293/76‑1 (MQ=255)
                                                   aGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACgtagt               >  2:1509171/1‑76 (MQ=255)
                                                   aGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACgtagt               <  1:978496/76‑1 (MQ=255)
                                                   aGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACgtagt               >  1:29919/1‑76 (MQ=255)
                                                   aGCGGCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACgtagt               <  1:218404/76‑1 (MQ=255)
                                                      ggCAGCCATGTACAGCAGATCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTagtagt            >  1:735179/1‑76 (MQ=255)
                                                      ggCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTagtagt            >  1:1269274/1‑76 (MQ=255)
                                                       gCAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAgtg           >  1:1036839/1‑76 (MQ=255)
                                                        cAGCCATGTACAGCAGATCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAgtgt          >  1:1611440/1‑76 (MQ=255)
                                                        cAGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAgtgt          >  2:825187/1‑76 (MQ=255)
                                                         aGCCATGTACAGCAGATCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGtt         >  1:1574031/1‑76 (MQ=255)
                                                         aGCCATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGtt         <  2:1211940/76‑1 (MQ=255)
                                                            cATGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGt      <  2:218380/76‑1 (MQ=255)
                                                             aTGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTa     <  2:1919488/76‑1 (MQ=255)
                                                             aTGTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTa     >  1:828433/1‑76 (MQ=255)
                                                               gTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTaaa   <  2:135568/76‑1 (MQ=255)
                                                               gTACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTaaa   <  2:90613/76‑1 (MQ=255)
                                                                tACAGCAGAGCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAAt  >  1:442730/1‑76 (MQ=255)
                                                                         |                                                                  
GAGGATGCCGATACCGATCGCAGCACCGATTGCCGCCAGACCCATCATCACAGCGGCAGCCATGTACAGCAGATCCATATTCAGGTTTTCCATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAAT  >  REL606/3881463‑3881602

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