Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	313,852	A→G	19.3%	R79R (CGT→CGC)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	313,852	0	A	G	19.3%	168.2 / 33.7	83	R79R (CGT→CGC)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (10/6);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGA  >  Exported/313783‑313921                                                                      |                                                                      cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc                                                                      <  1:2847420/71‑1 (MQ=255) cTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTAc                                                                      <  1:2489528/71‑1 (MQ=255)   gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc                                                                    >  1:5075756/1‑71 (MQ=255)   gATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGc                                                                    >  1:1979692/1‑71 (MQ=255)          cGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGa                                                             >  1:895848/1‑71 (MQ=255)           gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                            <  1:261144/71‑1 (MQ=255)           gCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                            <  1:772659/71‑1 (MQ=255)             aGGATCAATCTGAAAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGTTTCAGGATaa                                                          <  1:2256437/71‑1 (MQ=255)                 tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATccc                                                      >  1:3630095/1‑71 (MQ=255)                 tCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATcc                                                       >  1:88970/1‑70 (MQ=255)                  cAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCa                                                                >  1:637433/1‑60 (MQ=255)                       tGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAAc                                                >  1:5186021/1‑71 (MQ=255)                         aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTa                                                  <  1:472724/67‑1 (MQ=255)                         aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACtt                                              <  1:4302614/71‑1 (MQ=255)                         aTAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACt                                               <  1:4365912/70‑1 (MQ=255)                            gagagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAg                                           <  1:3355624/71‑1 (MQ=255)                             agagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAgg                                          <  1:2898845/71‑1 (MQ=255)                              gagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTa                                                                       >  1:2676447/1‑41 (MQ=255)                              gagagaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTa                                                                       >  1:409075/1‑41 (MQ=255)                                 agaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAAt                                                         >  1:4562670/1‑52 (MQ=255)                                 agaCACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                    >  1:3336319/1‑57 (MQ=255)                                   acacCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                    <  1:5298/71‑1 (MQ=255)                                      ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGAt                                                            <  1:688721/44‑1 (MQ=255)                                      ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                             <  1:320484/59‑1 (MQ=255)                                      ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                             <  1:755252/59‑1 (MQ=255)                                      ccTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACg                                       <  1:3472635/65‑1 (MQ=255)                                        tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                          >  1:2894055/1‑44 (MQ=255)                                        tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGa                                         >  1:1226219/1‑61 (MQ=255)                                        tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTc                                 <  1:268844/69‑1 (MQ=255)                                        tGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTc                                 <  1:390887/69‑1 (MQ=255)                                         gCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGaa                                     >  1:3313676/1‑64 (MQ=255)                                         gCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGaa                                     >  1:5434545/1‑64 (MQ=255)                                          cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                          >  1:403855/1‑42 (MQ=255)                                          cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                          >  1:5225794/1‑42 (MQ=255)                                          cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                               >  1:4422835/1‑69 (MQ=255)                                          cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                               >  1:1906437/1‑69 (MQ=255)                                          cTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTg                               >  1:3209556/1‑69 (MQ=255)                                             cGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATc                          >  1:4141687/1‑71 (MQ=255)                                               acatCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACt                                   <  1:4917131/57‑1 (MQ=255)                                               acatCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACt                                   <  1:5261606/57‑1 (MQ=255)                                               aTGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAg                                           <  1:3123355/52‑1 (MQ=255)                                                tGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                       <  1:4271291/71‑1 (MQ=255)                                                 gTCGCTACCTTCTTTGAAATTACGCTTCAGGATaa                                                          >  1:1187704/1‑35 (MQ=255)                                                 gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                    >  1:4953699/1‑41 (MQ=255)                                                 gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGc                                                    >  1:4344865/1‑41 (MQ=255)                                                 gTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCgg                       >  1:4521836/1‑70 (MQ=255)                                                 atCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGc                              <  1:3393494/62‑1 (MQ=255)                                                  tCGCTGCCTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACgat                      >  1:2845283/1‑68 (MQ=255)                                                    gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                  <  1:2650266/56‑1 (MQ=255)                                                    gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAgg                   >  1:4051084/1‑71 (MQ=255)                                                     cTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGc                  <  1:254615/71‑1 (MQ=255)                                                       aCCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                      <  1:5230191/65‑1 (MQ=255)                                                       aCCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                      <  1:1228423/65‑1 (MQ=255)                                                        ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCt                                                   <  1:991978/35‑1 (MQ=255)                                                        ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                    <  1:5564794/50‑1 (MQ=255)                                                        ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAAc                                    <  1:2626268/50‑1 (MQ=255)                                                        ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt               >  1:1869642/1‑71 (MQ=255)                                                        ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt               >  1:1876190/1‑71 (MQ=255)                                                        ccTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCatt               <  1:2516331/71‑3 (MQ=255)                                                          ttcttTGAAATTACGCTTCAGGATAATCCCGCTaa                                                 >  1:3772504/1‑35 (MQ=255)                                                          ttcttTGAAATTACGCTTCAGGATAATCCCGCTAACTTg                                             <  1:5387738/39‑1 (MQ=255)                                                           tcttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACt                                   >  1:1237771/1‑48 (MQ=255)                                                            cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                  >  1:3856700/1‑48 (MQ=255)                                                            cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACtt                                  >  1:4605907/1‑48 (MQ=255)                                                            cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTg           >  1:3984667/1‑71 (MQ=255)                                                             tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc          >  1:2760457/1‑71 (MQ=255)                                                             tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc          >  1:513687/1‑71 (MQ=255)                                                              ttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATcc                         <  1:2144178/55‑1 (MQ=255)                                                              gtGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCTGCTATACGATAAGCATTGGTGCg         <  1:2343094/70‑1 (MQ=255)                                                               tGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCgg        >  1:1375392/1‑70 (MQ=255)                                                                gAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGa        <  1:5551224/70‑1 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGt            >  1:1981234/1‑65 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg           >  1:1572776/1‑66 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg           >  1:3001702/1‑66 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg           >  1:536100/1‑66 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCggc       >  1:4205515/1‑68 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc      <  1:4455268/71‑1 (MQ=255)                                                                 aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc      >  1:2215143/1‑71 (MQ=255)                                                                 aaaTTACTCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGAc       <  1:5693475/70‑1 (MQ=255)                                                                 aaaTTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGAcc      >  1:2856905/1‑71 (MQ=255)                                                                   aTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg              >  1:5356662/1‑61 (MQ=255)                                                                   aTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg              >  1:5135467/1‑61 (MQ=255)                                                                   aTTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCg    <  1:2536126/71‑1 (MQ=255)                                                                    ttACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCgg   <  1:4605918/71‑1 (MQ=255)                                                                    ttACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCgg   <  1:2171877/71‑1 (MQ=255)                                                                      aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                        >  1:4701667/1‑48 (MQ=255)                                                                      aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                        >  1:4708280/1‑48 (MQ=255)                                                                      aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCg                        >  1:1544449/1‑48 (MQ=255)                                                                      aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGa  <  1:2308429/70‑1 (MQ=255)                                                                      |                                                                      CTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGA  >  Exported/313783‑313921

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