Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	485,028	T→G	9.9%	G278G (GGT→GGG)	tolB →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	485,028	0	T	G	9.9%	165.0 / 15.1	71	G278G (GGT→GGG)	tolB	periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (3/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

CAGCAAACTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCGGATAGCC  >  Exported/484959‑485092                                                                      |                                                                 cagcaAACTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGt                                                                  >  1:29784/1‑70 (MQ=255)  agcaAACTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCa                                                                <  1:1633100/71‑1 (MQ=255)    caAACTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGa                                                              <  1:2987171/71‑1 (MQ=255)    caAACTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGa                                                              <  1:40925/71‑1 (MQ=255)      aaCTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATc                                                            >  1:1849025/1‑71 (MQ=255)        cTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCg                                                          <  1:2429843/71‑1 (MQ=255)        cTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCg                                                          >  1:3703464/1‑71 (MQ=255)         tGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATc                                                            >  1:1958924/1‑68 (MQ=255)           gCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGc                                                         >  1:2440316/1‑69 (MQ=255)              ttCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGGCAGATCCg                                                          <  1:2462497/65‑1 (MQ=255)              ttCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGGCAGATCCg                                                          <  1:1987576/65‑1 (MQ=255)               tCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGt                                                    >  1:2041832/1‑70 (MQ=255)                cGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAg                                                               >  1:3382764/1‑58 (MQ=255)                cGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAg                                                               >  1:128846/1‑58 (MQ=255)                cGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAg                                                               >  1:1870266/1‑58 (MQ=255)                cGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGGCAGATCCGcc                                                        >  1:407975/1‑65 (MQ=255)                   cTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTg                                               <  1:1858693/71‑1 (MQ=255)                    ttGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAg                                                               <  1:2008353/54‑1 (MQ=255)                     tGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGAt                                             >  1:751003/1‑71 (MQ=255)                      gTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATg                                            >  1:3437555/1‑71 (MQ=255)                       tCGAAAACCGGTAGTCTGAACCTGTACGTAATTGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATgg                                           <  1:53932/71‑1 (MQ=255)                       tCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCa                                                                <  1:2854838/50‑1 (MQ=255)                           aaaCCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCg                                        >  1:3215278/1‑70 (MQ=255)                            aaCCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCTGATCCGCCAGGTGACTGATGGTCGCa                                      <  1:2874313/71‑1 (MQ=255)                            aaCCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCa                                      <  1:1076314/71‑1 (MQ=255)                            aaCCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCa                                      <  1:3258230/71‑1 (MQ=255)                            aaCCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCa                                      <  1:1029531/71‑1 (MQ=255)                              ccGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGt                                    <  1:2314820/71‑1 (MQ=255)                               cGGTAGTCTGAACCTGTACGTAATGGATTTGGGTTCTGGTCAGATCCGCCAGGTGACTGa                                              <  1:159940/60‑1 (MQ=255)                                 gTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGGCAGATCCg                                                          <  1:425666/46‑1 (MQ=255)                                   aGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAg                                                               >  1:3428056/1‑39 (MQ=255)                                   aGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAg                                                               >  1:2621795/1‑39 (MQ=255)                                       tGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATac                            >  1:2489845/1‑70 (MQ=255)                                        gAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTg                                                   <  1:3023350/46‑1 (MQ=255)                                        gAACCTGTACGTAATGGATTTGGCTTCTGGGCAGATCCg                                                          <  1:632593/39‑1 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGAc                                                 >  1:3327855/1‑47 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGAc                                                 >  1:2970251/1‑47 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGAc                                                 >  1:632062/1‑47 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGAc                                                 >  1:2294096/1‑47 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccga                         <  1:3432689/71‑1 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:689764/70‑1 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          >  1:2541741/1‑70 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:2865554/70‑1 (MQ=255)                                         aaCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:1723030/70‑1 (MQ=255)                                          aCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccgaa                        >  1:253307/1‑71 (MQ=255)                                          aCCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccgaa                        >  1:3351661/1‑71 (MQ=255)                                            cTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccgaacc                      >  1:1855269/1‑71 (MQ=255)                                              gTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCa                                                       >  1:956047/1‑36 (MQ=255)                                              gTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccgaaccga                    >  1:2101180/1‑71 (MQ=255)                                                 cGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCt                 <  1:2248863/71‑1 (MQ=255)                                                   tAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCg                                        <  1:1694324/46‑1 (MQ=255)                                                      tGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTc            <  1:2694509/71‑1 (MQ=255)                                                      tGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTc            >  1:1145882/1‑71 (MQ=255)                                                      tGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTc            >  1:527984/1‑71 (MQ=255)                                                      tGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTc            <  1:3670401/71‑1 (MQ=255)                                                      tGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTc            >  1:2343197/1‑71 (MQ=255)                                                      tGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTc            >  1:3564168/1‑71 (MQ=255)                                                        gATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTc                                         >  1:1321579/1‑40 (MQ=255)                                                        gATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTcc           <  1:3642929/70‑1 (MQ=255)                                                        gATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTcc           <  1:1342349/70‑1 (MQ=255)                                                         aTTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:527784/54‑1 (MQ=255)                                                         aTTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:2648697/54‑1 (MQ=255)                                                         aTTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:1348483/54‑1 (MQ=255)                                                         aTTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccg                          <  1:507337/54‑1 (MQ=255)                                                         aTTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCg         <  1:2321601/71‑1 (MQ=255)                                                          tttGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCg         <  1:1717416/70‑1 (MQ=255)                                                            tGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATaccgaa                        >  1:1717995/1‑53 (MQ=255)                                                            tGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCGGa       >  1:990719/1‑70 (MQ=255)                                                             ggCTTCTGGGCAGATCCGCCAGGTGACTGATGGTCGCAGTAAc                                 >  1:1467757/1‑43 (MQ=255)                                                             ggCTTCTGGGCAGATCCGCCAGGTGACTGATGGTCGCAGTAAc                                 >  1:1946131/1‑43 (MQ=255)                                                                  cTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCGGATAGcc  <  1:2300327/69‑1 (MQ=255)                                                                  cTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCGGATAGc   >  1:2150512/1‑68 (MQ=255)                                                                      |                                                                 CAGCAAACTGGCATTCGCCTTGTCGAAAACCGGTAGTCTGAACCTGTACGTAATGGATTTGGCTTCTGGTCAGATCCGCCAGGTGACTGATGGTCGCAGTAACAATACCGAACCGACCTGGTTCCCGGATAGCC  >  Exported/484959‑485092

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