Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,061,569 G→A 6.6% A26T (GCA→ACA) ‡ rpsU → 30S ribosomal subunit protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,061,5690GA 6.6% 244.8 / 5.2 75A26T (GCA→ACA) ‡rpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base G (40/30);  new base A (3/2);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGC  >  Exported/2061504‑2061634
                                                                 |                                                                 
ttAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAgg                                                               <  1:1864572/70‑1 (MQ=255)
ttAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgt                                                              >  1:711071/1‑71 (MQ=255)
 tAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtg                                                             >  1:1850508/1‑71 (MQ=255)
 tAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtg                                                             >  1:307214/1‑71 (MQ=255)
  aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            >  1:364969/1‑71 (MQ=255)
  aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            >  1:4416356/1‑71 (MQ=255)
  aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            >  1:3040003/1‑71 (MQ=255)
   aaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            <  1:2317958/70‑1 (MQ=255)
        cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc                                                                  >  1:127102/1‑59 (MQ=255)
        cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc                                                                  >  1:2404828/1‑59 (MQ=255)
        cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTgg                                                       <  1:1992416/70‑1 (MQ=255)
          tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            <  1:2518852/63‑1 (MQ=255)
          tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            <  1:1463528/63‑1 (MQ=255)
          tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            <  1:880993/63‑1 (MQ=255)
           gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCgg                                                    >  1:2811013/1‑70 (MQ=255)
            aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                            <  1:2022568/61‑1 (MQ=255)
            aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCgg                                                    >  1:285253/1‑69 (MQ=255)
            aaaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaa                                                  >  1:1019933/1‑71 (MQ=255)
            aaaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaa                                                  >  1:1449614/1‑71 (MQ=255)
            aaaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGa                                                   >  1:432435/1‑70 (MQ=255)
             aaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaac                                                 <  1:3206717/71‑2 (MQ=255)
             aaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaac                                                 <  1:1410787/71‑2 (MQ=255)
              aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                                <  1:987020/71‑1 (MQ=255)
              aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                                >  1:283093/1‑71 (MQ=255)
              aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                                <  1:4660032/71‑1 (MQ=255)
              aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                                <  1:4741940/71‑1 (MQ=255)
                      gTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:1479754/71‑1 (MQ=255)
                      gTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:1336043/71‑1 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCTTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                      >  1:4333830/1‑70 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTc                                                          >  1:375935/1‑50 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                                  >  1:544467/1‑58 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                                  >  1:353370/1‑58 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                                  >  1:946537/1‑58 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                                  >  1:2890048/1‑58 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                      >  1:2437457/1‑70 (MQ=255)
                         cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                      <  1:3036209/70‑1 (MQ=255)
                          gACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTg                                     <  1:4286524/70‑1 (MQ=255)
                              tAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGtt                                                           <  1:2920429/44‑1 (MQ=255)
                               aGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGt                                  >  1:3257521/1‑68 (MQ=255)
                               aGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTc                                >  1:1869581/1‑70 (MQ=255)
                                gCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCg                                                     <  1:3799372/48‑1 (MQ=255)
                                gCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTa                              <  1:1819439/71‑1 (MQ=255)
                                gCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGAGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCt                               <  1:4684452/70‑1 (MQ=255)
                                  tctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTAt                             <  1:718194/70‑1 (MQ=255)
                                   ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:1663780/58‑1 (MQ=255)
                                   ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:1307541/58‑1 (MQ=255)
                                   ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:867195/58‑1 (MQ=255)
                                   ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:2316137/58‑1 (MQ=255)
                                   ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                                        <  1:173323/58‑1 (MQ=255)
                                      cGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa                         >  1:2984343/1‑70 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                      >  1:2857210/1‑53 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTTGTGAGTTCTATGAAAAACCg                    <  1:4737592/71‑1 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa                         >  1:3080315/1‑66 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa                         >  1:735017/1‑66 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa                         >  1:4228349/1‑66 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa                         >  1:1500118/1‑66 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa                         >  1:4551479/1‑66 (MQ=255)
                                          gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGCAAAAcc                     >  1:102703/1‑70 (MQ=255)
                                            ttCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAcc                     <  1:1057783/68‑1 (MQ=255)
                                              cAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa                >  1:1024470/1‑71 (MQ=255)
                                              cAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa                >  1:422836/1‑71 (MQ=255)
                                                aGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                      >  1:4381512/1‑63 (MQ=255)
                                                aGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAcc              >  1:1270394/1‑71 (MQ=255)
                                                 gCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCg             <  1:926129/71‑1 (MQ=255)
                                                      ccTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTc                                >  1:2343183/1‑47 (MQ=255)
                                                      ccTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGa                           >  1:4254704/1‑52 (MQ=255)
                                                      ccTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACg         >  1:1597642/1‑70 (MQ=255)
                                                      ccTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACg         >  1:3127679/1‑70 (MQ=255)
                                                      ccTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACg         >  1:2985056/1‑70 (MQ=255)
                                                       cTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                      <  1:3887087/40‑1 (MQ=255)
                                                       cTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGt        <  1:1482646/70‑1 (MQ=255)
                                                       cTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGt        <  1:1339688/70‑1 (MQ=255)
                                                        tGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGGTCTATGAAAAACCGACTACCGAACGt        >  1:1083863/1‑69 (MQ=255)
                                                         gCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAg     >  1:4242072/1‑71 (MQ=255)
                                                            aaaaaGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAgcgc  >  1:2868685/1‑71 (MQ=255)
                                                                 |                                                                 
TTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGC  >  Exported/2061504‑2061634

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