Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	454,680	A→T	20.0%	R135R (CGT→CGA)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	454,680	0	A	T	20.0%	172.7 / 39.7	85	R135R (CGT→CGA)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base T (9/8);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCG  >  Exported/454612‑454749                                                                     |                                                                      caGGCGATAGGTGCGATAATAGACCCCACTTTACCCGCCGCATACATCAGCGAGAAACCGCCATCacgacg                                                                     <  1:1461168/71‑1 (MQ=255)       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:4880521/1‑70 (MQ=255)                                tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                       >  1:261339/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)                                                                     |                                                                      CAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCG  >  Exported/454612‑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: ‑