Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,902,269:1 +T 100% coding (14/342 nt) yebY ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,902,2691.T99.0% 326.1 / ‑4.2 95S5N (AGT→AAT) yebYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (35/59);  minor base C (0/1);  total (35/60)
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.17e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 123 total aligned reads displayed.

GTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCG  >  REL606/1902219‑1902314
                                                   |                                             
  gccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtt                                             <  1:10504730/51‑1 (MQ=255)
   ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttt                                            <  1:7658522/51‑1 (MQ=255)
   ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttt                                            <  1:6421618/51‑1 (MQ=255)
   ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttt                                            <  1:5890439/51‑1 (MQ=255)
     gccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttt                                          <  1:9900017/51‑1 (MQ=255)
     gccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttt                                          <  1:3910452/51‑1 (MQ=255)
      ccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtttttt                                         <  1:7225971/51‑1 (MQ=255)
      ccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtttttt                                         <  1:4801937/51‑1 (MQ=255)
      ccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtttttt                                         >  1:1439095/1‑51 (MQ=255)
       cAGCGCCGCTGCAGAACTGGTGAGTAATAGAAACGCCAGAATACttttttt                                        <  1:9096136/51‑1 (MQ=255)
       cAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttt                                        <  1:1292584/51‑1 (MQ=255)
        aGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtttttttt                                       <  1:2599079/51‑2 (MQ=255)
        aGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtttttttt                                       >  1:1541665/1‑50 (MQ=255)
        aGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACtttttttt                                       >  1:1982241/1‑50 (MQ=255)
          cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttca                                     <  1:10023668/51‑4 (MQ=255)
          cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttca                                     >  1:9022863/1‑48 (MQ=255)
          cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttca                                     <  1:8137111/51‑4 (MQ=255)
          cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttca                                     <  1:5589972/51‑4 (MQ=255)
          cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttca                                     <  1:4841127/51‑4 (MQ=255)
           gccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttcat                                    <  1:7371765/51‑5 (MQ=255)
           gccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttcat                                    <  1:3846242/51‑5 (MQ=255)
            ccgcTGCTGAACTGGTGAGTAACAGAAACGCCAGAATACttttttttcatc                                   <  1:9887133/51‑6 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   <  1:10147286/51‑1 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   <  1:10653807/51‑1 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   <  1:5324073/51‑1 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   >  1:1258195/1‑51 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   >  1:339478/1‑51 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   <  1:2917629/51‑1 (MQ=255)
             cgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                   <  1:8464140/51‑1 (MQ=255)
              gcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTcatcat                                  >  1:6362357/1‑51 (MQ=255)
              gcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTcatcat                                  <  1:7765391/51‑1 (MQ=255)
                tGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCAttt                                >  1:7658007/1‑51 (MQ=255)
                tGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCAttt                                <  1:8661122/51‑1 (MQ=255)
                tGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCAttt                                >  1:5881591/1‑51 (MQ=255)
                   aGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCg                             <  1:6788039/51‑1 (MQ=255)
                   aGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCg                             <  1:6578346/51‑1 (MQ=255)
                   aGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCg                             <  1:5814469/51‑1 (MQ=255)
                   aGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCg                             <  1:10325134/51‑1 (MQ=255)
                    gAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGt                            <  1:6523394/51‑1 (MQ=255)
                     aaCTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTa                           >  1:7990545/1‑51 (MQ=255)
                     aaCTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTa                           <  1:2605834/51‑1 (MQ=255)
                     aaCTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTa                           <  1:9653461/51‑1 (MQ=255)
                     aaCTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTa                           >  1:1807835/1‑51 (MQ=255)
                       cTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAAt                         >  1:561828/1‑51 (MQ=255)
                         ggTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATcc                       <  1:9552691/51‑1 (MQ=255)
                          gTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCa                      <  1:5678793/51‑1 (MQ=255)
                           tGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCAt                     >  1:1181193/1‑51 (MQ=255)
                           tGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCAt                     <  1:9337052/51‑1 (MQ=255)
                           tGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCAt                     <  1:49642/51‑1 (MQ=255)
                           tGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCAt                     >  1:10245804/1‑51 (MQ=255)
                            gAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATc                    >  1:2743537/1‑51 (MQ=255)
                            gAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATc                    <  1:5538854/51‑1 (MQ=255)
                            gAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATc                    <  1:2707882/51‑1 (MQ=255)
                             aGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCa                   <  1:7243921/51‑1 (MQ=255)
                             aGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCa                   <  1:3953788/51‑1 (MQ=255)
                              gTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaa                  <  1:7381486/51‑1 (MQ=255)
                              gTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaa                  <  1:7084155/51‑1 (MQ=255)
                              gTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaa                  <  1:3868198/51‑1 (MQ=255)
                               tAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaa                 <  1:4491969/51‑1 (MQ=255)
                                aaCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaa                <  1:8392796/51‑1 (MQ=255)
                                aaCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaa                <  1:7712898/51‑1 (MQ=255)
                                aaCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaa                <  1:2475240/51‑1 (MQ=255)
                                aaCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaa                <  1:6863515/51‑1 (MQ=255)
                                 aCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaaa               <  1:5589371/51‑1 (MQ=255)
                                 aCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaaa               <  1:1931002/51‑1 (MQ=255)
                                 aCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaaa               >  1:5331628/1‑51 (MQ=255)
                                 aCAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaaa               <  1:7968807/51‑1 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:3264141/51‑1 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              >  1:10779274/1‑51 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:524764/51‑1 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              >  1:2603473/1‑51 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              >  1:9831187/1‑51 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              <  1:4410063/51‑1 (MQ=255)
                                  cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg              >  1:5859620/1‑51 (MQ=255)
                                   aGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAgg             >  1:8741164/1‑51 (MQ=255)
                                   aGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAgg             >  1:7637228/1‑51 (MQ=255)
                                     aaaCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGtt           <  1:9528222/51‑1 (MQ=255)
                                     aaaCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGtt           >  1:4660513/1‑51 (MQ=255)
                                      aaCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          >  1:2050513/1‑51 (MQ=255)
                                      aaCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          <  1:244600/51‑1 (MQ=255)
                                      aaCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          <  1:4315714/51‑1 (MQ=255)
                                      aaCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTc          >  1:1851953/1‑51 (MQ=255)
                                       aCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg         <  1:7240776/51‑1 (MQ=255)
                                       aCGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg         >  1:5807370/1‑51 (MQ=255)
                                        cGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCga        <  1:541528/51‑1 (MQ=255)
                                        cGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCga        <  1:10127143/51‑1 (MQ=255)
                                         gCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:5752624/1‑51 (MQ=255)
                                         gCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:6264312/1‑51 (MQ=255)
                                         gCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:5562835/1‑51 (MQ=255)
                                         gCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:8935918/1‑51 (MQ=255)
                                         gCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgag       >  1:9065371/1‑51 (MQ=255)
                                          ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:801868/51‑1 (MQ=255)
                                          ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      >  1:3204119/1‑51 (MQ=255)
                                          ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:10488134/51‑1 (MQ=255)
                                          ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:4690758/51‑1 (MQ=255)
                                          ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCgaga      <  1:5415585/51‑1 (MQ=255)
                                           cAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGagag     >  1:9311387/1‑51 (MQ=255)
                                            aGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGt    <  1:7565328/51‑1 (MQ=255)
                                            aGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGt    >  1:1081148/1‑51 (MQ=255)
                                              aaTACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCg  >  1:5925815/1‑51 (MQ=255)
                                                   |                                             
GTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCG  >  REL606/1902219‑1902314

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