Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,396,897 A→T I295N (ATC→AAC)  yjfR ← predicted L‑ascorbate 6‑phosphate lactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,396,8970AT98.7% 586.6 / 0.1 149I295N (ATC→AAC) yjfRpredicted L‑ascorbate 6‑phosphate lactonase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (54/93);  total (54/95)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 152 total aligned reads displayed.

CCGCCAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCATTCAGG  >  REL606/4396825‑4396968
                                                                        |                                                                       
cgcccAACCTGCCAGATAAACGGCTCGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTtct                                                                      <  2:2044977/73‑1 (MQ=255)
ccgccACCCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTtct                                                                      <  1:1596178/76‑1 (MQ=255)
ccgccAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTtct                                                                      >  2:3053043/1‑76 (MQ=255)
cccccAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTtct                                                                      <  2:51700/73‑1 (MQ=255)
 cgccAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTtctc                                                                     >  1:4287505/1‑76 (MQ=255)
  gccAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTctct                                                                    >  1:3760728/1‑76 (MQ=255)
   ccAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCtt                                                                   >  1:4948424/1‑76 (MQ=255)
    cAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTg                                                                  >  2:408509/1‑76 (MQ=255)
    cAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTg                                                                  <  2:2631408/76‑1 (MQ=255)
      aCCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCg                                                                <  2:1922881/76‑1 (MQ=255)
      aCCTGCCAGAAAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCg                                                                <  1:860816/76‑1 (MQ=255)
         tGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGAt                                                             >  2:4495735/1‑76 (MQ=255)
          gCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATc                                                            >  1:1228762/1‑76 (MQ=255)
           ccAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCg                                                           <  1:1400413/76‑1 (MQ=255)
           ccAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCg                                                           <  1:4212544/76‑1 (MQ=255)
           ccAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCg                                                           <  2:2614418/76‑1 (MQ=255)
           ccAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCg                                                           <  1:3212366/76‑1 (MQ=255)
           ccAGATAAACGGCTTGAAGCCATAATTCAGGAGATGGTTTTTCATTTCCCCCAGCACGCGGTTCTCTTGCGGGTCg                                                           >  2:1504141/1‑76 (MQ=255)
            cAGAAAAAGGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCgg                                                          <  2:4731692/76‑1 (MQ=255)
             aGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGc                                                         >  1:2833423/1‑76 (MQ=255)
             aGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCCGc                                                         >  2:4329923/1‑76 (MQ=255)
             aGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCCGc                                                         >  2:362975/1‑76 (MQ=255)
               aTAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCt                                                       >  1:4367148/1‑76 (MQ=255)
                tAAGCGGCTTAAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTg                                                      <  2:4543423/76‑1 (MQ=255)
                tAAACGGTTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTg                                                      <  2:1644591/76‑1 (MQ=255)
                tAAACGGCTGGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTg                                                      <  2:3006006/76‑1 (MQ=255)
                aaaaCGGCTTAAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTg                                                      <  2:255723/75‑1 (MQ=255)
                 aaaCGTTTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTgg                                                     <  2:1934828/76‑1 (MQ=255)
                 aaaCGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTgg                                                     >  1:1552219/1‑76 (MQ=255)
                   ccGGCTTAAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGaa                                                   <  2:4218559/75‑1 (MQ=255)
                   aCGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGaa                                                   <  1:2960588/76‑1 (MQ=255)
                    cGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAg                                                  <  1:3462426/76‑1 (MQ=255)
                    cGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAg                                                  <  1:1399829/76‑1 (MQ=255)
                      gCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGtt                                                >  1:112981/1‑76 (MQ=255)
                      gCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGtt                                                <  1:2971206/76‑1 (MQ=255)
                        ttGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTg                                              <  2:4883984/76‑1 (MQ=255)
                         tGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGa                                             <  2:3650325/76‑1 (MQ=255)
                         tGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGa                                             >  1:1545845/1‑76 (MQ=255)
                          gAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGAc                                            >  2:1577746/1‑76 (MQ=255)
                           aaGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGAcc                                           <  1:1093281/76‑1 (MQ=255)
                             gCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAg                                         <  2:1732071/76‑1 (MQ=255)
                             gCCATACTACAGGCGATCTTTTTTCATCTCCCACAGCTCGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAg                                         <  2:5023093/76‑1 (MQ=255)
                              ccATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGa                                        >  1:5011130/1‑76 (MQ=255)
                               cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGat                                       >  1:1838689/1‑76 (MQ=255)
                               cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGat                                       <  1:4620297/76‑1 (MQ=255)
                               cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGat                                       >  1:2609439/1‑76 (MQ=255)
                               cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGat                                       <  1:3438398/76‑1 (MQ=255)
                               cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGat                                       <  2:1218918/76‑1 (MQ=255)
                                aTACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGata                                      <  2:3836161/76‑1 (MQ=255)
                                aTACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGata                                      <  2:225078/76‑1 (MQ=255)
                                aTACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGata                                      >  1:616807/1‑76 (MQ=255)
                                  aCTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATc                                    <  1:2206171/76‑1 (MQ=255)
                                     tCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCgtg                                 <  2:4126690/76‑1 (MQ=255)
                                      cAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCgtgg                                <  1:877144/76‑1 (MQ=255)
                                      cAGGCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCgtgg                                >  2:2393266/1‑76 (MQ=255)
                                         gCGATCTTTTTTCGTCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGtggtgg                             <  1:682995/76‑1 (MQ=255)
                                         gCGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGtggtgg                             <  2:4156010/76‑1 (MQ=255)
                                          cGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGa                            <  2:1503093/76‑1 (MQ=255)
                                          cGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGa                            >  2:4951952/1‑76 (MQ=255)
                                          cGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGa                            <  1:1670898/76‑1 (MQ=255)
                                          cGATCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGGGGa                            >  2:101707/1‑76 (MQ=255)
                                            aTCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAAc                          >  1:4839639/1‑76 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         <  1:4279649/76‑1 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         <  2:4234361/76‑1 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         <  1:2053284/76‑1 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         >  2:3212193/1‑76 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         <  2:210912/76‑1 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         <  1:4238723/76‑1 (MQ=255)
                                             tCTTTTTTCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACg                         <  2:1404413/76‑1 (MQ=255)
                                               ttttttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACggg                       <  2:2696313/76‑1 (MQ=255)
                                               ttttttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACggg                       <  1:2659841/76‑1 (MQ=255)
                                               ttttttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGGGGAACggg                       >  2:3200160/1‑76 (MQ=255)
                                                  tttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATa                    <  1:2495106/76‑1 (MQ=255)
                                                  tttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATa                    <  1:2213216/76‑1 (MQ=255)
                                                  tttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATa                    <  1:4760071/76‑1 (MQ=255)
                                                  tttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATa                    >  1:4990749/1‑76 (MQ=255)
                                                   ttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATaa                   <  1:1822447/76‑1 (MQ=255)
                                                   ttCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATaa                   >  1:3664451/1‑76 (MQ=255)
                                                    tCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCCGATAACGTGGTGGAACGGGATAAc                  >  2:4382046/1‑76 (MQ=255)
                                                    tCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAAc                  <  2:2075456/76‑1 (MQ=255)
                                                    tCATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAAc                  >  1:3375320/1‑76 (MQ=255)
                                                     cATCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAAcc                 <  1:4482219/76‑1 (MQ=255)
                                                      aTCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCa                <  2:3890328/76‑1 (MQ=255)
                                                      aTCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCa                <  2:2582498/76‑1 (MQ=255)
                                                      aTCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCa                <  2:4152902/76‑1 (MQ=255)
                                                      aTCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCa                <  2:114544/76‑1 (MQ=255)
                                                      aTCTCCCACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCa                <  2:737896/76‑1 (MQ=255)
                                                          cccACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAAACACttt            >  2:1338961/1‑76 (MQ=255)
                                                           ccACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTc           <  1:3267604/76‑1 (MQ=255)
                                                           ccACAGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTc           >  2:727825/1‑76 (MQ=255)
                                                               aGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCAt       <  2:2836372/76‑1 (MQ=255)
                                                               aGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCAt       <  2:1790292/76‑1 (MQ=255)
                                                               aGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCAt       <  1:2321595/76‑1 (MQ=255)
                                                               aGCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCAt       <  1:1258312/76‑1 (MQ=255)
                                                                gCACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCAtt      <  2:921342/76‑1 (MQ=255)
                                                                 cACGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCATTc     >  2:4302549/1‑76 (MQ=255)
                                                                  aCGCGGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCATTCa    >  1:2898818/1‑76 (MQ=255)
                                                                    gcgGTTCTCTTGCGGATCGGCCTGTAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCATTCAgg  <  2:302187/76‑1 (MQ=255)
                                                                    gcgGTTCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCATTCAgg  <  1:2116239/76‑1 (MQ=255)
                                                                    gcgGTTCTCTCGCGGATCGGCCTGGAAATTTTAGCGGCTATCGTGGTGGGACGGGATAACCCCTTTAGGATTCacg  >  2:2432612/1‑74 (MQ=255)
                                                                        |                                                                       
CCGCCAACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATAACCACTTTCGCATTCAGG  >  REL606/4396825‑4396968

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