Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,746,079 C→A 11.5% G166G (GGC→GGA srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,746,0790CA11.5% 210.4 / 9.9 69G166G (GGC→GGAsrmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (3/5);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGC  >  Exported/1746015‑1746147
                                                                |                                                                    
aCTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGc                                                                 <  1:1265895/70‑1 (MQ=255)
   tCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAg                                                             >  1:3588607/1‑71 (MQ=255)
    cGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCa                                                              <  1:4779973/69‑1 (MQ=255)
        tGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGatat                                                        >  1:2686935/1‑71 (MQ=255)
        tGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGatat                                                        >  1:1042860/1‑71 (MQ=255)
         gCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATc                                                       >  1:674286/1‑71 (MQ=255)
            gcgcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCtt                                                                    <  1:5432963/55‑1 (MQ=255)
            gcgcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCtt                                                                    <  1:2949589/55‑1 (MQ=255)
             cgcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGa                                                     >  1:113931/1‑69 (MQ=255)
              gcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACa                                                  <  1:3492741/71‑1 (MQ=255)
               cgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACa                                                  >  1:4811583/1‑70 (MQ=255)
                  ttGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAAc                                                   <  1:5374915/66‑1 (MQ=255)
                  ttGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAAc                                                   <  1:4353710/66‑1 (MQ=255)
                   tGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATAtt                                              >  1:2074983/1‑70 (MQ=255)
                     aaaCGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCt                                           <  1:1608963/71‑1 (MQ=255)
                         gCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTgg                                         >  1:5661286/1‑69 (MQ=255)
                         gCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTgg                                         >  1:1894554/1‑69 (MQ=255)
                         gCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTgg                                         >  1:2237582/1‑69 (MQ=255)
                            gATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGaaa                                    >  1:4574408/1‑71 (MQ=255)
                             aTCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGaa                                     <  1:2200044/69‑1 (MQ=255)
                              tCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAcg                                  >  1:13016/1‑71 (MQ=255)
                              tCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAc                                   >  1:4128484/1‑70 (MQ=255)
                               ccTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACa                                                  <  1:4507421/54‑1 (MQ=255)
                                   gacgaAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCt                              <  1:1464438/70‑1 (MQ=255)
                                    atgaggCCGACCGCATGCTGGATATGGGATTtagcgatgccattgatgat                                                 <  1:3035133‑M2/48‑20 (MQ=255)
                                    acgaAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTgg                            >  1:3941/1‑71 (MQ=255)
                                        aGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTa                        <  1:2352681/71‑1 (MQ=255)
                                        aGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTa                        <  1:820739/71‑1 (MQ=255)
                                         gCCGACCGCATGCTGGATATGGGATTtagcgatgccattgatgatg                                                >  1:5502308‑M2/1‑26 (MQ=255)
                                         gCCGACCGCATGCTGGATATGGGATTtagcgatgccattgatg                                                   >  1:604080‑M2/1‑26 (MQ=255)
                                         gCCGACCGCATGCTGGATATGGGATTtagcgatgccattgatg                                                   >  1:953336‑M2/1‑26 (MQ=255)
                                         gCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCt                                           >  1:2356918/1‑51 (MQ=255)
                                          cAGACCGTATGCTGGATATGGGCGTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCt                              <  1:4293497/63‑1 (MQ=255)
                                            gACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAc                                   >  1:2222886/1‑56 (MQ=255)
                                            gACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAc                                   >  1:397465/1‑56 (MQ=255)
                                             aCCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTaaa                      >  1:5340464/1‑68 (MQ=255)
                                             aCCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTaaa                      >  1:513643/1‑68 (MQ=255)
                                             aCCGCATGCTGGATATGGGATTtagcgatgccattgatgat                                                 <  1:1296963‑M2/41‑20 (MQ=255)
                                             aCCGCATGCTGGATATGGGATTtagcgatgccattgatgat                                                 <  1:1430596‑M2/41‑20 (MQ=255)
                                             aCCGCATGCTGGATATGGGATTtagcgatgccattgatga                                                  <  1:4334062‑M2/40‑19 (MQ=255)
                                              ccGCATGCTGGATATGGGATTtagcgatgccattgatga                                                  <  1:2077166‑M2/39‑19 (MQ=255)
                                               cGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGAc                 <  1:338415/71‑1 (MQ=255)
                                               cGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTCGCGAAACGCGCTGGCGTaaa                      <  1:5158084/66‑1 (MQ=255)
                                                 tATGCTGGATATGGTCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGAccc               >  1:1292817/1‑71 (MQ=255)
                                                 tATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGAccc               >  1:2181596/1‑71 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTgg                                         >  1:5094537/1‑43 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTaaa                      <  1:1765459/62‑1 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTa                        <  1:829784/60‑1 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAg                   >  1:1493264/1‑65 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             <  1:3271073/71‑1 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             >  1:2495703/1‑71 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             <  1:224483/71‑1 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             <  1:2347464/71‑1 (MQ=255)
                                                   tGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             <  1:5646632/71‑1 (MQ=255)
                                                    gCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACgcgc                               <  1:1273878/52‑1 (MQ=255)
                                                     cTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGa                                      >  1:95597/1‑44 (MQ=255)
                                                     cTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGa                                      >  1:4886729/1‑44 (MQ=255)
                                                     cTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTaaa                      >  1:1606350/1‑60 (MQ=255)
                                                     cTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             >  1:3303003/1‑69 (MQ=255)
                                                     cTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg             >  1:3228107/1‑69 (MQ=255)
                                                        gATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCtt        <  1:267169/71‑1 (MQ=255)
                                                        gATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCtt        <  1:888634/71‑1 (MQ=255)
                                                          tatGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCtttt      >  1:2138441/1‑71 (MQ=255)
                                                              ggCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGc  <  1:2696519/71‑1 (MQ=255)
                                                               gCTTCGCTCAGGATATCTAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGc  <  1:4076971/70‑1 (MQ=255)
                                                               gCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAg                   <  1:4566419/53‑1 (MQ=255)
                                                               gCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGAc                 <  1:5438494/55‑1 (MQ=255)
                                                               gCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGAc                 <  1:342319/55‑1 (MQ=255)
                                                               gCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGAc                 <  1:2021068/55‑1 (MQ=255)
                                                                |                                                                    
ACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGC  >  Exported/1746015‑1746147

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: