Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 983,805 C→G 27.2% G31G (GGC→GGG msbA → fused lipid transporter subunits of ABC superfamily: membrane component/ATP‑binding component

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606983,8050C→G27.2% 2.4 ‑3.3 0.01893msbAfused lipid transporter subunits of ABC superfamily: membrane component/ATP‑binding component
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (5/8):  total (10/8)
Fisher's exact test strand distribution p-value = 3.59e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGCGATACCT  >  REL606/983770‑983839
                                   |                                  
aTTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggc                                    <  1:5319253/36‑1 (MQ=255)
   gcgcCTTTCAAAGCGGGTCTGATCGTGGCGGGGGtg                                 >  1:4069479/1‑35 (MQ=255)
    cgcCTTTCAAAGCGGGTCTGACCGTGGCGGGCGTAg                                <  1:3211561/36‑1 (MQ=255)
     gcCTTTCAAAGCGGGTCTGATCGTGGCGGGGGttgg                               >  1:867432/1‑33 (MQ=38)
     gcCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGc                               <  1:5670593/36‑1 (MQ=255)
        tttCAAAGCGGGCCGGATCGTGGCGGGGGTGGCGtt                            >  1:5741420/1‑36 (MQ=19)
         ttCAAAGCGGGTATGATCGTGGCGGGCGTAGCGTTa                           <  1:4191518/36‑1 (MQ=255)
            aaaGCGGGTCTGATCGTGGCGGGGGGAGGGTTTATc                        >  1:2893955/1‑36 (MQ=19)
              aGCGGGTCTGATCGTGGTGGGCGTAGCGTTAATcat                      >  1:3040706/1‑34 (MQ=37)
              aGCGGGTCTGATCGTGGCGGGCGGTGCGGTAATTCt                      >  1:5024340/1‑36 (MQ=16)
                cGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCa                    <  1:3833836/36‑1 (MQ=255)
                 gggTCTGATCGTGGCGGGGGGGGCGTTTAGCCTcca                   >  1:5180415/1‑34 (MQ=12)
                 gggTCTGATCGTGGCGGGCGTGGCGTTTATCCTcca                   >  1:1561875/1‑34 (MQ=21)
                  ggTCTGATCGTGGCGGGCGTAGCGTTAATCCTCCAc                  >  1:962467/1‑36 (MQ=255)
                  ggTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAAc                  <  1:1434928/36‑1 (MQ=255)
                    tCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGc                <  1:862292/36‑1 (MQ=255)
                     cTGATCGTGGCGGGCGTCGCGTTAATCCTCAACGca               <  1:1477100/36‑1 (MQ=255)
                          cGTGGCGGGCGTAGGGTTTATCCTCAAAGcagccag          >  1:5250567/1‑36 (MQ=37)
                           gTGGCGGGCGTAGCGTTAATCCTCAACGcagccagc         <  1:3989869/36‑1 (MQ=255)
                              gcgggcgTAGCGTTAATCCTCAACGCAGCCAGCGAt      <  1:672832/36‑1 (MQ=255)
                              gcgggcgTAGCGTTAATCCTCAACGCAGCCAGCGAt      <  1:1766005/36‑1 (MQ=255)
                               cgggcgTAGCGTTAATCCTCAACGCAGCCAGCGATa     >  1:213256/1‑36 (MQ=255)
                                gggcgTAGCGTTAATCCTCAACGCAGCCAGCGATAc    >  1:3177227/1‑36 (MQ=255)
                                gggcgTAGCGTTAATCCTCAACGCAGCCAGCGATAc    <  1:2578584/36‑1 (MQ=255)
                                 ggcgTAGCGTTAATCCTCAACGCAGCCAGCGATAcc   <  1:2376790/36‑1 (MQ=255)
                                 ggcgTAGCGTTAATCCTCAACGCAGCCAGCGATAcc   <  1:1802061/36‑1 (MQ=255)
                                 ggcgTAGCGTTAATCCTCAACGCAGCCAGCGATAcc   >  1:1295048/1‑36 (MQ=255)
                                  gcgTAGCGTTAATCCTCAACGCAGCCAGCGATACCt  >  1:5095146/1‑36 (MQ=255)
                                   |                                  
ATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGCGATACCT  >  REL606/983770‑983839

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 25 ≤  ATCG  < 37 ≤  ATCG