Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,694,173 C→G 7.3% V15V (GTG→GTC)  iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,694,1730CG 7.3% 262.7 / 7.9 82V15V (GTG→GTC) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base G (4/2);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

GAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTC  >  Exported/1694105‑1694243
                                                                     |                                                                      
gAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCAc                                                                      <  1:4295352/71‑1 (MQ=255)
gAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCAc                                                                      <  1:167116/71‑1 (MQ=255)
       ggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTg                                                               >  1:987110/1‑71 (MQ=255)
         ttACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGt                                                             <  1:3312584/71‑1 (MQ=255)
           aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGtt                                                            >  1:557156/1‑70 (MQ=255)
           aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTg                                                           >  1:2765422/1‑71 (MQ=255)
            ccAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGc                                                          >  1:144185/1‑71 (MQ=255)
            ccAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGc                                                          >  1:2876430/1‑71 (MQ=255)
             cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         <  1:1483016/71‑1 (MQ=255)
             cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         <  1:693309/71‑1 (MQ=255)
             cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         <  1:4684166/71‑1 (MQ=255)
             cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         <  1:3036191/71‑1 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:2568256/1‑64 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         >  1:1885585/1‑69 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         >  1:2421011/1‑69 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         >  1:2262037/1‑69 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         >  1:3713322/1‑69 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         >  1:2843325/1‑69 (MQ=255)
               aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         >  1:1376285/1‑69 (MQ=255)
                       gTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCga                                               <  1:2664001/71‑1 (MQ=255)
                       gTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCga                                               <  1:949945/71‑1 (MQ=255)
                        tCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCcgcg                                                                 <  1:348113/53‑1 (MQ=17)
                        tCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCcgcg                                                                 <  1:1324200/53‑1 (MQ=17)
                          cATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTg                                                               >  1:3347546/1‑53 (MQ=17)
                          cATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTg                                                               >  1:3992115/1‑53 (MQ=17)
                          cATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTg                                                               >  1:4303097/1‑53 (MQ=17)
                          cATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCAGGTGGACCCGCGTg                                                               >  1:454780/1‑53 (MQ=18)
                          cATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGa                                                       >  1:3993151/1‑60 (MQ=255)
                          cATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGa                                                       >  1:2123433/1‑60 (MQ=255)
                           aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCg                                                                 <  1:1123912/49‑1 (MQ=255)
                           aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCg                                                                 <  1:2484044/49‑1 (MQ=255)
                           aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCg                                                                 <  1:2066527/49‑1 (MQ=255)
                            tCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATaa                                          <  1:1300003/71‑1 (MQ=255)
                            tCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATa                                           >  1:3073359/1‑70 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGt                                                                >  1:1956745/1‑47 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATaa                                          >  1:351449/1‑69 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        <  1:2167152/71‑1 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        <  1:1932629/71‑1 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        <  1:3000198/71‑1 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        <  1:4159018/71‑1 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        >  1:3092883/1‑71 (MQ=255)
                              gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTACACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        <  1:3825822/71‑1 (MQ=255)
                                 aTAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCgg                                     >  1:1994963/1‑71 (MQ=255)
                                 aTAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCgg                                     >  1:4112090/1‑71 (MQ=255)
                                  tAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGt                                    <  1:3757865/71‑1 (MQ=255)
                                  tAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGt                                    <  1:3397106/71‑1 (MQ=255)
                                  tAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGt                                    <  1:1217204/71‑1 (MQ=255)
                                   aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCa                                                                       >  1:1523454/1‑35 (MQ=255)
                                   aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:1478470/1‑44 (MQ=255)
                                   aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:2911189/1‑44 (MQ=255)
                                   aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:3759844/1‑44 (MQ=255)
                                   aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:391077/1‑44 (MQ=255)
                                    aaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:438666/1‑43 (MQ=255)
                                    aaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                              >  1:1411993/1‑43 (MQ=255)
                                    aaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTaa                                  <  1:3434808/71‑1 (MQ=255)
                                      cTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGt                                                             >  1:694269/1‑42 (MQ=255)
                                      cTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGt                                                             >  1:1485243/1‑42 (MQ=255)
                                        gCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCg                                                         <  1:787599/44‑1 (MQ=255)
                                        gCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                        <  1:3041108/61‑1 (MQ=255)
                                         catcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:4165728/69‑1 (MQ=255)
                                         catcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             <  1:2188294/71‑1 (MQ=255)
                                             atTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGc                         <  1:1935625/71‑1 (MQ=255)
                                                ttCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGagtag                                                   >  1:3720751/1‑42 (MQ=255)
                                                ttCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                       >  1:3169639/1‑70 (MQ=255)
                                                  ctcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCta                     <  1:772784/70‑1 (MQ=255)
                                                     ggCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGc                         <  1:1556782/63‑1 (MQ=255)
                                                     ggCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCtata                   <  1:1805488/69‑1 (MQ=255)
                                                     ggCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaaa                 >  1:2502021/1‑71 (MQ=255)
                                                        aa‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTg                          >  1:989232/1‑59 (MQ=255)
                                                        aa‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTc              >  1:3567120/1‑71 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:1676248/50‑1 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:4226358/50‑1 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:2022115/50‑1 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:582013/50‑1 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:3073211/50‑1 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                               <  1:3297253/50‑1 (MQ=255)
                                                             cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             >  1:2303063/1‑52 (MQ=255)
                                                              gcgGGTCCACCGGCGTTGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             <  1:1572957/51‑1 (MQ=255)
                                                              gcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             <  1:526466/51‑1 (MQ=255)
                                                              gcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             <  1:1489171/51‑1 (MQ=255)
                                                              gcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             <  1:2615394/51‑1 (MQ=255)
                                                              gcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                             <  1:3392036/51‑1 (MQ=255)
                                                                     cACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc  <  1:4242561/71‑1 (MQ=255)
                                                                     cACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc  <  1:2008304/71‑1 (MQ=255)
                                                                     cACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc  <  1:2999959/71‑1 (MQ=255)
                                                                     |                                                                      
GAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTC  >  Exported/1694105‑1694243

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: ‑