Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,176,730	T→C	Y349C (TAC→TGC)	uxaC ←	glucuronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,176,730	0	T	C	98.7%	295.7 / ‑4.7	79	Y349C (TAC→TGC)	uxaC	glucuronate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (34/44);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TTACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGTAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGC  >  REL606/3176682‑3176777                                                 |                                                ttACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGcagcag                                               >  1:5597554/1‑51 (MQ=255) ttACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGcagcag                                               <  1:2513683/51‑1 (MQ=255) ttACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGcagcag                                               <  1:1619420/51‑1 (MQ=255)  tACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCaga                                              <  1:5095983/51‑1 (MQ=255)  tACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCaga                                              <  1:2637289/51‑1 (MQ=255)   aCCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCagag                                             <  1:6428291/51‑1 (MQ=255)    ccGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAgaga                                            <  1:2914545/51‑1 (MQ=255)     cGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGaa                                           <  1:3169263/51‑1 (MQ=255)      gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:6452193/51‑1 (MQ=255)      gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:2337881/51‑1 (MQ=255)      gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:1765262/51‑1 (MQ=255)      gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          >  1:3980043/1‑51 (MQ=255)      gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:5166842/51‑1 (MQ=255)        tCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATgg                                        <  1:3389183/51‑1 (MQ=255)        tCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATgg                                        >  1:2222926/1‑51 (MQ=255)        tCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATgg                                        >  1:4082928/1‑51 (MQ=255)         cATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGt                                       >  1:1171104/1‑51 (MQ=255)          aTGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtt                                      <  1:5993719/51‑1 (MQ=255)          aTGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtt                                      >  1:3892888/1‑51 (MQ=255)          aTGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtt                                      <  1:6815832/51‑1 (MQ=255)           tGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGTAGAGAATGGttt                                     <  1:6096369/51‑1 (MQ=255)           tGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGttt                                     <  1:4500008/51‑1 (MQ=255)           tGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGttt                                     >  1:1338028/1‑51 (MQ=255)            ggTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtttt                                    >  1:2102316/1‑51 (MQ=255)            ggTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtttt                                    <  1:66350/51‑1 (MQ=255)             gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   <  1:6722022/51‑1 (MQ=255)             gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   >  1:1412901/1‑51 (MQ=255)             gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   <  1:5284360/51‑1 (MQ=255)             gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   <  1:3229702/51‑1 (MQ=255)             gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   >  1:3483954/1‑51 (MQ=255)                cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:6653531/1‑51 (MQ=255)                cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:2361622/1‑51 (MQ=255)                cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:1309288/1‑51 (MQ=255)                 ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCTGAGAATGGTTTTCGgca                               <  1:5545965/51‑1 (MQ=38)                 ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:4232747/1‑51 (MQ=255)                 ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:6668000/1‑51 (MQ=255)                   aGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagc                             <  1:5471092/51‑1 (MQ=255)                    gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagca                            <  1:4745076/51‑1 (MQ=255)                    gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagca                            <  1:3298278/51‑1 (MQ=255)                    gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagca                            <  1:4496063/51‑1 (MQ=255)                    gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:5007329/1‑48 (MQ=255)                     cACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGcagcag                           >  1:3188497/1‑51 (MQ=255)                      aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          <  1:2764111/51‑1 (MQ=255)                      aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          <  1:6474601/51‑1 (MQ=255)                      aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          <  1:5448173/51‑1 (MQ=255)                      aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:1042065/1‑51 (MQ=255)                      aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:4408224/1‑51 (MQ=255)                       cTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGt                         >  1:2244285/1‑51 (MQ=255)                        ttCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:5339435/1‑49 (MQ=255)                        ttCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGtt                        <  1:2257041/51‑1 (MQ=255)                         tCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGcagcag                           <  1:582316/47‑1 (MQ=255)                         tCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGttt                       >  1:6874736/1‑51 (MQ=255)                          cGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:2566394/1‑47 (MQ=255)                          cGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGtttt                      >  1:1855230/1‑51 (MQ=255)                           gTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttc                     <  1:6594317/51‑1 (MQ=255)                           gTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttc                     <  1:5154441/51‑1 (MQ=255)                            ttATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttct                    >  1:4255353/1‑51 (MQ=255)                             tATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttctt                   >  1:6323834/1‑51 (MQ=255)                                cgcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttct                    <  1:5197235/47‑1 (MQ=255)                                cgcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCAt                <  1:3219023/51‑1 (MQ=255)                                 gcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCAtt               >  1:3705320/1‑51 (MQ=255)                                  cgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTc              <  1:2920526/51‑1 (MQ=255)                                  cgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTc              >  1:3286005/1‑51 (MQ=255)                                   gcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCt             <  1:3847803/51‑1 (MQ=255)                                   gcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCt             <  1:1164899/51‑1 (MQ=255)                                    cgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtg            >  1:3892246/1‑51 (MQ=255)                                    cgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtg            <  1:308756/51‑1 (MQ=255)                                    cgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtg            <  1:2704524/51‑1 (MQ=255)                                     gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           >  1:4550342/1‑51 (MQ=255)                                     gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           <  1:3035602/51‑1 (MQ=255)                                     gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           <  1:3592645/51‑1 (MQ=255)                                     gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           >  1:391559/1‑51 (MQ=255)                                      ggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgct          >  1:1664651/1‑51 (MQ=255)                                       gTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgctt         <  1:101007/51‑1 (MQ=255)                                       gTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgctt         <  1:5712666/51‑1 (MQ=255)                                        ttCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttg        <  1:153430/51‑1 (MQ=255)                                         tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       >  1:3324381/1‑51 (MQ=255)                                         tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       >  1:404963/1‑51 (MQ=255)                                         tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       <  1:6388064/51‑1 (MQ=255)                                         tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       <  1:2038108/51‑1 (MQ=255)                                          cAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgct      >  1:7281051/1‑51 (MQ=255)                                           aGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTc     <  1:234703/51‑1 (MQ=255)                                             gcagcaGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAg   >  1:2727651/1‑51 (MQ=255)                                              cagcagAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGc  >  1:4273599/1‑51 (MQ=255)                                              cagcagAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGc  >  1:4274854/1‑51 (MQ=255)                                                 |                                                TTACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGTAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGC  >  REL606/3176682‑3176777

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑