Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,775,879 Δ1 bp intergenic (‑49/+200) insJ‑3 ← / ← cysH IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,775,8780G.98.3% 277.2 / ‑1.5 59intergenic (‑48/+201)insJ‑3/cysHIS150 hypothetical protein/phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (30/28);  minor base A (0/1);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGC  >  REL606/2775836‑2775924
                                          |                                              
gggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGa                                       <  1:1760044/51‑1 (MQ=255)
 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaa                                      <  1:6226600/51‑1 (MQ=255)
 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaa                                      <  1:2934350/51‑1 (MQ=255)
  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaaa                                     >  1:574784/1‑51 (MQ=255)
  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaaa                                     >  1:2693330/1‑51 (MQ=255)
  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGaaa                                     <  1:250300/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:2342345/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:3826564/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:838421/51‑1 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    >  1:2393717/1‑51 (MQ=255)
   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAat                                    <  1:1232261/51‑1 (MQ=255)
    tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAata                                   >  1:2700231/1‑51 (MQ=255)
     aTTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAatat                                  <  1:5198260/51‑1 (MQ=255)
      ttCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATc                                 >  1:3428358/1‑51 (MQ=255)
       tCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCg                                >  1:5669934/1‑51 (MQ=255)
       tCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCg                                <  1:4638243/51‑1 (MQ=255)
        ccATTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGc                               >  1:6568564/1‑51 (MQ=255)
          aTTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCt                             >  1:586746/1‑51 (MQ=255)
          aTTTCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCt                             <  1:3294411/51‑1 (MQ=255)
           tttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTg                            <  1:1358552/51‑1 (MQ=255)
            ttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgg                           >  1:2790543/1‑51 (MQ=255)
            ttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgg                           >  1:1171428/1‑51 (MQ=255)
            ttCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgg                           <  1:2481093/51‑1 (MQ=255)
             tCATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTggg                          >  1:183625/1‑51 (MQ=255)
              cATCGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTgggg                         >  1:3762330/1‑51 (MQ=255)
                 cGTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtt                      >  1:1920997/1‑51 (MQ=255)
                 cGTCCAACAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtt                      <  1:3124501/50‑1 (MQ=25)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     <  1:1974002/51‑1 (MQ=255)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     >  1:1586727/1‑51 (MQ=255)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     >  1:6204760/1‑51 (MQ=255)
                  gTCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCttt                     <  1:6017716/51‑1 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    >  1:6553812/1‑51 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    <  1:7280371/51‑1 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    <  1:2852418/51‑1 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    >  1:6002299/1‑51 (MQ=255)
                   tCCAACAAAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCtttt                    >  1:6384997/1‑51 (MQ=255)
                        caaAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa               <  1:5259812/51‑1 (MQ=255)
                        caaAATGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa               <  1:4624461/51‑1 (MQ=255)
                         aaaaTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt              >  1:5655888/1‑51 (MQ=255)
                          aaaTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc             <  1:2934459/51‑1 (MQ=255)
                           aaTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGt                  >  1:3755827/1‑45 (MQ=38)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           >  1:1008149/1‑51 (MQ=255)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           <  1:2382496/51‑1 (MQ=255)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           >  1:2196692/1‑51 (MQ=255)
                            aTGGGTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg           <  1:813699/51‑1 (MQ=255)
                               ggTGCAGTAAAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa        <  1:1396976/51‑1 (MQ=38)
                                gTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc       >  1:890665/1‑51 (MQ=255)
                                gTGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc       >  1:5432522/1‑51 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      >  1:4539840/1‑51 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      <  1:5757739/51‑1 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      >  1:1171085/1‑51 (MQ=255)
                                 tGCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc      <  1:6535078/51‑1 (MQ=255)
                                  gCAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt     >  1:2878867/1‑51 (MQ=255)
                                   cAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc    >  1:324949/1‑51 (MQ=255)
                                   cAGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc    <  1:806368/51‑1 (MQ=255)
                                    aGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg   >  1:4874273/1‑51 (MQ=255)
                                    aGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg   <  1:5839725/51‑1 (MQ=255)
                                    aGTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg   <  1:5787157/51‑1 (MQ=255)
                                     gTACAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc  >  1:6077982/1‑51 (MQ=255)
                                          |                                              
GGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGC  >  REL606/2775836‑2775924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: