Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,111,966 (G)7→8 intergenic (‑235/+167) insE‑1 ← / ← serX putative transposase‑related protein/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,111,9591.G97.1% 133.8 / ‑3.0 35intergenic (‑228/+174)insE‑1/serXputative transposase‑related protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (0/34);  total (0/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAA  >  REL606/1111885‑1112025
                                                                            |                                                                     
aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                      <  2:2651229/76‑1 (MQ=255)
aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                      <  2:4923373/76‑1 (MQ=255)
aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                      <  1:4275125/76‑1 (MQ=255)
aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                      <  2:3630741/76‑1 (MQ=255)
aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                                                      <  2:1026092/76‑1 (MQ=255)
 tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgg                                                                     <  1:4538726/76‑1 (MQ=255)
 tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgg                                                                     <  1:4078245/76‑1 (MQ=255)
 tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgg                                                                     <  2:3263075/76‑1 (MQ=255)
  ccccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAggg                                                                   >  1:2886423/1‑76 (MQ=255)
  ccccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggg                                                                   >  2:6050340/1‑76 (MQ=255)
   cccccGGGTTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGCGTTGAggggg                                                                <  2:5249587/76‑1 (MQ=14)
     cccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgggggg                                                               <  1:3004755/76‑1 (MQ=255)
      ccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAggggggg                                                             <  2:2920566/76‑1 (MQ=255)
       cGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGggggggggg                                                            >  2:1094050/1‑75 (MQ=255)
       cGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgggggggg                                                            <  2:989926/76‑2 (MQ=255)
        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAggggaggga                                                           <  1:3645484/76‑3 (MQ=255)
        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggggggggg                                                           >  1:2560648/1‑74 (MQ=255)
         gtgtgtGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGTggggggggcg                                                          >  1:2660840/1‑73 (MQ=14)
                cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcg                                                    <  2:3735561/76‑1 (MQ=255)
                  atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcgct                                                  <  1:5356601/76‑1 (MQ=255)
                   tataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcgctt                                                 <  1:1184446/76‑1 (MQ=255)
                     taCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcgcttcg                                               <  1:4150706/76‑1 (MQ=255)
                             aaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGccca                                       <  2:4763269/76‑2 (MQ=255)
                             aaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCt                                       <  1:4053595/76‑1 (MQ=255)
                              aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTg                                      <  2:4323798/76‑1 (MQ=255)
                              aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTg                                      <  2:3804040/76‑1 (MQ=255)
                              aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTg                                      <  2:1988062/76‑1 (MQ=255)
                               aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGGTGAATCGCTTCGCccccccctc                                     >  1:4505050/1‑73 (MQ=14)
                                aGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCg                                    <  2:3759515/76‑1 (MQ=255)
                                 gCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGGGTGACTCCCTTCCCTCCCCCCGccg                                   >  2:6258318/1‑74 (MQ=9)
                                  cccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCggg                                  <  1:3958091/76‑1 (MQ=255)
                                    cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                <  2:988304/76‑1 (MQ=255)
                                    cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                <  1:1905939/76‑1 (MQ=255)
                                    cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                                <  1:881981/76‑1 (MQ=255)
                                     gTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                               <  2:5004222/76‑1 (MQ=255)
                                     gTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                               <  2:4497730/76‑1 (MQ=255)
                                     gTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                               <  2:2883827/76‑1 (MQ=255)
                                     gTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                               <  2:3766453/76‑1 (MQ=255)
                                         tttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCg                           <  2:5392910/76‑1 (MQ=9)
                                           tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCt                         <  1:5986046/76‑1 (MQ=255)
                                           tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCt                         <  2:431894/76‑1 (MQ=255)
                                              tACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc                      <  1:3579330/76‑1 (MQ=255)
                                                cGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTg                    <  2:1483641/76‑1 (MQ=255)
                                                 gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGc                   <  2:5234303/76‑1 (MQ=255)
                                                    ctctTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGtt                <  2:58295/76‑1 (MQ=255)
                                                     tcttctTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc               <  1:824248/76‑1 (MQ=255)
                                                     tcttctTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc               <  1:2148659/76‑1 (MQ=255)
                                                     tcttctTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCa              <  1:3100383/76‑1 (MQ=35)
                                                      cttcttTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCa              <  1:4318196/76‑1 (MQ=255)
                                                        tcttTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACg            <  1:4152219/76‑1 (MQ=255)
                                                        tcttTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACg            <  1:4889492/76‑1 (MQ=255)
                                                        tcttTAAACACGGCGGCGGGGGGGGAGCTTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACg            <  2:1964623/76‑1 (MQ=14)
                                                            tAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCt        <  2:5328111/76‑1 (MQ=255)
                                                            tAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCt        <  2:584086/76‑1 (MQ=255)
                                                             aaaTATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTg       <  2:4391840/76‑1 (MQ=255)
                                                             aaaTATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTg       <  1:6258318/76‑1 (MQ=255)
                                                              actatGGCGGTGAGGTAGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGt      <  2:3024169/74‑1 (MQ=14)
                                                                catGGCGGTGAGGGGGGAGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTcc    <  2:269816/75‑1 (MQ=255)
                                                                 atGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCa   <  1:5487089/76‑1 (MQ=255)
                                                                 atGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCa   <  2:2277229/76‑1 (MQ=255)
                                                                 atGGCGGCTAGGGAGGGTATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCa   <  2:164344/76‑1 (MQ=11)
                                                                  tGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCaa  <  1:3790309/76‑1 (MQ=255)
                                                                            |                                                                     
ATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAA  >  REL606/1111885‑1112025

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