Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	3,244,017	G→A	9.7%	Q86Q (CAG→CAA)	exuT →	hexuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	3,244,017	0	G	A	9.7%	285.8 / 15.8	93	Q86Q (CAG→CAA)	exuT	hexuronate transporter
Reads supporting (aligned to +/- strand):  ref base G (47/37);  new base A (6/3);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

ACCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTA  >  AP009048/3243947‑3244083                                                                       |                                                                   aCCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAg                                                                    >  1:1510329/1‑71 (MQ=255)   cTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTa                                                                  <  1:190289/71‑1 (MQ=255)     gACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAtt                                                                >  1:2182533/1‑71 (MQ=255)     gACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                 >  1:3929040/1‑70 (MQ=255)     gACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                 >  1:2267731/1‑70 (MQ=255)       cgcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                              >  1:1384097/1‑71 (MQ=255)       cgcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                              >  1:3299404/1‑71 (MQ=255)        gcgTAACACTGTGGCGGCAGCTGCGCCCACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                              >  1:1558590/1‑70 (MQ=255)        gcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                 <  1:3498315/67‑1 (MQ=255)         cgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCt                                                             <  1:2897286/70‑1 (MQ=255)          gTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCta                                                            >  1:1159015/1‑70 (MQ=255)             acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                         >  1:1765036/1‑70 (MQ=255)             acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                         >  1:1953122/1‑70 (MQ=255)              cacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatca                                                       <  1:8382/71‑1 (MQ=255)              cacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatca                                                       <  1:330893/71‑1 (MQ=255)                 tgtgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCg                                                    >  1:3386397/1‑71 (MQ=255)                   tgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGc                                                   <  1:4176175/70‑1 (MQ=255)                     gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAg                                                 <  1:983918/70‑1 (MQ=255)                      cggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAg                                                 >  1:4529237/1‑69 (MQ=255)                      cggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                <  1:1797487/70‑1 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                             <  1:3648174/71‑1 (MQ=255)                          aGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTAtt                                           >  1:549410/1‑71 (MQ=255)                           gCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTc                                          <  1:2578318/71‑1 (MQ=255)                            cTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatc                                                        <  1:498826/56‑1 (MQ=255)                               cgcCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgct                                      <  1:1772660/71‑1 (MQ=255)                               cgcCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTctgct                                      >  1:1554730/1‑71 (MQ=255)                                gcCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgctg                                     >  1:3679331/1‑71 (MQ=255)                                 ccAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCtgctgc                                    <  1:748722/71‑1 (MQ=255)                                 ccAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATCCtgctgc                                    >  1:1699423/1‑71 (MQ=255)                                  cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcat                                                      >  1:815844/1‑52 (MQ=255)                                  cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                   >  1:3558234/1‑71 (MQ=255)                                  cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                   >  1:4578149/1‑71 (MQ=255)                                  cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTgctgct                                   >  1:2767791/1‑71 (MQ=255)                                   aaCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                   >  1:1713374/1‑70 (MQ=255)                                    aCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTta                                 <  1:3354395/71‑1 (MQ=255)                                     ctctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtat                                <  1:2184373/71‑1 (MQ=255)                                      tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                               >  1:465283/1‑71 (MQ=255)                                       ctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                              >  1:4469410/1‑71 (MQ=255)                                        tGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                               <  1:1988414/69‑1 (MQ=255)                                        tGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                              >  1:1442434/1‑70 (MQ=255)                                        tGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACg                             >  1:115077/1‑71 (MQ=255)                                         gATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACgg                            >  1:2730187/1‑71 (MQ=255)                                          aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCta                                                            >  1:4069065/1‑38 (MQ=255)                                          aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCa                                                  <  1:4546062/48‑1 (MQ=255)                                          aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACgg                            <  1:837087/70‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                   <  1:3437051/62‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGt                           >  1:3971664/1‑70 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                          <  1:2905006/71‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                          <  1:2389897/71‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGt                           >  1:185352/1‑70 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                          <  1:3955196/71‑1 (MQ=255)                                            ggAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctg                                        <  1:2868956/56‑1 (MQ=255)                                            ggAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTta                                 >  1:2904965/1‑63 (MQ=255)                                             gAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCAt                        <  1:489451/71‑1 (MQ=255)                                             gAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCAt                        >  1:2524276/1‑71 (MQ=255)                                              aaGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgct                                      >  1:2682499/1‑56 (MQ=255)                                               agagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGc                      <  1:2536592/71‑1 (MQ=255)                                               agagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCCGCTTATACGGTCATGc                      <  1:1990455/71‑1 (MQ=255)                                               agagTTAAACATCTACACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATg                       >  1:3438182/1‑70 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                >  1:4370456/1‑44 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                              >  1:2024525/1‑62 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                              >  1:1925774/1‑62 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACTGTCATGCa                     <  1:1618197/71‑1 (MQ=255)                                                   ttAAACATCTCCACCCAACAGTATTCCTATATCATCGc                                                   <  1:4443507/38‑1 (MQ=255)                                                   ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctg                                        <  1:4263/49‑1 (MQ=255)                                                   ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACgg                            <  1:333241/61‑1 (MQ=255)                                                   ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACgg                            <  1:4033590/61‑1 (MQ=255)                                                   ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCAt                        <  1:3432789/65‑1 (MQ=255)                                                   ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAAcc                  <  1:2761103/71‑1 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCa                                                  >  1:697844/1‑37 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                             <  1:2083264/42‑1 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACg                             >  1:1733541/1‑58 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:1657388/1‑71 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:1983414/1‑71 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:254069/1‑71 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:2612723/1‑71 (MQ=255)                                                       aCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:3234187/1‑69 (MQ=255)                                                       aCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:2116465/1‑69 (MQ=255)                                                       aCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGt               <  1:1084478/70‑1 (MQ=255)                                                         aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCaa                    <  1:457215/63‑1 (MQ=255)                                                         aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                >  1:4206573/1‑67 (MQ=255)                                                         aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGc            <  1:2566494/71‑1 (MQ=255)                                                         aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACAGGTa              >  1:2978696/1‑69 (MQ=255)                                                          tctcCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgct                                      >  1:1731284/1‑44 (MQ=255)                                                            tcCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAg             >  1:2435735/1‑67 (MQ=255)                                                              cACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt        >  1:2861598/1‑70 (MQ=255)                                                              cACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                          <  1:1564177/52‑1 (MQ=255)                                                               aCCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTa      <  1:4409861/71‑1 (MQ=255)                                                               aCCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTa      >  1:2925107/1‑71 (MQ=255)                                                                 ccAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCtt                                  >  1:3586831/1‑41 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa  >  1:2590367/1‑71 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa  >  1:4238494/1‑71 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa  <  1:313024/71‑1 (MQ=255)                                                                      aGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGtt       >  1:171354/1‑63 (MQ=255)                                                                       |                                                                   ACCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTA  >  AP009048/3243947‑3244083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑