Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,976,997 G→A 6.7% L83L (TTG→TTA sprT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,976,9970GA 6.7% 241.0 / 5.5 75L83L (TTG→TTAsprTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (40/30);  new base A (2/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGC  >  Exported/1976928‑1977067
                                                                     |                                                                      
tCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGc                                                                       >  1:341530/1‑71 (MQ=255)
         gCTGTTGGAAAACAGTGAAGCTTTTATTTAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATgg                                                               >  1:192099/1‑70 (MQ=255)
          cTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGa                                                              <  1:823497/70‑1 (MQ=255)
          cTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGa                                                              <  1:793841/70‑1 (MQ=255)
           tGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGaaa                                                            >  1:3529479/1‑71 (MQ=255)
            gTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGaaa                                                            >  1:628347/1‑70 (MQ=255)
             ttGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGt                                                                   <  1:2402843/62‑1 (MQ=255)
               ggAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACAtt                                                        >  1:1752393/1‑71 (MQ=255)
               ggAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACAt                                                         <  1:961574/70‑1 (MQ=255)
                 aaaaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACAtt                                                        <  1:3114092/69‑1 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGc                                                                       >  1:2069925/1‑52 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGc                                                                       >  1:1761578/1‑52 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGc                                                                       >  1:530829/1‑52 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGc                                                                       >  1:3550052/1‑52 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAAc                                                           >  1:2204580/1‑64 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCgg                                                    >  1:696517/1‑71 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCgg                                                    >  1:86155/1‑71 (MQ=255)
                   aaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCgg                                                    >  1:1183451/1‑71 (MQ=255)
                    aCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGc                                                   >  1:2082764/1‑71 (MQ=255)
                    aCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGc                                                   >  1:2705401/1‑71 (MQ=255)
                        tGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCcgcg                                               >  1:3550337/1‑71 (MQ=255)
                        tGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACAGTTGCTGGTATGGAAACATTTCGGCcgcg                                               >  1:2378894/1‑71 (MQ=255)
                          aaGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCcgcg                                               >  1:889215/1‑69 (MQ=255)
                          aaGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTACTGGTATGGAAACATTTCGGCCGCGt                                              <  1:2741821/70‑1 (MQ=255)
                           aGCTTTTATTGTAGAAGTGGTACCGCACGAACTGGCTCATTTACTGGTAt                                                                 <  1:1520907/50‑1 (MQ=255)
                            gCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCcgc                                                >  1:2048602/1‑66 (MQ=255)
                            gCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAgc                                           <  1:1211778/71‑1 (MQ=255)
                             cTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTgg                                                                    <  1:453098/45‑1 (MQ=255)
                             cTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTgg                                                                    <  1:3306780/45‑1 (MQ=255)
                             cTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAgcg                                          <  1:2797147/71‑1 (MQ=255)
                              ttttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAgcgc                                         >  1:830225/1‑71 (MQ=255)
                                   ttGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCcaca                                     <  1:791382/70‑1 (MQ=255)
                                    tGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACAt                                    <  1:493973/70‑1 (MQ=255)
                                       agaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGc                                 >  1:1285206/1‑70 (MQ=255)
                                       agaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGc                                 >  1:1044504/1‑70 (MQ=255)
                                       agaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGc                                 >  1:3388384/1‑70 (MQ=255)
                                       agaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGc                                 >  1:132761/1‑70 (MQ=255)
                                       agaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGc                                 >  1:1527548/1‑70 (MQ=255)
                                        gaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCaa                               <  1:2927198/71‑1 (MQ=255)
                                         aagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCaaa                              <  1:451684/71‑1 (MQ=255)
                                           gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGt                                              <  1:1013744/53‑1 (MQ=255)
                                           gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGt                                              <  1:996270/53‑1 (MQ=255)
                                           gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAaga                            <  1:2579457/71‑1 (MQ=255)
                                           gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAaga                            <  1:1509750/71‑1 (MQ=255)
                                            tGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAgc                                           <  1:1928924/55‑1 (MQ=255)
                                              gTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtg                         >  1:1722643/1‑71 (MQ=255)
                                              gTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagt                          <  1:2905113/70‑1 (MQ=255)
                                               tACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgg                        <  1:2710432/71‑1 (MQ=255)
                                               tACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgg                        <  1:2303808/71‑1 (MQ=255)
                                               tACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgg                        <  1:829249/71‑1 (MQ=255)
                                               tACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtg                         >  1:953264/1‑70 (MQ=255)
                                                 ccGCACGAACTGGCACATTTACTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaa                      <  1:709858/71‑1 (MQ=255)
                                                  cGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaag                     <  1:2868789/71‑1 (MQ=255)
                                                  cGCACGAACTGGCACATTTACTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaag                     >  1:1750774/1‑71 (MQ=255)
                                                   gCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAgcg                                          >  1:1441273/1‑49 (MQ=255)
                                                   gCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagt                    >  1:270088/1‑71 (MQ=255)
                                                   gCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagt                    <  1:2135255/71‑1 (MQ=255)
                                                     aCGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAgcg                                          >  1:2108471/1‑47 (MQ=255)
                                                     aCGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagtg                   >  1:99183/1‑70 (MQ=255)
                                                     aCGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagtg                   >  1:751459/1‑70 (MQ=255)
                                                      cGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtg                         <  1:1433358/63‑1 (MQ=255)
                                                      cGAACTGGCACATTTACTGGTATGGAAACATTTCg                                                     >  1:1420592/1‑35 (MQ=255)
                                                       gAACTGGCACATTTGCTGGTATGTAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtg                         <  1:1392595/62‑1 (MQ=255)
                                                       gAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgat                >  1:1136555/1‑71 (MQ=255)
                                                        aaCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaa                      <  1:2499333/64‑1 (MQ=255)
                                                        aaCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgat                >  1:3139263/1‑70 (MQ=255)
                                                         aCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgatga              <  1:2248028/71‑1 (MQ=255)
                                                              cacaTTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGa          >  1:2567263/1‑70 (MQ=255)
                                                               acaTTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgg                        >  1:3619149/1‑55 (MQ=255)
                                                                  tttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACg     >  1:1441333/1‑71 (MQ=255)
                                                                  tttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACg     <  1:1952809/71‑1 (MQ=255)
                                                                  tttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACg     <  1:1351010/71‑1 (MQ=255)
                                                                   ttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGt    >  1:1288551/1‑71 (MQ=255)
                                                                    tGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTg   >  1:547075/1‑71 (MQ=255)
                                                                    tGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTg   >  1:92536/1‑71 (MQ=255)
                                                                     gCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGc  >  1:1499473/1‑71 (MQ=255)
                                                                     |                                                                      
TCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGC  >  Exported/1976928‑1977067

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: