Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	313,903	G→A	22.4%	A62A (GCC→GCT)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	313,903	0	G	A	22.4%	124.2 / 31.5	67	A62A (GCC→GCT)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base A (12/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCC  >  Exported/313834‑313968                                                                      |                                                                  gCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAgg                                                                  >  1:4051084/1‑71 (MQ=255)  cTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGc                                                                 <  1:254615/71‑1 (MQ=255)     ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt                                                              >  1:1876190/1‑71 (MQ=255)     ccTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACt                                                              >  1:1869642/1‑71 (MQ=255)     ccTTCAGTGAAATTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCatt                                                              <  1:2516331/71‑3 (MQ=255)         cttTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTg                                                          >  1:3984667/1‑71 (MQ=255)          tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                                                         >  1:2760457/1‑71 (MQ=255)          tttGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                                                         >  1:513687/1‑71 (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:3001702/1‑66 (MQ=255)              aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                                                          >  1:1572776/1‑66 (MQ=255)              aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTg                                                          >  1:536100/1‑66 (MQ=255)              aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCggc                                                      >  1:4205515/1‑68 (MQ=255)              aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc                                                     >  1:2215143/1‑71 (MQ=255)              aaaTTGCGATTCAGTACGATCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTGGTGCGGcc                                                     <  1:4455268/71‑1 (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)                   aCGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGa                                                 <  1:2308429/70‑1 (MQ=255)                     gCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGc                                              >  1:139169/1‑71 (MQ=255)                      cTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCt                                             <  1:1138572/71‑1 (MQ=255)                        tCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTc                                            >  1:4285273/1‑70 (MQ=255)                         cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCgg                                          >  1:3156673/1‑71 (MQ=255)                         cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCgg                                          >  1:723581/1‑71 (MQ=255)                         cAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCg                                           <  1:5147570/70‑1 (MQ=255)                          aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                                                             >  1:887874/1‑51 (MQ=255)                          aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                                                             >  1:1081018/1‑51 (MQ=255)                          aGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTg                                                             >  1:2580969/1‑51 (MQ=255)                           ggATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGtt                                        >  1:1133632/1‑71 (MQ=255)                            gATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtg                                       <  1:663073/71‑1 (MQ=255)                              tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                                                                >  1:5552838/1‑44 (MQ=255)                              tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                                                                >  1:480781/1‑44 (MQ=255)                              tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:2803890/71‑1 (MQ=255)                              tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:157456/71‑1 (MQ=255)                              tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:4601654/71‑1 (MQ=255)                              tAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:4781529/71‑1 (MQ=255)                                aTCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGa                                    <  1:28518/70‑1 (MQ=255)                                aTCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                                             >  1:1602933/1‑45 (MQ=255)                                aTCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                                             >  1:3751927/1‑45 (MQ=255)                                aTCCCGCTAACCTGTGGGCGAACTTCGGCTATACGATAAGCATTg                                                             >  1:180262/1‑45 (MQ=255)                                   ccGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTg                                >  1:4630114/1‑71 (MQ=255)                                     gCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCa                              <  1:5687358/71‑1 (MQ=255)                                      cTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag                             <  1:2953582/71‑1 (MQ=255)                                      cTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag                             <  1:4744225/71‑1 (MQ=255)                                      cTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag                             <  1:2961360/71‑1 (MQ=255)                                       tAACTTGAGGAGGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:4120954/62‑1 (MQ=255)                                       tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCa                                                                >  1:2529336/1‑35 (MQ=255)                                       tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGtt                                        >  1:464566/1‑59 (MQ=255)                                       tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:1502636/62‑1 (MQ=255)                                       tAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:2541091/62‑1 (MQ=255)                                        aaCTTGAGGACGAACTTCTGCGGTCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCagag                           >  1:5040791/1‑71 (MQ=255)                                        aaCTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:2107542/61‑1 (MQ=255)                                          cTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGAt                                   >  1:4652778/1‑61 (MQ=255)                                           ttGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTtgtg                                     <  1:1048885/58‑1 (MQ=255)                                             gAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGAt                                   >  1:4345115/1‑58 (MQ=255)                                                gACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCag                             <  1:2579096/61‑1 (MQ=255)                                                         cTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGt           <  1:3090748/70‑1 (MQ=255)                                                           gCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGtttt                 <  1:1925149/62‑1 (MQ=255)                                                            cGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTtac        <  1:4439050/70‑1 (MQ=255)                                                                 cGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGcc  <  1:4939065/71‑1 (MQ=255)                                                                 cGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGcc  <  1:3288105/71‑1 (MQ=255)                                                                      ggCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGtt                   <  1:108045/49‑1 (MQ=255)                                                                      |                                                                  GCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTACTACGCC  >  Exported/313834‑313968

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