Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,249,096 C→G 5.4% P96A (CCC→GCC) ‡ thiF → thiamin (thiazole moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,249,0960CG 5.4% 303.2 / 5.1 93P96A (CCC→GCC) ‡thiFthiamin (thiazole moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base G (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.25e-01

CTCTTTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGC  >  Exported/2249027‑2249166
                                                                     |                                                                      
ctctTTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATc                                                                        <  1:108623/70‑1 (MQ=255)
    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:2776724/67‑1 (MQ=38)
              aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                             <  1:1262891/67‑1 (MQ=38)
              aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                             <  1:2955225/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:2066369/1‑71 (MQ=255)
                              aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                         <  1:2914453/71‑1 (MQ=255)
                              aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                         <  1:155787/71‑1 (MQ=255)
                               aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                        >  1:2568380/1‑71 (MQ=255)
                               aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                        >  1:2633378/1‑71 (MQ=255)
                               aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                        >  1:2329542/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: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:2972876/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)
                                          aGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGa                             <  1:1009428/71‑1 (MQ=255)
                                           gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                             >  1:2267382/1‑38 (MQ=255)
                                           gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                             >  1:1781061/1‑38 (MQ=255)
                                           gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                             >  1:2376599/1‑38 (MQ=255)
                                           gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                             >  1:3092991/1‑38 (MQ=255)
                                           gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                             >  1:3606443/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)
                                                  gCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaa                     <  1:991825/71‑1 (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:1282952/71‑1 (MQ=255)
                                                             aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt          <  1:1157989/71‑1 (MQ=255)
                                                               ttGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                        >  1:2437659/1‑55 (MQ=255)
                                                               ttGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                        >  1:2441946/1‑55 (MQ=255)
                                                                 gAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa       <  1:2431431/70‑1 (MQ=255)
                                                                 gAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa       <  1:1188539/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)
                                                                     cccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACggg   >  1:2650473/1‑70 (MQ=255)
                                                                     cccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGc  >  1:3343496/1‑71 (MQ=255)
                                                                     |                                                                      
CTCTTTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGC  >  Exported/2249027‑2249166

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: