Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,746,064 T→C 14.3% R161R (CGT→CGC srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,746,0640TC14.3% 129.3 / 13.2 56R161R (CGT→CGCsrmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base C (3/5);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

ACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCC  >  Exported/1746000‑1746134
                                                                |                                                                      
aCATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTg                                                                  >  1:342939/1‑71 (MQ=255)
  aTAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATTCTCGACGAAGCAGACCGTATGCTGGa                                                                <  1:5552704/71‑1 (MQ=255)
   tAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGa                                                                <  1:4745681/70‑1 (MQ=255)
     aaaGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGatat                                                             <  1:5657354/71‑1 (MQ=255)
     aaaGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGatat                                                             >  1:390100/1‑71 (MQ=255)
      aagaagAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATg                                                            <  1:2586786/71‑1 (MQ=255)
        gaagaGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTa                                                                       >  1:2326723/1‑58 (MQ=255)
               aCTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGc                                                    <  1:1265895/70‑1 (MQ=255)
                gttCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGata                                                              >  1:5682657/2‑59 (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)
                        gCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTa                                                                       >  1:5675419/1‑42 (MQ=255)
                        gCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATc                                          >  1:674286/1‑71 (MQ=255)
                           gcgcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCtt                                                       <  1:5432963/55‑1 (MQ=255)
                           gcgcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCtt                                                       <  1:2949589/55‑1 (MQ=255)
                            cgcgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATg                                                                     >  1:4674973/1‑40 (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)
                                  tGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATgg                                                           >  1:2401611/1‑44 (MQ=255)
                                  tGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATAtt                                 >  1:2074983/1‑70 (MQ=255)
                                    aaaCGCTGATCCTCGACGAAGCAGACCGTATGCTTGata                                                              <  1:931081/39‑1 (MQ=255)
                                    aaaCGCTGATCCTCGACGAAGCAGACCGTATGCTGGata                                                              <  1:3608051/39‑1 (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:820739/71‑1 (MQ=255)
                                                       aGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTa           <  1:2352681/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:397465/1‑56 (MQ=255)
                                                           gACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAc                      >  1:2222886/1‑56 (MQ=255)
                                                            aCCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTaaa         >  1:5340464/1‑68 (MQ=255)
                                                            aCCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTaaa         >  1:513643/1‑68 (MQ=255)
                                                            aCCGCATGCTGGATATGGGATTtagcgatgccattgatgat                                    <  1:1430596‑M2/41‑20 (MQ=255)
                                                            aCCGCATGCTGGATATGGGATTtagcgatgccattgatgat                                    <  1:1296963‑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)
                                                                |                                                                      
ACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCC  >  Exported/1746000‑1746134

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