Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,468,889	C→T	6.5%	E945K (GAA→AAA)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,468,889	0	C	T	6.5%	194.0 / 3.3	61	E945K (GAA→AAA)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGA  >  Exported/1468825‑1468958                                                                 |                                                                      cgcAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTTcgggtc                                                                 <  1:2443621/71‑1 (MQ=255)    aGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCcgggtcggg                                                              <  1:1788999/71‑1 (MQ=255)           taAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGa                                                       >  1:625085/1‑71 (MQ=255)            aaGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGa                                                       <  1:1131064/70‑1 (MQ=255)            aaGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGAt                                                      >  1:2298319/1‑71 (MQ=255)             aGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGAt                                                      >  1:735990/1‑70 (MQ=255)              gATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCa                                                    <  1:2970723/71‑1 (MQ=255)                  aCCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGgtgt                                                           >  1:1491961/1‑60 (MQ=255)                  aCCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGAt                                                >  1:3346566/1‑71 (MQ=255)                      ccAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCgg                                            >  1:2476549/1‑71 (MQ=255)                      ccAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCgg                                            >  1:239840/1‑71 (MQ=255)                        aaCTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCcgggtcgggt                                                             <  1:1060700/52‑1 (MQ=255)                        aaCTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCg                                          <  1:1980325/71‑1 (MQ=255)                        aaCTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCg                                          <  1:3165158/71‑1 (MQ=255)                         aCTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCcgggtcggg                                                              >  1:3462874/1‑50 (MQ=255)                          cTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGt                                         <  1:2374538/70‑1 (MQ=255)                          cTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGt                                         <  1:16863/70‑1 (MQ=255)                          cTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGt                                         <  1:318715/70‑1 (MQ=255)                           tttGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCg                                          <  1:42450/68‑1 (MQ=255)                               gCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGAt                                    >  1:3391349/1‑70 (MQ=255)                                cATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCcgggtc                                                                 <  1:2079595/40‑1 (MQ=255)                                cATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCa                                  <  1:2130107/71‑1 (MQ=255)                                  ttCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGc                                           <  1:3620741/60‑1 (MQ=255)                                  ttCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAg                                 >  1:707176/1‑70 (MQ=255)                                      acacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGTGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCg                            <  1:552856/71‑1 (MQ=255)                                       cacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGc                                           >  1:435044/1‑55 (MQ=255)                                       cacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATc                                   <  1:3695395/63‑1 (MQ=255)                                       cacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCgg                           >  1:3374354/1‑71 (MQ=255)                                       cacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCgg                           >  1:1620023/1‑71 (MQ=255)                                       cacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCgg                           >  1:1667569/1‑71 (MQ=255)                                       cacCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCgg                           >  1:14706/1‑71 (MQ=255)                                         ccGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGa                         <  1:2527114/71‑1 (MQ=255)                                         ccGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTCGGATCAGGATCCGGCGTGGGATCAGGATCGGGa                         <  1:366209/71‑1 (MQ=255)                                          cGGCTGGTAAGCAGGCGTAGGTTTCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAg                        >  1:2738809/1‑71 (MQ=255)                                          cGGCTGGTAAGCAGGCGTAGGTTTCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAg                        <  1:2787961/71‑1 (MQ=255)                                           ggCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGt                       <  1:928304/71‑1 (MQ=255)                                           ggCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGt                       <  1:3307274/71‑1 (MQ=255)                                            gCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTc                      <  1:2810424/71‑1 (MQ=255)                                             cTGGTAAGCAGGCGTAGGTTTCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCg                     >  1:3553062/1‑71 (MQ=255)                                             cTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGAt                              <  1:856691/62‑1 (MQ=255)                                              tGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATc                                   <  1:2436731/56‑1 (MQ=255)                                              tGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATc                                   <  1:2528575/56‑1 (MQ=255)                                                gTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATcc                                              >  1:3457646/1‑43 (MQ=255)                                                gTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATc                                   >  1:2196087/1‑54 (MQ=255)                                                gTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATc                                   >  1:537831/1‑54 (MQ=255)                                                 tAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTc                 <  1:502635/71‑1 (MQ=255)                                                  aaGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCa                                                    >  1:1604908/1‑35 (MQ=255)                                                  aaGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCt                >  1:1298417/1‑71 (MQ=255)                                                   aGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCt                <  1:3526444/70‑1 (MQ=255)                                                   aGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTg               >  1:1857640/1‑71 (MQ=255)                                                    gCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCt                <  1:594359/69‑1 (MQ=255)                                                      aGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGAt            <  1:122968/71‑1 (MQ=255)                                                       ggCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGa             >  1:2756146/1‑69 (MQ=255)                                                        gCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCgg                                            >  1:1810324/1‑37 (MQ=255)                                                        gCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCgg                                            >  1:542443/1‑37 (MQ=255)                                                          gTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCggag                         >  1:3632564/1‑52 (MQ=255)                                                              gTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGc    <  1:3733043/71‑1 (MQ=255)                                                               ttCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGc    >  1:2630496/1‑70 (MQ=255)                                                                tCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGa  >  1:1587617/1‑71 (MQ=255)                                                                tCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGa  >  1:3605637/1‑71 (MQ=255)                                                                tCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGa  >  1:973039/1‑71 (MQ=255)                                                                 ccGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGa  >  1:643112/1‑70 (MQ=255)                                                                 ccGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGa  >  1:820947/1‑70 (MQ=255)                                                                 |                                                                      CGCAGGTTATTAAGATAACCGCCAACTTTGGCATTCAACACCGGCTGGTAAGCAGGCGTAGGTTCCGGGTCGGGTGTTGGATCAGGATCCGGCGTGGGATCAGGATCGGGAGTCGGGTCTGGATCAGGTGGCGA  >  Exported/1468825‑1468958

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