Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,249,097 C→T 5.4% P96L (CCC→CTC) ‡ thiF → thiamin (thiazole moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,249,0970CT 5.4% 305.7 / 5.4 93P96L (CCC→CTC) ‡thiFthiamin (thiazole moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base T (2/3);  total (46/47)
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.06e-01

Only 100 of 102 total aligned reads displayed.

TTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCC  >  Exported/2249031‑2249167
                                                                  |                                                                      
ttACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGa                                                                     >  1:310457/1‑70 (MQ=255)
ttACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGa                                                                     >  1:3728011/1‑70 (MQ=255)
 tACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                   >  1:1931028/1‑71 (MQ=255)
 tACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                   >  1:1533246/1‑71 (MQ=255)
  aCCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                   <  1:3486575/70‑1 (MQ=255)
  aCCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACAt                                                                  >  1:2225614/1‑71 (MQ=255)
   ccACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                   <  1:377496/69‑1 (MQ=255)
        gAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCg                                                                      >  1:3006873/1‑61 (MQ=255)
         aaGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAAc                                                             <  1:2678762/69‑1 (MQ=255)
          aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                              <  1:2955225/67‑1 (MQ=38)
          aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                              <  1:2776724/67‑1 (MQ=38)
          aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                              <  1:1262891/67‑1 (MQ=38)
             tatCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                       <  1:2109898/71‑1 (MQ=255)
             tatCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                       <  1:3253787/71‑1 (MQ=255)
             tatCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCAGACATTCAACTGACAg                                                       <  1:3601901/71‑1 (MQ=255)
              atcgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                       >  1:2002520/1‑70 (MQ=255)
               tcgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCa                                                               <  1:1193030/61‑1 (MQ=255)
               tcgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCa                                                               <  1:1612525/61‑1 (MQ=255)
                cgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAt                                                    >  1:1006459/1‑71 (MQ=255)
                cgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCAGACATTCAACTGACAGCAt                                                    >  1:1439610/1‑71 (MQ=255)
                 gatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                       <  1:511637/67‑1 (MQ=255)
                 gatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAtt                                                   <  1:2838668/71‑1 (MQ=255)
                  atcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                       <  1:3313781/66‑1 (MQ=255)
                    cgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATccc                                                                       <  1:468279/48‑1 (MQ=255)
                    cgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAtt                                                   >  1:1704172/1‑68 (MQ=255)
                    cgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAt                                                    >  1:3391501/1‑67 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATccc                                                                       >  1:3596327/1‑43 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACAtt                                                                 >  1:1407061/1‑49 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTacaa                                               >  1:1053957/1‑67 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTacaa                                               >  1:2209311/1‑67 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACaacaac                                           >  1:2571391/1‑71 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACaacaac                                           >  1:59185/1‑71 (MQ=255)
                          aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                          <  1:2914453/71‑1 (MQ=255)
                          aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                          >  1:2066369/1‑71 (MQ=255)
                          aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                          <  1:155787/71‑1 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:2633378/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:3552400/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:3510594/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         <  1:2202934/71‑1 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:2329542/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:2568380/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:761964/1‑71 (MQ=255)
                            aTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCAACt                                                            >  1:471242/1‑51 (MQ=21)
                             tCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGtt                                       <  1:221539/71‑1 (MQ=255)
                             tCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGt                                        >  1:3226309/1‑70 (MQ=255)
                                cAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAAc                                    <  1:2550780/71‑1 (MQ=255)
                                cAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAAc                                    <  1:2637188/71‑1 (MQ=255)
                                 aGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACg                                   >  1:76835/1‑71 (MQ=255)
                                  ggTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACgg                                  >  1:3213508/1‑71 (MQ=255)
                                  ggTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACgg                                  >  1:1984309/1‑71 (MQ=255)
                                   gTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTAcaacaa                                            <  1:1062202/60‑1 (MQ=255)
                                   gTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTAcaacaa                                            <  1:2836354/60‑1 (MQ=255)
                                   gTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         <  1:1235658/63‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACTGTTAACGGGt                                <  1:1535744/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                <  1:1666646/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                <  1:3089535/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                <  1:259178/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                <  1:2972876/71‑1 (MQ=255)
                                      aGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGa                              <  1:1009428/71‑1 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                              >  1:3606443/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                              >  1:1781061/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                              >  1:2267382/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                              >  1:3092991/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                              >  1:2376599/1‑38 (MQ=255)
                                        ccAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGa                              >  1:746357/1‑69 (MQ=255)
                                             aGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGtt                        >  1:3680797/1‑70 (MQ=255)
                                             aGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                         >  1:2039817/1‑69 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                     <  1:1750771/71‑1 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                     <  1:849242/71‑1 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                     <  1:1236793/71‑1 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaa                      >  1:403044/1‑70 (MQ=255)
                                                gACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                     <  1:469085/70‑1 (MQ=255)
                                                 aCTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAg                   <  1:2266739/71‑1 (MQ=255)
                                                 aCTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAg                   <  1:1206511/71‑1 (MQ=255)
                                                   tGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                <  1:893385/56‑1 (MQ=255)
                                                   tGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                <  1:2833306/56‑1 (MQ=255)
                                                    gACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGtt                                       <  1:2544291/48‑1 (MQ=255)
                                                      cacaGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCg              <  1:3322120/71‑1 (MQ=255)
                                                      cacaGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCg              >  1:2780328/1‑71 (MQ=255)
                                                      cacaGTTGAATCCCGACAGTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCg              >  1:3370693/1‑71 (MQ=255)
                                                        caGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                         >  1:2991207/1‑42 (MQ=255)
                                                        caGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGt            >  1:3383977/1‑71 (MQ=255)
                                                        caGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGt            >  1:1798135/1‑71 (MQ=255)
                                                         aGTTGAATGTCGGGATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt           >  1:2863709/1‑71 (MQ=38)
                                                         aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAg                   >  1:3271254/1‑63 (MQ=255)
                                                         aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt           <  1:1157989/71‑1 (MQ=255)
                                                         aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt           <  1:1282952/71‑1 (MQ=255)
                                                           ttGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                         >  1:2437659/1‑55 (MQ=255)
                                                           ttGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                         >  1:2441946/1‑55 (MQ=255)
                                                             gAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa        <  1:1188539/70‑1 (MQ=255)
                                                             gAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa        <  1:2431431/70‑1 (MQ=255)
                                                              aaTCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCAc       <  1:671049/70‑1 (MQ=255)
                                                                tCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaaa                    <  1:2404064/55‑1 (MQ=255)
                                                                tCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCgg             <  1:2973394/62‑1 (MQ=255)
                                                                 cccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGc   >  1:3343496/1‑71 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt           >  1:1560025/1‑62 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt           >  1:374197/1‑62 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt           >  1:3087601/1‑62 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGcc  <  1:2241826/71‑1 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGcc  >  1:862905/1‑71 (MQ=255)
                                                                  |                                                                      
TTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCC  >  Exported/2249031‑2249167

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