Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,749,029 G→A V292M (GTG→ATG)  hypE → carbamoyl phosphate phosphatase, hydrogenase 3 maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,749,0290GA97.9% 173.9 / ‑3.1 48V292M (GTG→ATG) hypEcarbamoyl phosphate phosphatase, hydrogenase 3 maturation protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (26/21);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

CAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCC  >  REL606/2748956‑2749100
                                                                         |                                                                       
caacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatg                                                                       >  1:1491915/1‑76 (MQ=255)
 aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatgg                                                                      >  1:1278274/1‑76 (MQ=255)
 aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatgg                                                                      >  2:224956/1‑76 (MQ=255)
 aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatgg                                                                      >  2:1424833/1‑76 (MQ=255)
    gcgGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAAtggtgg                                                                   >  2:232556/1‑76 (MQ=255)
     cgGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGa                                                                  >  1:1679015/1‑76 (MQ=255)
      ggCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGaa                                                                 <  2:988730/76‑1 (MQ=255)
      ggCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGaa                                                                 >  2:1713352/1‑76 (MQ=255)
         agagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGt                                                              >  1:932958/1‑76 (MQ=255)
          gagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTa                                                             >  1:116170/1‑76 (MQ=255)
          gagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTa                                                             >  2:1136209/1‑76 (MQ=255)
              caGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAgg                                                         <  1:210592/75‑1 (MQ=255)
              aaGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAgg                                                         >  1:1701017/1‑76 (MQ=255)
              aaGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAgg                                                         >  2:1345377/1‑76 (MQ=255)
               aGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGgt                                                        <  2:176261/76‑1 (MQ=255)
                 tGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGgtgt                                                      <  1:1401654/76‑1 (MQ=255)
                     ggCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTcac                                                  <  2:1396502/76‑3 (MQ=255)
                       cAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCt                                                <  2:470005/76‑1 (MQ=255)
                       cAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCt                                                >  1:1367315/1‑76 (MQ=255)
                        ggCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCtt                                               <  2:856945/75‑1 (MQ=255)
                          cGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGc                                             <  2:1545913/76‑1 (MQ=255)
                          cGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGc                                             >  1:545120/1‑76 (MQ=255)
                          cGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGc                                             <  2:1179158/76‑1 (MQ=255)
                           gTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGcc                                            >  1:1346567/1‑76 (MQ=255)
                           gTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGcc                                            >  1:1451582/1‑76 (MQ=255)
                            ttACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCg                                           <  1:1713352/76‑1 (MQ=255)
                            ttACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCg                                           >  2:136215/1‑76 (MQ=255)
                              aCATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGt                                         >  2:1349832/1‑76 (MQ=255)
                                aTTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCt                                       >  1:1826425/1‑76 (MQ=255)
                                aTTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCt                                       >  1:312612/1‑76 (MQ=255)
                                  tCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGt                                     >  1:1021840/1‑76 (MQ=255)
                                     cATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATg                                  <  1:1371756/76‑1 (MQ=255)
                                     cATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATg                                  <  2:1056926/76‑1 (MQ=255)
                                      gtCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTAGTCTTGCCGGCCGGTATgg                                 >  1:50311/2‑76 (MQ=255)
                                      aTCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATgg                                 >  1:764827/1‑76 (MQ=255)
                                      aTCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATgg                                 <  2:1712366/76‑1 (MQ=255)
                                        ccACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCg                               <  2:951470/76‑1 (MQ=255)
                                           cTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTgg                            >  1:973321/1‑75 (MQ=255)
                                           cTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGa                            >  2:1452155/1‑76 (MQ=255)
                                            tGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGaa                           >  1:1732215/1‑76 (MQ=255)
                                                 aaaGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGaa                      >  2:257765/1‑76 (MQ=255)
                                                  aaGACGCGGCGCTGATTGGGGAAAGGGTGGAACGTAAAGGTGTTCGTCTTGCTGGTCTGTATGGCGTGAAACGAAc                     >  2:918551/1‑76 (MQ=255)
                                                      cgcgGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTc                 >  2:603302/1‑76 (MQ=255)
                                                         ggcgCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGAt              <  1:379524/76‑1 (MQ=255)
                                                            gcTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTa           <  2:87239/76‑1 (MQ=255)
                                                            gcTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTa           <  2:1293057/76‑1 (MQ=255)
                                                              tGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTAcc         <  1:317984/76‑1 (MQ=255)
                                                                aTTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac       >  1:196307/1‑76 (MQ=255)
                                                                aTTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac       <  1:1476606/76‑1 (MQ=255)
                                                                aTTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac       <  1:1575413/76‑1 (MQ=255)
                                                                   ggTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACg    <  2:495705/76‑1 (MQ=255)
                                                                     tGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGcc  >  1:55441/1‑76 (MQ=255)
                                                                         |                                                                       
CAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCC  >  REL606/2748956‑2749100

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