Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,893,551 +G intergenic (+6/‑50) kup → / → insJ‑5 potassium transporter/IS150 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,893,5501.G95.3% 155.9 / 0.4 43intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (24/17);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACC  >  REL606/3893475‑3893619
                                                                             |                                                                     
tGCGTGCACCAGATCGATTTGAAATCCCGCCAATCAGGGTTATCGAACTGGGTACTCAGGTNGAAATCTAACGCCa                                                                        >  1:1153291/1‑76 (MQ=255)
tGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                                                        >  1:863095/1‑76 (MQ=255)
 gCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAg                                                                      <  1:857365/76‑1 (MQ=255)
   gTGCACCAGATCAATTTGACATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAgga                                                                    >  2:749953/1‑74 (MQ=255)
    tGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAATNCCAggac                                                                   <  2:1084454/76‑4 (MQ=255)
               aaTTTGAAATCCCGCCACACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAGCGCCAGGACTGTACTGCAcc                                                         >  1:1325205/1‑76 (MQ=255)
                     aaaTCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCAtttt                                                   >  2:93334/1‑76 (MQ=255)
                         cccGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTttgttg                                               <  1:150977/76‑1 (MQ=255)
                           cGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGa                                             <  1:882792/76‑1 (MQ=255)
                            gCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGAc                                            >  1:337498/1‑76 (MQ=255)
                            gCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGCACTGCACCCATTTTGTTGGAc                                            >  1:421594/1‑76 (MQ=39)
                             ccAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACg                                           <  1:693866/76‑1 (MQ=255)
                              cAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGa                                          >  1:671915/1‑76 (MQ=255)
                               aaaCAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGAt                                         >  2:631543/1‑76 (MQ=255)
                                   aGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAAGGCCCGGACTGTACTGCACTCATTTTGTTGGACGATGaaa                                     >  2:1059051/1‑76 (MQ=25)
                                   aGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGaaa                                     >  1:740134/1‑76 (MQ=255)
                                     ggTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                                   <  2:244661/76‑1 (MQ=255)
                                       ttATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGa                                 <  2:301496/76‑1 (MQ=255)
                                       ttATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGa                                 >  1:824522/1‑76 (MQ=255)
                                       ttATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa                                >  2:135209/1‑76 (MQ=255)
                                        tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa                                >  2:502182/1‑76 (MQ=255)
                                         aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                               <  1:1222774/76‑1 (MQ=255)
                                          tCGAATGGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGCTGAAATGGAATa                              <  2:56562/76‑1 (MQ=21)
                                          tCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                              <  1:1109948/76‑1 (MQ=255)
                                          tCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                              <  2:1220117/76‑1 (MQ=255)
                                          tCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                              <  2:1220123/76‑1 (MQ=255)
                                          tCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                              <  2:384021/76‑1 (MQ=255)
                                           cGAACTGGGTACTCATGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                             >  2:719792/1‑76 (MQ=21)
                                           cGAACTGGGTACTCACGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                             >  2:719776/1‑76 (MQ=21)
                                            gAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGc                            >  2:1044491/1‑76 (MQ=255)
                                              aCTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc                          >  2:749273/1‑76 (MQ=255)
                                               cTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc                         >  2:777881/1‑76 (MQ=255)
                                                tGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                        >  1:140537/1‑76 (MQ=255)
                                                tGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                        <  2:1174168/76‑1 (MQ=255)
                                                  ggTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa                      >  1:864774/1‑76 (MQ=255)
                                                    tACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                    <  2:741055/76‑1 (MQ=255)
                                                     aCTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                   >  2:245831/1‑76 (MQ=255)
                                                            tCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg            <  2:159366/76‑1 (MQ=255)
                                                            tCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg            >  2:841337/1‑76 (MQ=255)
                                                            tCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg            >  1:944495/1‑76 (MQ=255)
                                                             cGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGc           >  2:32832/1‑76 (MQ=255)
                                                               naaTCTAACGCCAGGACTGTACTGCACCCATTTTCTTNGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         >  1:1308885/2‑76 (MQ=14)
                                                               aaaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         <  2:415223/76‑1 (MQ=255)
                                                               aaaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         <  1:380245/76‑1 (MQ=255)
                                                               aaaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         >  2:1007810/1‑76 (MQ=255)
                                                               aaaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         >  2:781267/1‑76 (MQ=255)
                                                                  tCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaaa      >  1:439395/1‑76 (MQ=255)
                                                                      aCGCCAGGACTGTACTGTAGCCTTTTTGTTGGACGATGAAATGGAATAGCCGCTAGTATGTCAAAGCCAAAATAcc  <  1:64603/76‑1 (MQ=11)
                                                                             |                                                                     
TGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACC  >  REL606/3893475‑3893619

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