Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	454,683	C→A	17.1%	P134P (CCG→CCT) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	454,683	0	C	A	17.1%	268.2 / 29.3	99	P134P (CCG→CCT) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (38/44);  new base A (9/8);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCG  >  Exported/454618‑454749                                                                  |                                                                   aTAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATc                                                               <  1:504591/71‑1 (MQ=255)     gTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGtt                                                           >  1:2602373/1‑71 (MQ=255)           taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     >  1:4584737/1‑71 (MQ=255)           taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     <  1:1288495/71‑1 (MQ=255)           taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     >  1:4852968/1‑71 (MQ=255)            aataGAACCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGCCGCGGATCGGt                                                            >  1:733942/1‑63 (MQ=255)             ataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCAt                                                    <  1:5206795/70‑1 (MQ=255)                gACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAg                                                >  1:1597261/1‑71 (MQ=255)                   cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGc                                                        <  1:2548731/60‑1 (MQ=255)                   cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                                               <  1:3234855/69‑1 (MQ=255)                    ccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCgg                                                                  >  1:4515173/1‑49 (MQ=255)                    ccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTc                                             <  1:5219523/70‑1 (MQ=255)                     cACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGAt                                                                >  1:5310487/1‑50 (MQ=255)                       cGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGcc                                          >  1:347294/1‑70 (MQ=255)                       cGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCg                                         <  1:5280789/71‑1 (MQ=255)                        gTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGtt                                                           <  1:1217606/52‑1 (MQ=255)                         ttACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     <  1:4590272/57‑1 (MQ=255)                         ttACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       <  1:2071261/71‑1 (MQ=255)                          tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTg                                                          >  1:98821/1‑51 (MQ=255)                          tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                                                   >  1:309591/1‑58 (MQ=255)                          tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       >  1:261339/1‑70 (MQ=255)                          tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       >  1:4880521/1‑70 (MQ=255)                           gcccGCAGCATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGc                                      >  1:1873060/2‑70 (MQ=255)                            cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                     <  1:3495535/70‑1 (MQ=255)                            cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                     <  1:3502401/70‑1 (MQ=255)                            cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                     <  1:2095070/70‑1 (MQ=255)                            cccGCAGCATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCa                                     <  1:4836805/70‑1 (MQ=255)                             ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     <  1:3817972/53‑1 (MQ=255)                                 agcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAAc                                  <  1:5091729/67‑1 (MQ=255)                                 agcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAg                                >  1:4744271/2‑70 (MQ=255)                                   cATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCt                              >  1:4100189/1‑70 (MQ=255)                                    aTACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtga                            >  1:49635/1‑69 (MQ=255)                                     tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     <  1:4205409/45‑1 (MQ=255)                                     tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                     <  1:2830171/45‑1 (MQ=255)                                     tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       <  1:5309415/59‑1 (MQ=255)                                     tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       <  1:3852371/59‑1 (MQ=255)                                     tACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtgat                           >  1:296372/1‑68 (MQ=255)                                     tACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtgat                           >  1:3753874/1‑68 (MQ=255)                                       cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                     >  1:1126806/1‑59 (MQ=255)                                       cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                         >  1:4004165/1‑71 (MQ=255)                                       cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                         >  1:1176694/1‑71 (MQ=255)                                        aTCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                        <  1:3866004/71‑1 (MQ=255)                                        aTCAGCGAGAAAACGCCATCACGACGCGGATCGGTTGGCTCa                                                     <  1:5260130/42‑1 (MQ=255)                                         tCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa                                     >  1:5491708/1‑57 (MQ=255)                                          cAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGttt                       >  1:4326694/1‑69 (MQ=255)                                           aGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                                                   >  1:1199174/1‑41 (MQ=255)                                           aGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                                                   >  1:3734168/1‑41 (MQ=255)                                           aGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAg                                                >  1:593407/1‑44 (MQ=255)                                           aGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGa                     <  1:1171842/71‑1 (MQ=255)                                            gCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                                               >  1:4334390/1‑44 (MQ=255)                                              gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                   >  1:4849648/1‑70 (MQ=255)                                              gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                   >  1:4574586/1‑70 (MQ=255)                                              gagaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                   >  1:2078218/1‑70 (MQ=255)                                              gaaaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAttt                  <  1:1086908/68‑1 (MQ=255)                                               agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAg                                                <  1:1334661/40‑1 (MQ=255)                                               agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 <  1:3817507/71‑1 (MQ=255)                                               agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 <  1:742179/71‑1 (MQ=255)                                               agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 <  1:2163150/71‑1 (MQ=255)                                               agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 <  1:1109709/71‑1 (MQ=255)                                               agaAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 <  1:1105808/71‑1 (MQ=255)                                               aaaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtttg                 <  1:1687386/69‑2 (MQ=255)                                                aaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtttga                <  1:3393136/70‑3 (MQ=255)                                                aaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtttga                <  1:1109761/70‑3 (MQ=255)                                                aaaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGAtttg                 >  1:674426/2‑69 (MQ=255)                                                 aaaCCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtga                            >  1:1477802/1‑56 (MQ=255)                                                 aaaCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 >  1:2577328/1‑69 (MQ=255)                                                  aaCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa               <  1:1567786/70‑1 (MQ=255)                                                  aaCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAAc              <  1:5627258/71‑1 (MQ=255)                                                  aaCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAAc              <  1:5510105/71‑1 (MQ=255)                                                   aCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa               >  1:3771252/1‑69 (MQ=255)                                                   aCCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa               >  1:2965054/1‑69 (MQ=255)                                                    ccGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCTGATGTTTGATTTgaataa            <  1:3078987/71‑4 (MQ=255)                                                      gCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAttt                  >  1:1372152/1‑63 (MQ=255)                                                       ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTAc                           >  1:1110962/1‑53 (MQ=255)                                                       ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTAc                           >  1:5176542/1‑53 (MQ=255)                                                       ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTAc                           >  1:5221943/1‑53 (MQ=255)                                                       ccATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAtt                   <  1:3610407/61‑1 (MQ=255)                                                        cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       >  1:3695262/1‑40 (MQ=255)                                                        cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa               >  1:2182490/1‑64 (MQ=255)                                                        cATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaaa               >  1:696834/1‑64 (MQ=255)                                                          tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGAc                                  <  1:4230463/43‑1 (MQ=255)                                                          tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                        <  1:4060548/53‑1 (MQ=255)                                                          tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                        <  1:3970061/53‑1 (MQ=255)                                                          tCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTa                 <  1:3263689/60‑1 (MQ=255)                                                           cACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGAc                                  <  1:3938897/42‑1 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       <  1:2997829/36‑1 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGa                                   >  1:2797306/1‑40 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                >  1:2279261/1‑43 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAg                                >  1:2915267/1‑43 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                              <  1:321950/45‑1 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                              <  1:3711023/45‑1 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                              <  1:3996830/45‑1 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                              <  1:2100564/45‑1 (MQ=255)                                                            acgacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGc    >  1:3275336/1‑71 (MQ=255)                                                              gacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg  >  1:3473341/1‑71 (MQ=255)                                                              gacgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg  <  1:1205343/71‑1 (MQ=255)                                                               acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                         <  1:3498253/47‑1 (MQ=255)                                                               acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt                         <  1:2854758/47‑1 (MQ=255)                                                               acgCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaa                <  1:2667699/56‑1 (MQ=255)                                                                  |                                                                   ATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCG  >  Exported/454618‑454749

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