Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	67,598	C→T	12.3%	G247D (GGC→GAC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	67,598	0	C	T	12.3%	182.1 / 14.5	65	G247D (GGC→GAC)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (4/4);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

TTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCG  >  W3110S.gb/67531‑67666                                                                    |                                                                     ttGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGt                                                                   >  1:77962/1‑71 (MQ=255)   gTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTg                                                                  <  1:472557/70‑1 (MQ=255)    ttCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGa                                                                <  1:3932722/71‑1 (MQ=255)       cAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGAtt                                                              <  1:2365930/70‑1 (MQ=255)        aGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtg                                                            <  1:3248534/71‑1 (MQ=255)        aGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtg                                                            <  1:277077/71‑1 (MQ=255)        aGGAAACGCTTCATCCCCAGCTCAATACGCGCAGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtg                                                            <  1:4734526/71‑1 (MQ=255)         ggAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtg                                                            <  1:2996184/70‑1 (MQ=255)          gAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtgtg                                                          >  1:3953515/1‑71 (MQ=255)          gAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCAGTGGATTtgtg                                                          >  1:4716129/1‑71 (MQ=255)            aaCGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtgtgtg                                                        <  1:147401/71‑1 (MQ=255)            aaCGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtgtgtg                                                        <  1:4980520/71‑1 (MQ=255)            aaCGCTTCATCCCCAGCTCAATACGCGCAGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtgtgtg                                                        <  1:3200196/71‑1 (MQ=255)                cTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAg                                                    <  1:1111958/71‑1 (MQ=255)                 ttCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAgg                                                   >  1:3874627/1‑71 (MQ=255)                  tCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGc                                                  >  1:3559819/1‑71 (MQ=255)                  tCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGc                                                  >  1:5123322/1‑71 (MQ=255)                   cATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCg                                                 >  1:84678/1‑71 (MQ=255)                   cATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCg                                                 >  1:4102502/1‑71 (MQ=255)                       cccAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCa                                              <  1:4546100/70‑1 (MQ=255)                          aGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtg                                                            <  1:4532431/53‑1 (MQ=255)                          aGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGGGTGGCAGGCGt                                                <  1:184867/65‑1 (MQ=255)                              cAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTag                                      >  1:4282675/1‑71 (MQ=255)                               aaTACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtgtgt                                                         <  1:1036782/51‑1 (MQ=255)                                 tACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGTCGTGGATTtg                                                            <  1:2863306/46‑1 (MQ=255)                                 tACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGTCGTGGATTtg                                                            <  1:284902/46‑1 (MQ=255)                                 tACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTgg                                                                 >  1:3703052/1‑41 (MQ=255)                                 tACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGa                                                                >  1:4295547/1‑42 (MQ=255)                                 tACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGa                                                                >  1:1647501/1‑42 (MQ=255)                                 tACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGAtt                                                              <  1:4358194/44‑1 (MQ=255)                                  aCGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTagcagc                                  >  1:63992/1‑71 (MQ=255)                                  aCGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTagcagc                                  >  1:4039732/1‑71 (MQ=255)                                  aCGCGCCGCTTACAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTagca                                    >  1:322932/1‑69 (MQ=255)                                   cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTagcagc                                  >  1:461196/1‑70 (MQ=255)                                   cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCt                                 >  1:835141/1‑71 (MQ=255)                                   cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCt                                 >  1:4216014/1‑71 (MQ=255)                                      gccgcTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTtg                                                            <  1:4889270/41‑1 (MQ=255)                                      gccgcTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTc                              <  1:3009984/71‑1 (MQ=255)                                       ccgcTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCg                             <  1:1098034/71‑1 (MQ=255)                                         gcTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATg                                            <  1:1192582/54‑1 (MQ=255)                                            tCCAGCACGTTCTGTCGTTTTTTGTCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                         >  1:5098996/1‑70 (MQ=255)                                            tCCAGCACGTTCTGTCGTTTTTTGTCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                         >  1:2001236/1‑70 (MQ=255)                                            tCCAGCACGTTCTGTCGTTTTTTGTCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                         >  1:2890916/1‑70 (MQ=255)                                            tCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTagcagc                                  >  1:5098701/1‑61 (MQ=255)                                            tCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                         >  1:413909/1‑70 (MQ=255)                                            tCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                         >  1:28319/1‑70 (MQ=255)                                              cAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCa                       <  1:2175106/70‑1 (MQ=255)                                               aGCACGTTCTGTCGTTTTTTGTCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATc                     <  1:807392/71‑1 (MQ=255)                                               aGCACGTTCTGTCGTTTTTTGTCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATc                     <  1:2489281/71‑1 (MQ=255)                                                gCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                         <  1:2641182/66‑1 (MQ=255)                                                 cACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCTGGCGTCATGGTGTAGCAGCTTTCGTACTCATCg                    <  1:944868/70‑1 (MQ=255)                                                        tGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc                        <  1:2295423/59‑1 (MQ=255)                                                        tGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc                        <  1:123869/59‑1 (MQ=255)                                                        tGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc                        <  1:3992373/59‑1 (MQ=255)                                                         gTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCt                                 >  1:2290064/1‑49 (MQ=255)                                                          tCGTTTTTTGTCGTGGATTTGTGTGGCAGGCGTCATGgt                                          >  1:3964070/1‑39 (MQ=255)                                                            gTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGtt         >  1:3304174/1‑70 (MQ=255)                                                             ttttttGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTa        <  1:4754894/70‑1 (MQ=255)                                                              tttttGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTaa       >  1:3551680/1‑70 (MQ=255)                                                               ttttGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAAc      <  1:1147194/70‑1 (MQ=255)                                                               ttttGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACa     >  1:1191001/1‑71 (MQ=255)                                                               ttttGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACa     >  1:41824/1‑71 (MQ=255)                                                               ttttGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACa     >  1:1133662/1‑71 (MQ=255)                                                                   gCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCg  <  1:4852675/70‑1 (MQ=255)                                                                   gCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCg  <  1:2990435/70‑1 (MQ=255)                                                                    |                                                                     TTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCG  >  W3110S.gb/67531‑67666

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: ‑