Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,769,363 G→T 7.9% intergenic (‑246/+171) ydiV ← / ← nlpC hypothetical protein/predicted lipoprotein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,769,3630G→T 7.9% 2.3 ‑2.7 ‑0.038intergenic (‑246/+171)ydiV/nlpChypothetical protein/predicted lipoprotein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (13/22):  total (16/22)
Fisher's exact test strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.51e-01

TAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACT  >  REL606/1769328‑1769398
                                   |                                   
tAAAAGGGATGAGGCCATAAGACATAACAATTACAg                                     >  1:5045110/1‑36 (MQ=255)
tAAAAGGGATGAGGCCATAAGACATAACAATTACAg                                     >  1:3421144/1‑36 (MQ=255)
  aaaGGGATGAGGCAATAAGACATAACAATTACAGaa                                   <  1:6961199/36‑1 (MQ=255)
    aGGGATGAGGCCATAAGACATAACAATTACAGAAgg                                 >  1:5948961/1‑36 (MQ=255)
       gATGAGGCCATAAGACATAACAATTACAGAAGGAGt                              >  1:5624273/1‑36 (MQ=255)
        aTGAGGCCATAAGACATAACAATTACAGAAGGAGTa                             <  1:6686597/36‑1 (MQ=255)
         tGAGGACATAAGACATAACAATTACAGAAGGAGTaa                            <  1:207538/36‑1 (MQ=255)
          gAGGCCATAAGACATAACAATTACAGAAGGAGTAAc                           <  1:1338040/36‑1 (MQ=255)
           aGGCCATAAGACATAACAATTACAGAAGGAGTAACt                          >  1:4007867/1‑36 (MQ=255)
           aGGCCACAAGACATAACAATTACAGAAGGAGTAACt                          <  1:1363309/36‑1 (MQ=255)
            ggCCATAAGACATAACAATTACAGAAGGAGTAACtt                         <  1:6735618/36‑1 (MQ=255)
            ggCCATAAGACATAACAATTACAGAAGGAGTAACtt                         <  1:5609309/36‑1 (MQ=255)
             gCCATAAGACATAACAATTACAGAAGGAGTAACttt                        <  1:1408070/36‑1 (MQ=255)
             gCCATAAGACATAACAATTACAGAAGGAGTAACttt                        <  1:1502607/36‑1 (MQ=255)
                aTAAGACATAACAATTACAGAAGGAGTAACTTTCAt                     >  1:594939/1‑36 (MQ=255)
                  aaGACATAACAATTACAGAAGGAGTAACTTTCAttt                   >  1:7235124/1‑36 (MQ=255)
                     aCATAACAATTACAGAAGGAGTAACTTTCATTTGtt                <  1:8158934/36‑1 (MQ=255)
                     aCATAACAATTACAGAAGGAGTAACTTTCATTTGtt                >  1:5951648/1‑36 (MQ=255)
                       aTAACAATTACAGAAGGAGTAACTTTCATTTGTTcc              <  1:678272/36‑1 (MQ=255)
                       aTAACAATTACAGAAGGAGTAACTTTCATTTGTTcc              <  1:34354/36‑1 (MQ=255)
                        tAACAATTACAGAAGGAGTAACTTTCATTTGTTCCa             <  1:3387709/36‑1 (MQ=255)
                        tAACAATTACAGAAGGAGTAACTTTCATTTGTTCCa             >  1:6604743/1‑36 (MQ=255)
                         aaCAATTACAGAAGGAGTAACTTTCATTTGTTCCAt            <  1:5263296/36‑1 (MQ=255)
                         aaCAATTACAGAAGGAGTAACTTTCATTTGTTCCAt            <  1:842658/36‑1 (MQ=255)
                         aaCAATTACAGAAGGAGTAACTTTAAATTCTTCCAt            >  1:6122889/1‑36 (MQ=37)
                          aCAATTACAGAAGGAGTAACTTTCATTTGTTCCATg           >  1:3802091/1‑36 (MQ=255)
                           cAATTCCAGAAGGAGTAACTTTCATTTGTTCCATGt          <  1:8064564/36‑1 (MQ=255)
                           cAATTACAGAAGGAGTAACTTTCATTTGTTCCATGt          <  1:5173619/36‑1 (MQ=255)
                           cAATTACAGAAGGAGTAACTTTCATTTGTTCCATGt          <  1:4488040/36‑1 (MQ=255)
                            aaTTACAGAAGGAGT.ACTTTCATTTGTTCCATGTTa        >  1:668741/1‑36 (MQ=38)
                             aTTACAGAAGGAGTAACTTTCATTTGTTCCATGTTa        <  1:873868/36‑1 (MQ=255)
                               tACATAAGGAGTAACTTTCATTTGTTCCATGTTAAc      >  1:5637944/1‑36 (MQ=255)
                               tACATAAGGAGTAACTTTCATTTGTTCCATGTTAAc      >  1:8395178/1‑36 (MQ=255)
                               tACATAAGGAGTAACTTTCATTTGTTCCATGTTAAc      >  1:4595815/1‑36 (MQ=255)
                               tACAGAAGGAGTAACTTTCATTTGTTCCATGTTAAc      <  1:8366575/36‑1 (MQ=255)
                               tACAGAAGGAGTAACTTTCATTTGTTCCATGTTAAc      >  1:5487900/1‑36 (MQ=255)
                                aCAGAAGGAGTAACTTTCATTTGTTCCATGTTAAcc     <  1:352179/36‑1 (MQ=255)
                                aCAGAAGGAGTAACTTGCATTTGTTCCATGTTAAcc     >  1:88127/1‑36 (MQ=255)
                                 cAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCa    <  1:6426486/36‑1 (MQ=255)
                                  aGAAGGAGTAACTTTCATTTGTTCCATGTTAACCAc   >  1:419523/1‑36 (MQ=255)
                                   gAAGGAGTAACTTTCATTTGTTCCATGTTAACCACt  <  1:507496/36‑1 (MQ=255)
                                   |                                   
TAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACT  >  REL606/1769328‑1769398

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG