Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,902,276 (T)7→8 coding (7/342 nt) yebY ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,902,2691.T98.4% 212.4 / ‑4.0 63S5N (AGT→AAT) yebYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (23/39);  minor base C (0/1);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCG  >  REL606/1902222‑1902314
                                                |                                              
ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttt                                            <  1:5890439/51‑1 (MQ=255)
ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttt                                            >  1:2509183/1‑51 (MQ=255)
ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttt                                            <  1:6421618/51‑1 (MQ=255)
  gccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttt                                          <  1:3910452/51‑1 (MQ=255)
   ccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑tttttt                                         <  1:4801937/51‑1 (MQ=255)
   ccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑tttttt                                         >  1:1439095/1‑51 (MQ=255)
   ccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑tttttt                                         <  1:7225971/51‑1 (MQ=255)
    cAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttttt                                        <  1:1292584/51‑1 (MQ=255)
     aGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑tttttttt                                       >  1:1982241/1‑50 (MQ=255)
     aGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑tttttttt                                       >  1:1541665/1‑50 (MQ=255)
     aGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑tttttttt                                       <  1:2599079/51‑2 (MQ=255)
      gcgcCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttttttc                                      >  1:1955068/1‑49 (MQ=255)
       cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttttttca                                     <  1:5589972/51‑4 (MQ=255)
       cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttttttca                                     <  1:4841127/51‑4 (MQ=255)
        gccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttttttcat                                    <  1:3846242/51‑5 (MQ=255)
        gccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑ttttttttcat                                    <  1:5657918/51‑5 (MQ=255)
          cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTtcatca                                   <  1:2917629/51‑1 (MQ=255)
          cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTtcatca                                   <  1:5324073/51‑1 (MQ=255)
          cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTtcatca                                   >  1:1258195/1‑51 (MQ=255)
          cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTtcatca                                   >  1:339478/1‑51 (MQ=255)
           gcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTcatcat                                  >  1:6362357/1‑51 (MQ=255)
             tGCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCAttt                                >  1:5881591/1‑51 (MQ=255)
              gCAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCAtttt                               <  1:2915078/51‑1 (MQ=255)
               cAGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTc                              <  1:1443537/51‑1 (MQ=255)
                aGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCg                             <  1:1039758/51‑1 (MQ=255)
                aGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCg                             <  1:5814469/51‑1 (MQ=255)
                aGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCg                             <  1:6578346/51‑1 (MQ=255)
                aGAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCg                             <  1:6788039/51‑1 (MQ=255)
                 gAACTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGt                            <  1:6523394/51‑1 (MQ=255)
                  aaCTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTa                           >  1:1807835/1‑51 (MQ=255)
                  aaCTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTa                           <  1:2605834/51‑1 (MQ=255)
                  aaCTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTa                           <  1:6720982/51‑1 (MQ=255)
                    cTGGTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAAt                         >  1:561828/1‑51 (MQ=255)
                       gTGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCa                      <  1:5678793/51‑1 (MQ=255)
                        tGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCAt                     <  1:49642/51‑1 (MQ=255)
                        tGAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCAt                     >  1:1181193/1‑51 (MQ=255)
                         gAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATc                    >  1:2743537/1‑51 (MQ=255)
                         gAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATc                    <  1:5538854/51‑1 (MQ=255)
                         gAGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATc                    <  1:2707882/51‑1 (MQ=255)
                          aGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCa                   <  1:7243921/51‑1 (MQ=255)
                          aGTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCa                   <  1:3953788/51‑1 (MQ=255)
                           gTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaa                  >  1:3478643/1‑51 (MQ=255)
                           gTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaa                  <  1:3868198/51‑1 (MQ=255)
                           gTAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaa                  <  1:7084155/51‑1 (MQ=255)
                            tAACAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaa                 <  1:4491969/51‑1 (MQ=255)
                             aaCAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaaa                <  1:6863515/51‑1 (MQ=255)
                             aaCAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaaa                <  1:2475240/51‑1 (MQ=255)
                             aaCAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaaa                >  1:6139170/1‑51 (MQ=255)
                              aCAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaaaa               <  1:5589371/51‑1 (MQ=255)
                              aCAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaaaa               >  1:5331628/1‑51 (MQ=255)
                              aCAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCaaaaa               <  1:1931002/51‑1 (MQ=255)
                               cAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:4410063/51‑1 (MQ=255)
                               cAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:524764/51‑1 (MQ=255)
                               cAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAg              >  1:2603473/1‑51 (MQ=255)
                               cAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:4727577/51‑1 (MQ=255)
                               cAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAg              >  1:5859620/1‑51 (MQ=255)
                               cAGAAACGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:3264141/51‑1 (MQ=255)
                                  aaaCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGtt           >  1:4660513/1‑51 (MQ=255)
                                   aaCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          <  1:244600/51‑1 (MQ=255)
                                   aaCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          <  1:4315714/51‑1 (MQ=255)
                                   aaCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          >  1:2050513/1‑51 (MQ=255)
                                   aaCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          >  1:1851953/1‑51 (MQ=255)
                                    aCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg         >  1:5807370/1‑51 (MQ=255)
                                    aCGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg         <  1:7240776/51‑1 (MQ=255)
                                     cGCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCga        <  1:541528/51‑1 (MQ=255)
                                     cGCCAGAATACCTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg         <  1:52517/51‑1 (MQ=255)
                                      gCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:6264312/1‑51 (MQ=255)
                                      gCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:5562835/1‑51 (MQ=255)
                                      gCCAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:5752624/1‑51 (MQ=255)
                                       ccAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:5828896/51‑1 (MQ=255)
                                       ccAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:5415585/51‑1 (MQ=255)
                                       ccAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      >  1:3204119/1‑51 (MQ=255)
                                       ccAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:4690758/51‑1 (MQ=255)
                                       ccAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:3993156/51‑1 (MQ=255)
                                       ccAGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:801868/51‑1 (MQ=255)
                                         aGAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGt    >  1:1081148/1‑51 (MQ=255)
                                          gAATACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTc   <  1:5459400/51‑1 (MQ=255)
                                           aaTACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCg  <  1:7234563/51‑1 (MQ=255)
                                           aaTACT‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCg  >  1:5925815/1‑51 (MQ=255)
                                                |                                              
CCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCG  >  REL606/1902222‑1902314

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