Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,716,825 C→A 9.1% intergenic (+12/‑1) atpH → / → atpA F1 sector of membrane‑bound ATP synthase, delta subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,716,8250CA 9.1% 199.9 / 4.1 65intergenic (+12/‑1)atpH/atpAF1 sector of membrane‑bound ATP synthase, delta subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (24/35);  new base A (4/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GATGGCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAATGTTGTGAGTG  >  W3110S.gb/3716760‑3716895
                                                                 |                                                                      
gATGGCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGc                                                                    <  1:1936508/70‑1 (MQ=255)
gATGGCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGc                                                                    >  1:1122272/1‑70 (MQ=255)
gATGGCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGc                                                                    >  1:7151/1‑70 (MQ=255)
    gCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTg                                                               <  1:3756242/71‑1 (MQ=255)
    gCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTg                                                               <  1:3786164/71‑1 (MQ=255)
     cAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGa                                                              <  1:624223/71‑1 (MQ=255)
     cAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGa                                                              <  1:3488867/71‑1 (MQ=255)
        cGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGc                                                                        >  1:1782795/1‑58 (MQ=255)
            cgcgGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCAc                                                       >  1:1286178/1‑71 (MQ=255)
            cgcgGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCAc                                                       >  1:750283/1‑71 (MQ=255)
            cgcgGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCAc                                                       >  1:2203786/1‑71 (MQ=255)
             gcgGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGTGCATGCAACTGAATTCCAcc                                                      <  1:3433536/71‑1 (MQ=255)
                 tcgtcTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGaaa                                                  <  1:2871582/71‑1 (MQ=255)
                  cgtcTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAAt                                                 >  1:2532122/1‑71 (MQ=255)
                   gtcTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAAt                                                 <  1:1698061/70‑1 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGc                                             >  1:1674101/1‑71 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGc                                             >  1:3431046/1‑71 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGc                                             >  1:2314370/1‑71 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGAATGCAACTGAATTCCACCGAAATCAGc                                             >  1:2453651/1‑71 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGAATGCAACTGAATTCCACCGAAATCAGc                                             >  1:1796068/1‑71 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGAATGCAACTGAATTCCACCGAAATCAGc                                             >  1:2030320/1‑71 (MQ=255)
                      ttGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGAATGCAACTGAATTCCACCGAAATCAGc                                             >  1:2025794/1‑71 (MQ=255)
                      ttGAGCGCCTTCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGc                                             >  1:771856/1‑70 (MQ=255)
                        gAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCg                                            >  1:1760173/1‑70 (MQ=255)
                         atcgcCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGaa                                          <  1:2159385/69‑1 (MQ=255)
                           cgcCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACt                                        <  1:2348151/71‑1 (MQ=255)
                             ccTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGa                                      <  1:1048830/71‑1 (MQ=255)
                                tGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCg                                            <  1:2380276/62‑1 (MQ=255)
                                tGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCg                                            <  1:2393069/62‑1 (MQ=255)
                                      cGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGa                                                              <  1:2333306/38‑1 (MQ=255)
                                       gTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcg                            <  1:3826833/71‑1 (MQ=255)
                                       gTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcg                            <  1:3102805/71‑1 (MQ=255)
                                         cTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTc                                                          >  1:1225705/1‑39 (MQ=255)
                                         cTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTc                                                          >  1:1936633/1‑39 (MQ=255)
                                            gCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGaa                                          >  1:3005070/1‑52 (MQ=255)
                                              aGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCt                     <  1:2231111/71‑1 (MQ=255)
                                               gTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcg                            <  1:2218024/63‑1 (MQ=255)
                                               gTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTc                    <  1:538378/71‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAAt                                                 <  1:1073146/41‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAAt                                                 <  1:1308953/41‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAAt                                                 <  1:58672/41‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAAt                                                 <  1:2000966/41‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcgc                           <  1:2574275/63‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcgc                           <  1:1631962/63‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcg                            <  1:3812092/62‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcg                            <  1:2105408/62‑1 (MQ=255)
                                                tCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAgcg                            <  1:3650327/62‑1 (MQ=255)
                                                  ttAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGt                 >  1:683944/1‑71 (MQ=255)
                                                  ttAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAACGAACTGATCAAGCAGCGCATTGCTCAGt                 >  1:320740/1‑71 (MQ=255)
                                                    aaGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGaa                                          >  1:2861509/1‑44 (MQ=255)
                                                    aaGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGaa                                          >  1:2048493/1‑44 (MQ=255)
                                                    aaGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAg                  <  1:1350014/68‑1 (MQ=255)
                                                    aaGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGtt                >  1:289446/1‑70 (MQ=255)
                                                    aaGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGtt                >  1:2937775/1‑70 (MQ=255)
                                                    aaGGGGACTGGAGAATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGt                 <  1:665129/69‑1 (MQ=255)
                                                    aaGGGGACTGGAGAATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGt                 <  1:2527386/69‑1 (MQ=255)
                                                     aGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGa                                      >  1:1172835/1‑47 (MQ=255)
                                                     aGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTc                    >  1:2271053/1‑65 (MQ=255)
                                                     aGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTc               <  1:2294440/70‑1 (MQ=255)
                                                     aGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTc               <  1:1485970/70‑1 (MQ=255)
                                                     aGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTc               <  1:3014294/70‑1 (MQ=255)
                                                      ggggACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCaa             <  1:3360352/71‑1 (MQ=255)
                                                      ggggACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCaa             <  1:70316/71‑1 (MQ=255)
                                                         gACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAAt            >  1:2646188/1‑69 (MQ=255)
                                                           cTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAAtgttg        <  1:2377696/71‑1 (MQ=255)
                                                               aGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAATGTTgtgag    >  1:234490/1‑71 (MQ=255)
                                                                 cATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAATGTTgtgagtg  <  1:481781/71‑1 (MQ=255)
                                                                 cATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAATGTTgtgagtg  <  1:374873/71‑1 (MQ=255)
                                                                 cATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAATGTTgtgagtg  <  1:1919609/71‑1 (MQ=255)
                                                                 |                                                                      
GATGGCAGCGTACGCGGTCGTCTTGAGCGCCTTGCAGACGTCTTGCAGTCTTAAGGGGACTGGAGCATGCAACTGAATTCCACCGAAATCAGCGAACTGATCAAGCAGCGCATTGCTCAGTTCAATGTTGTGAGTG  >  W3110S.gb/3716760‑3716895

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