Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,067,926 G→A 100% A40V (GCA→GTA)  yghS ← predicted protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,067,9260GA99.2% 376.8 / ‑5.3 110A40V (GCA→GTA) yghSpredicted protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (69/40);  total (69/41)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 116 total aligned reads displayed.

ATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAA  >  REL606/3067878‑3067973
                                                |                                               
aTTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtact                                               <  1:6112183/51‑1 (MQ=255)
 ttCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactg                                              <  1:6090288/51‑1 (MQ=255)
 ttCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactg                                              <  1:3409742/51‑1 (MQ=255)
 ttCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactg                                              >  1:4097589/1‑51 (MQ=255)
  tCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactgc                                             >  1:9303749/1‑51 (MQ=255)
    gCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCa                                           >  1:1099011/1‑51 (MQ=255)
    gCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCa                                           >  1:3181213/1‑51 (MQ=255)
     cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          <  1:2713634/51‑1 (MQ=255)
     cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          >  1:5331681/1‑51 (MQ=255)
     cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          >  1:1256494/1‑51 (MQ=255)
     cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          <  1:7747212/51‑1 (MQ=255)
      ccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGt                                         >  1:1260505/1‑51 (MQ=255)
       cNGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGtt                                        >  1:5197768/1‑51 (MQ=255)
       cGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGtt                                        >  1:5170635/1‑51 (MQ=255)
        ggACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTc                                       <  1:6509822/51‑1 (MQ=255)
         gACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCa                                      <  1:7291065/51‑1 (MQ=255)
          aCGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAt                                     >  1:10760738/1‑51 (MQ=255)
           cGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAtt                                    <  1:10522513/51‑1 (MQ=255)
           cGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAtt                                    <  1:6112565/51‑1 (MQ=255)
           cGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAtt                                    >  1:9687036/1‑51 (MQ=255)
            gATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAttt                                   <  1:8438637/51‑1 (MQ=255)
            gATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAttt                                   <  1:7783435/51‑1 (MQ=255)
            gATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAttt                                   >  1:10879107/1‑51 (MQ=255)
            gATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAttt                                   <  1:8096030/51‑1 (MQ=255)
             aTTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTa                                  >  1:8486787/1‑51 (MQ=255)
             aTTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTa                                  >  1:6480181/1‑51 (MQ=255)
              ttGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAc                                 >  1:2299604/1‑51 (MQ=255)
              ttGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAc                                 >  1:8130621/1‑51 (MQ=255)
               tGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAcc                                >  1:9412936/1‑51 (MQ=255)
               tGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAcc                                >  1:32116/1‑51 (MQ=255)
               tGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAcc                                <  1:10904287/51‑1 (MQ=255)
                gCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCa                               <  1:5144940/51‑1 (MQ=255)
                gCCCGAGATAAATGTGTTCTGTTGGCATTCTTACCGCCAGTTCATTTACCa                               >  1:6536822/1‑51 (MQ=255)
                 cccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAg                              <  1:10855558/51‑1 (MQ=255)
                  ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             >  1:503045/1‑51 (MQ=255)
                  ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             <  1:3779506/51‑1 (MQ=255)
                  ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             >  1:3076228/1‑51 (MQ=255)
                  ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             <  1:7007762/51‑1 (MQ=255)
                  ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             <  1:3183008/51‑1 (MQ=255)
                   cGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCa                               >  1:685817/1‑48 (MQ=255)
                   cGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGc                            >  1:10422702/1‑51 (MQ=255)
                   cGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGc                            >  1:1581594/1‑51 (MQ=255)
                    gagagaAATGTGTTCTGTTGGCATTCTGGCTGCCAGTTCa                                      <  1:9893996/40‑1 (MQ=37)
                    gagaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCt                           >  1:2068505/1‑51 (MQ=255)
                     agaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCtt                          >  1:2582798/1‑51 (MQ=255)
                     agaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCtt                          <  1:5465507/51‑1 (MQ=255)
                     agaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCtt                          <  1:7603619/51‑1 (MQ=255)
                      gaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTg                         >  1:9689779/1‑51 (MQ=255)
                      gaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTg                         <  1:7594841/51‑1 (MQ=255)
                      gaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTg                         <  1:8638556/51‑1 (MQ=255)
                       aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:1595859/1‑51 (MQ=255)
                       aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:4931556/1‑51 (MQ=255)
                       aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:3221740/1‑51 (MQ=255)
                       aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:3887086/1‑51 (MQ=255)
                        tAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGcc                       >  1:7518813/1‑51 (MQ=255)
                        tAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGcc                       >  1:4945320/1‑51 (MQ=255)
                        tAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGcc                       >  1:10122305/1‑51 (MQ=255)
                        tAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGcc                       >  1:8104537/1‑51 (MQ=255)
                         aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGc                            <  1:4740010/45‑1 (MQ=255)
                         aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCg                      >  1:3973470/1‑51 (MQ=255)
                         aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCAATTACCAGGCTTGCCg                      >  1:2366290/1‑51 (MQ=255)
                          aaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGt                     <  1:7911683/51‑1 (MQ=255)
                          aaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGt                     >  1:2053874/1‑51 (MQ=255)
                          aaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGt                     <  1:9827743/51‑1 (MQ=255)
                           aTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTg                    >  1:5213986/1‑51 (MQ=255)
                            tgtgTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGa                   >  1:5172924/1‑51 (MQ=255)
                            tgtgTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGa                   <  1:8248674/51‑1 (MQ=255)
                             gtgtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAg                  >  1:6137849/1‑51 (MQ=255)
                             gtgtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAg                  >  1:2437602/1‑51 (MQ=255)
                               gtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAggg                >  1:1112455/1‑51 (MQ=255)
                                ttCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGt               <  1:7928887/51‑1 (MQ=255)
                                 tCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTc              <  1:9238176/51‑1 (MQ=255)
                                 tCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTc              <  1:4784104/51‑1 (MQ=255)
                                  cTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCg             >  1:4503790/1‑51 (MQ=255)
                                  cTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCg             >  1:5392741/1‑51 (MQ=255)
                                   tGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGa            >  1:6962258/1‑51 (MQ=255)
                                    gTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAt           >  1:6046415/1‑51 (MQ=255)
                                    gTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAt           >  1:4266834/1‑51 (MQ=255)
                                    gTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAt           <  1:3829046/51‑1 (MQ=255)
                                     ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:9573124/1‑51 (MQ=255)
                                     ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:3859866/1‑51 (MQ=255)
                                     ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:906839/1‑51 (MQ=255)
                                     ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:1819669/1‑51 (MQ=255)
                                     ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:9327317/1‑51 (MQ=255)
                                     ttGGCATTCTTACTGCCAATTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:3284694/1‑51 (MQ=255)
                                      tGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt         >  1:5585135/1‑51 (MQ=255)
                                      tGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt         <  1:10419536/51‑1 (MQ=255)
                                      tGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt         >  1:489772/1‑51 (MQ=255)
                                      tGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt         >  1:8006266/1‑51 (MQ=255)
                                        gCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTAc       >  1:7627935/1‑51 (MQ=255)
                                         cATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTAcc      <  1:5770520/51‑1 (MQ=255)
                                          aTTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCt     >  1:10013092/1‑51 (MQ=255)
                                          aTTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCt     <  1:10249110/51‑1 (MQ=255)
                                           ttCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTg    <  1:9613505/51‑1 (MQ=255)
                                           ttCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTg    >  1:6439857/1‑51 (MQ=255)
                                            tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTg    <  1:7647688/50‑1 (MQ=255)
                                            tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGa   >  1:4084602/1‑51 (MQ=255)
                                            tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGa   <  1:1400330/51‑1 (MQ=255)
                                            tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGa   >  1:10566843/1‑51 (MQ=255)
                                             cTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGaa  <  1:4054430/51‑1 (MQ=255)
                                                |                                               
ATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAA  >  REL606/3067878‑3067973

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: