Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,732,014 C→A 34.6% R91L (CGC→CTC)  hypF ← carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,732,0140C→A34.6% 38.8 0.0 ‑0.353272hypFcarbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  new base (9/9):  ref base (16/18):  total (26/27)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.80e-01

AATTTGCGTATTCATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAGATA  >  REL606/2731980‑2732049
                                  |                                   
aaTTTGCGTATTCATGGTGCCGCCTGTGCTCTGGCg                                    >  1:3980559/1‑36 (MQ=255)
 aTTTGCGTATTCATGGTGCCGCCTGTGCTCTGGCgg                                   >  1:4237599/1‑36 (MQ=255)
 aTTTGCGTATTCATGGTGCCGCCTGTGCTCTGGAgg                                   <  1:2947962/36‑1 (MQ=255)
 aTTTGCGTATACCTGGTGCCGCCTGTGCTCTGGCgg                                   <  1:5885652/36‑1 (MQ=38)
  tttGCGTATTCATGGTGCCGCCTGTGCTCTGGAGGa                                  <  1:2104377/36‑1 (MQ=255)
    tGCGTATTCATGGTGCCGCCTGTGCTCTGGCGGGTa                                >  1:5894268/1‑36 (MQ=255)
    tGCGTATTCATGGTGCCGCCTGTGCTCTGGCGGATa                                <  1:2658030/36‑1 (MQ=255)
     gCGTATTCATGGTGCCGCCTGTGCTCTGGCGGATAg                               >  1:1866489/1‑36 (MQ=255)
       gTATTCATGGTGCCGCCTGTGCTCTGGAGGATAGTg                             >  1:3488759/1‑36 (MQ=255)
          ttCATGGTGCCGCCTGTGCTCTGGAGGATAGTGAAc                          <  1:5096116/36‑1 (MQ=255)
          ttCATGGTGCCGCCTGTGCTCTGGAGGATAGTGAAc                          <  1:3786518/36‑1 (MQ=255)
          ttCATGCTGCCGCCTGTGCTCTTGCGGATAGTGAAc                          <  1:3797970/36‑1 (MQ=37)
           tCATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACt                         >  1:533321/1‑36 (MQ=255)
            cATGGTGGCGCCTGTGCTCTGGCGGATAGTGAACTc                        <  1:1355717/36‑1 (MQ=255)
            cATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTc                        <  1:1131084/36‑1 (MQ=255)
            cATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTc                        <  1:1068600/36‑1 (MQ=255)
             attgTGCCGCCTGTGCTCTGGATGATAGTGAACTCa                       <  1:2333864/33‑1 (MQ=25)
             aTGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTCa                       >  1:840449/1‑36 (MQ=255)
             aTGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTCa                       >  1:331364/1‑36 (MQ=255)
              tGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTCAg                      >  1:1661207/1‑36 (MQ=255)
               ggcgccgccTGTGCTCTGGCGGATAGTGAACTCAGt                     >  1:4871029/4‑36 (MQ=255)
               ggTGCCGCCTGTGCTCTGGAGGATAGTGAACTCAGt                     <  1:1814769/36‑1 (MQ=255)
               ggTACCGCCTGTGCTCTGGAGGATAGTGAACTCAGt                     <  1:3448461/36‑1 (MQ=38)
                gTGTCGCTTGTGCTCTGGCGGATAGTGAACTCAGTg                    <  1:5531149/36‑1 (MQ=38)
                  gccgccTGTGCTCTGGCGGATAGTGAACTCAGTggg                  >  1:3863485/1‑36 (MQ=255)
                   ccgccTGTGCTCTGGAGGATAGTGAACTCAGTGGGc                 >  1:2227422/1‑36 (MQ=255)
                    cgccTGTGCTCTGGCGGATAGTGAACTCAGTGGGCa                <  1:3140052/36‑1 (MQ=255)
                    cgccTGTGCTCTGGCGGATAGTGAACTCAGTGGGCa                <  1:3773138/36‑1 (MQ=255)
                    cgccTGTGCTCTGGAGGATAGTGAACTCAGTGGGCa                >  1:3216363/1‑36 (MQ=255)
                      ccTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGt              <  1:609668/36‑1 (MQ=255)
                      ccTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGt              >  1:2421001/1‑36 (MQ=255)
                      ccTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGt              <  1:1453193/36‑1 (MQ=255)
                      ccTGTGCTCTGGAGGATAGTGAACTCAGTGGGCAGt              <  1:3120960/36‑1 (MQ=255)
                      ccTGTGCTCTGGAGGATAGTGAACTCAGTGGGCAGt              >  1:4862675/1‑36 (MQ=255)
                        tgtgCTCTGGTGGATGGTGAACTCAGTGTGCTGTtg            >  1:5717799/1‑36 (MQ=16)
                        tgtgCTCTGGCGGATAGTGAACTCAGTGGGCAGTtg            <  1:5886001/36‑1 (MQ=255)
                        tgtgCTCTGGCGGATAGTGAACTCAGTGGGCAGTtg            >  1:256906/1‑36 (MQ=255)
                         gtgCTCTGGCGGATAGTGAACTCAGTGGGCAGTtgt           <  1:184606/36‑1 (MQ=255)
                         gtgCTCTGGCGGATAGTGAACTCAGTGGGCAGTtgt           <  1:1891702/36‑1 (MQ=255)
                          tgCTCTGGCGGATAGTGAACTCAGTGGGCTGTtgtg          >  1:5201182/1‑36 (MQ=255)
                          tgCTCTGGCGGATAGTGAACTCAGTGGGCAGTtgtg          >  1:4028307/1‑36 (MQ=255)
                           gCTCTGGCGGATAGTGAACTCAGTGGGCAGTTGTGa         >  1:14604/1‑36 (MQ=255)
                           gCTCTGGCGGATAGTGAACTCAGTGGGCAGTTGTGa         <  1:3591689/36‑1 (MQ=255)
                           gCTCTGGAGGATAGTGAACTAAGTGGGCAGTTGTGa         >  1:834341/1‑36 (MQ=255)
                             tctGGAGGATAGTGAACTCAGTGGGCAGTTGTGAcc       >  1:4844731/1‑36 (MQ=255)
                             tatGGAGGATAGTGAACTCAGTGGGCAGTTGTGAcc       >  1:647392/3‑36 (MQ=38)
                              ctGGAGGATAGTGAACTCAGTGGGCAGTTGTGACCa      <  1:4151208/36‑1 (MQ=255)
                               tGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAg     <  1:707126/36‑1 (MQ=255)
                               tGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAg     >  1:1101046/1‑36 (MQ=255)
                               tGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAg     >  1:1850886/1‑36 (MQ=255)
                                ggCGGATAGTGAACTCAGTGGGCAGTTGTGACCAGa    <  1:3601633/36‑1 (MQ=255)
                                 gCGGATAGTGAACTCAGTGGGCAGTTGTGACCAGAt   <  1:571573/36‑1 (MQ=255)
                                 gCGGATAGTGAACTCAGTGGGCAGTTGTGACCAGAt   >  1:4104306/1‑36 (MQ=255)
                                  cGGATAGTGAACTCAGTGGGCAGTTGTGACCGGATa  <  1:4546855/36‑1 (MQ=255)
                                  |                                   
AATTTGCGTATTCATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAGATA  >  REL606/2731980‑2732049

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG