Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,775,878 Δ1 bp 11.1% intergenic (‑48/+201) insJ‑3 ← / ← cysH IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,775,8780G→.11.1% 18.5 ‑2.6 ‑0.154intergenic (‑48/+201)insJ‑3/cysHIS150 hypothetical protein/phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (20/28):  total (20/34)
Fisher's exact test strand distribution p-value = 7.40e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.91e-01

CCATTTCATCGTCCAACAAAATGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCTTTTGTCCAT  >  REL606/2775844‑2775912
                                   |                                    
ccATTTCATCGTCCAACAAAATGGGTGCAGTAC.Agg                                     <  1:2917255/36‑1 (MQ=32)
ccATTTCATCGTCCAACAAAATGGGTGCAGTAC.Agg                                     <  1:2686804/36‑1 (MQ=32)
ccATTTCATCGTCCAACAAAATGGGTGCAGTAC.Agg                                     >  1:6351287/1‑36 (MQ=32)
 cATTTCATCGTCCAACAAAATGGGTGCAGTAC.AGGt                                    <  1:7415775/36‑1 (MQ=35)
 cATTTCATCGTCCAACAAAATGGGTGCAGTAC.AGGt                                    >  1:5461320/1‑36 (MQ=35)
  aTTTCATCGTCCAACAAAATGGGTGCAGTAC.AGGtt                                   <  1:6301714/36‑1 (MQ=38)
   tttCATCGTCCAACAAAATGGGTGCAGTAC.AGGTTg                                  <  1:270052/36‑1 (MQ=40)
    ttCATAGTCCAAAAAAATGGGTGCAGTAC.AGGTTGa                                 >  1:1414446/1‑36 (MQ=17)
     tCATCGTCCAACAAAATGGGTGCAGTAC.AGGTTGAt                                >  1:2184627/1‑36 (MQ=255)
      cATCGTCCAACAAAATGGGTGCAGTAC.AGGTTGATa                               <  1:3695949/36‑1 (MQ=255)
       aTCGTCCAACAAAATGGGTGCAGTAC.AGGTTGATAg                              >  1:603805/1‑36 (MQ=255)
       aTCGTCCAACAAAATGGGTGCAGTAC.AGGTTGATAg                              <  1:3478375/36‑1 (MQ=255)
       aTCGTCCAACAAAATGGGTGCAGTAC.AGGTTGATAg                              >  1:2218754/1‑36 (MQ=255)
       aTCGTCCAACAAAATGGGTGCAGTAC.A.GTTGATAG..a                           <  1:1793419/36‑1 (MQ=11)
        tCGTCCAACAAAATGGGTGCAGTAC.A.GTTGATAG..aa                          <  1:6341575/36‑1 (MQ=14)
        tCGTCCAACAAAATGGGTGCAGTAC.A.GTTGATAG..aa                          <  1:3062395/36‑1 (MQ=14)
        tCGTCCAACAAAATGGGTGCAGTAC.A.GTTGATAG..aa                          <  1:7841047/36‑1 (MQ=14)
         cGTCCAACAAAATGGGTGCAGTAC.AGGTTGATAG..aa                          >  1:7235741/1‑36 (MQ=255)
          gTCCAACAAAATGGGTGCAGTAC.AGGTTGATAG..aaa                         <  1:5364957/36‑1 (MQ=255)
          gTCCAACAAAATGGGTGCAGTAC.AGGTTGATAG..aaa                         <  1:1864805/36‑1 (MQ=255)
          gTCCAACAAAATGGGTGCAGTAC.AGGTTGATAG..aaa                         <  1:1092837/36‑1 (MQ=255)
            ccAACAAAATGGGTGCAGTAC.AGGTTGATAG..AAata                       <  1:4531070/36‑1 (MQ=255)
            ccAACAAAATGAGTGCAGTAC.AGGTTGATAG..AAata                       <  1:3869308/36‑1 (MQ=255)
             caacaaAATGGGTGCAGTAC.AGGTTGATAG..TAaaat                      <  1:147504/36‑1 (MQ=18)
              aacaaAATGGGTGCAGTAC.AGGTTGATAG..AAATATc                     >  1:5482093/1‑36 (MQ=255)
              aacaaAATGGGTGCAGTAC.AGGTTGATAG..AAATATc                     >  1:1755429/1‑36 (MQ=255)
              aacaaAATGGGTGCAGTAC.AGGTTGATAG..AAATATc                     <  1:3594722/36‑1 (MQ=255)
              aacaaAATGGGTGCAGTAC.A.GTTGATAG..AAATATCg                    <  1:410732/36‑1 (MQ=18)
               acaaAATGGGTGCAGTAC.AGGTTGATAG..AAATATCg                    <  1:1872656/36‑1 (MQ=255)
                caaAATGGGTGCAGTAC.AGGTTGATAG..AAATATCGc                   <  1:7805265/36‑1 (MQ=255)
                caaAATGGGTGCAGTAC.AGGTTGATAG..AAATATCGc                   >  1:3467592/1‑36 (MQ=255)
                 aaaaTGGGTGCAGTAC.A.GTTGATAG..AAATATCggtc                 <  1:7726538/36‑4 (MQ=18)
                  aaaTGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCt                 >  1:4851825/1‑36 (MQ=255)
                  aaaTGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCt                 <  1:669316/36‑1 (MQ=255)
                   aaTGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCTg                <  1:123425/36‑1 (MQ=255)
                   aaTGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCTg                >  1:14174/1‑36 (MQ=255)
                    aTGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCTgg               <  1:6606959/36‑1 (MQ=255)
                      gggTGCAGTAC.AGGTTGATAG..AAATATCGCCTgggg             <  1:788543/36‑1 (MQ=255)
                      gggTGCAGTAC.AGGTTGATAG..AAATATCGCCTgggg             <  1:6852583/36‑1 (MQ=255)
                       ggTGCAGTAC.AGGTTGATAG..AAATATCGCCTGGGGc            >  1:995394/1‑36 (MQ=255)
                       ggTGCAGTAC.AGGTTGATAG..AAATATCGCCTGGGGc            >  1:3527456/1‑36 (MQ=255)
                        gTGCAGTAC.AGGTTTATAG..AAATATCGCCTGGGGCt           <  1:3760599/36‑1 (MQ=255)
                         tGCAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCtt          >  1:5110309/1‑36 (MQ=255)
                          gCAGTACAAGGTTGATAG..AAATATCGCCTGGGGCtt          <  1:2007250/36‑1 (MQ=25)
                          gCAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCttt         >  1:5103364/1‑36 (MQ=255)
                           cAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCtttt        >  1:5130693/1‑36 (MQ=255)
                           cAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCtttt        <  1:509843/36‑1 (MQ=255)
                           cAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCtttt        <  1:1115675/36‑1 (MQ=255)
                            aGTAC.AGGTTGATAGTAAAATATCACTCGGGGCttt         <  1:1193869/36‑1 (MQ=25)
                            aGTAC.AGGTTGATAG..AAATATCGCCTGGGGCTTTTg       <  1:4374519/36‑1 (MQ=255)
                             gTAC.AGGTTGATAG..AAATATCGCCTGGGGCTTTTGt      >  1:5578499/1‑36 (MQ=255)
                             gTAC.AGGTTGATAG..AAATATCGCCTGGGGCTTTTGt      >  1:16828/1‑36 (MQ=255)
                               aC.AGGTTGATAG..AAATATCGCCTGGGGCTTTTTTcc    >  1:3376607/1‑36 (MQ=255)
                                c.AGGTTGATAG..AAATATCGCCTGGGGCTTTTGTCCa   <  1:7378032/36‑1 (MQ=255)
                                c.AGGTTGATAG..AAATATCGCCTGGGGCTTTTGTCCa   >  1:7766192/1‑36 (MQ=255)
                                  aGGTTGATAG..AAATATCGCCTGGGGCTTTTGTCCAt  <  1:156915/36‑1 (MQ=255)
                                  aGGTTGATAG..AAATATCGCCTGGGGCTTTTGTCCAt  <  1:4652710/36‑1 (MQ=255)
                                   |                                    
CCATTTCATCGTCCAACAAAATGGGTGCAGTAC.AGGTTGATAG..AAATATCGCCTGGGGCTTTTGTCCAT  >  REL606/2775844‑2775912

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG