Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,920,708 G→A A151V (GCG→GTG)  znuA ← high‑affinity zinc transporter periplasmic component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,920,7080GA97.4% 139.9 / ‑6.5 38A151V (GCG→GTG) znuAhigh‑affinity zinc transporter periplasmic component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/18);  minor base C (1/0);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTT  >  REL606/1920635‑1920772
                                                                         |                                                                
ggCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACt                                                                <  2:1076492/76‑1 (MQ=255)
  cGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTAt                                                              <  1:1154716/76‑1 (MQ=255)
  cGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTAt                                                              <  2:822234/76‑1 (MQ=255)
     cAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTAtctc                                                           <  2:610241/76‑1 (MQ=255)
     cAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTAtctc                                                           <  2:610235/76‑1 (MQ=255)
     cAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTAtctc                                                           <  1:277622/76‑1 (MQ=255)
     cAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTAtctc                                                           >  1:687564/1‑76 (MQ=255)
      aaGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATctct                                                          <  2:782236/76‑1 (MQ=255)
        gtgtGGCTCGACTTTGTGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTgg                                                        <  1:201474/73‑1 (MQ=255)
        gTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTgg                                                        <  1:201490/76‑1 (MQ=255)
         tttGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTggg                                                       <  2:1065980/76‑1 (MQ=255)
            ggCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGaa                                                    <  2:898378/76‑1 (MQ=255)
            ggCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGaa                                                    <  2:1010468/76‑1 (MQ=255)
                cGACTTTGCGGCATAAGTTCCACTAATTTTCCATGAATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGcc                                                >  1:124725/1‑76 (MQ=255)
                cGACTTTGCGGCATAAGTTCCACTAATTTTCCATGAATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGcc                                                >  1:124717/1‑76 (MQ=255)
                    tttGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAg                                            <  1:647097/76‑1 (MQ=255)
                      tGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCCCTATCTCTGGGGAAAGCCAAAGAt                                          >  2:824033/1‑76 (MQ=255)
                      tGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGAt                                          >  2:824038/1‑76 (MQ=255)
                      tGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGAt                                          <  1:986332/76‑1 (MQ=255)
                        cGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGc                                        >  2:535237/1‑76 (MQ=255)
                             tAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGtt                                   <  1:630702/76‑1 (MQ=255)
                             tAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGtt                                   >  1:470763/1‑76 (MQ=255)
                              aaGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTg                                  >  1:281213/1‑76 (MQ=255)
                              aaGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTg                                  >  2:1086574/1‑76 (MQ=255)
                               aGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGa                                 >  1:1175480/1‑76 (MQ=255)
                                 ttCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGaaa                               >  1:280643/1‑76 (MQ=255)
                                    cACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCg                            <  1:155554/76‑1 (MQ=255)
                                       tAATTTTCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCg                         >  2:1058270/1‑76 (MQ=255)
                                          ttttCCATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGa                      >  1:995897/1‑76 (MQ=255)
                                              ccATGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGt                  >  1:221443/1‑76 (MQ=255)
                                                aTGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGa                <  1:276776/76‑1 (MQ=255)
                                                 tGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGAt               >  2:654600/1‑76 (MQ=255)
                                                 tGGATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGAt               <  1:1042495/76‑1 (MQ=255)
                                                   gATTGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCt             >  1:647747/1‑76 (MQ=255)
                                                      tGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCg          >  2:1216375/1‑76 (MQ=255)
                                                      tGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCg          >  1:340641/1‑76 (MQ=255)
                                                      tGCAACCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCg          >  1:340633/1‑76 (MQ=255)
                                                         aaCCGCTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCa       >  1:462778/1‑76 (MQ=255)
                                                              cTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACtttt  >  2:738636/1‑76 (MQ=255)
                                                              cTGTAGCCCGCACTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACtttt  >  1:577049/1‑76 (MQ=255)
                                                                         |                                                                
GGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTT  >  REL606/1920635‑1920772

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