Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 690,095 A→G 47.1% D403D (GAT→GAC lnt ← apolipoprotein N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb690,0950AG47.1% ‑1.8 / 46.7 34D403D (GAT→GAClntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (6/12);  new base G (8/8);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTC  >  W3110S.gb/690034‑690163
                                                             |                                                                    
gCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTGTCCGGGCg                                                              >  1:1350057/1‑70 (MQ=255)
      gACCAATAGATTTACGAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTa                                                       >  1:3539825/1‑71 (MQ=255)
      gACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTa                                                       >  1:3631908/1‑71 (MQ=255)
      gACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCAGGGCGGAAGTTa                                                       >  1:3143508/1‑71 (MQ=255)
           aTAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATcgcg                                                  <  1:375419/71‑1 (MQ=255)
           aTAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATcgcg                                                  <  1:21484/71‑1 (MQ=255)
            tAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCgcgc                                                 <  1:2426452/71‑1 (MQ=255)
               aTTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCgcgc                                                 <  1:2235084/68‑1 (MQ=255)
                 ttACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCAc                                               >  1:3117198/1‑68 (MQ=255)
                 ttACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCAc                                               >  1:3857927/1‑68 (MQ=255)
                              cATCGTTGGAGATAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATg                               <  1:2579453/71‑1 (MQ=255)
                               aTCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATg                               <  1:851623/70‑1 (MQ=255)
                                 cGTTGGAGATACTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATg                               <  1:1883575/68‑1 (MQ=255)
                                   ttGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtctc                          <  1:3856738/71‑1 (MQ=255)
                                        gaTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAt                             >  1:2464141/1‑63 (MQ=255)
                                        gaTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAt                             >  1:2850470/1‑63 (MQ=255)
                                        gaTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAt                             >  1:505880/1‑63 (MQ=255)
                                        gaTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAt                             >  1:3939987/1‑63 (MQ=255)
                                         aTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCa                                                <  1:2164113/43‑1 (MQ=255)
                                         aTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCa                                                <  1:476015/43‑1 (MQ=255)
                                         aTAGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCa                    <  1:2389998/71‑1 (MQ=255)
                                         aTAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGc                     >  1:3208880/1‑70 (MQ=255)
                                         aTAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCa                    <  1:3570688/71‑1 (MQ=255)
                                         aTAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCa                    <  1:4084558/71‑1 (MQ=255)
                                          tAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCa                    <  1:1935991/70‑1 (MQ=255)
                                           aGTCAGCAGATAGTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCaaa                  >  1:3826225/1‑71 (MQ=255)
                                             tCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCa                    <  1:222491/67‑1 (MQ=255)
                                               agcagATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGc              >  1:2191390/1‑71 (MQ=255)
                                               agcagATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGc              >  1:1611137/1‑71 (MQ=255)
                                                       gTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTa       <  1:1295379/70‑1 (MQ=255)
                                                       gTCGGTGTCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTa       <  1:2665795/70‑1 (MQ=255)
                                                       gTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATg                               <  1:611010/46‑1 (MQ=255)
                                                       gTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTa       <  1:1476868/70‑1 (MQ=255)
                                                           gTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTc  <  1:3752148/71‑1 (MQ=255)
                                                             |                                                                    
GCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTC  >  W3110S.gb/690034‑690163

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: