Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,838,721	T→G	Y587S (TAC→TCC)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,838,721	0	T	G	97.0%	122.5 / ‑3.9	33	Y587S (TAC→TCC)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (19/13);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTG  >  REL606/3838649‑3838790                                                                         |                                                                      tagataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACt                                                                    >  2:58767/1‑76 (MQ=255)  agataAGCTCTTTCNGCATTGCTTTCGGATAACGNCGCTCCATNCGATTGATCATTTTCTGCGTCGCGTTGGACTc                                                                   <  2:1175191/76‑1 (MQ=255)  agataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTc                                                                   <  1:844775/76‑1 (MQ=255)    ataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAg                                                                 <  2:379891/76‑1 (MQ=255)    ataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAg                                                                 >  2:305285/1‑76 (MQ=255)      aaGCCCTTTCAGCATTGCTTTCGGATAACGACGCCCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGAt                                                               >  2:1395339/1‑76 (MQ=255)         ctctTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATAcc                                                            <  1:330485/76‑1 (MQ=255)         ctctTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATAcc                                                            >  2:1245445/1‑76 (MQ=255)          tctTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCa                                                           >  2:258583/1‑76 (MQ=255)           ctTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAg                                                          <  1:354394/76‑1 (MQ=255)              tCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGttt                                                       >  1:278033/1‑76 (MQ=255)                   aTTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTa                                                  <  1:626845/76‑1 (MQ=255)                   aTTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTa                                                  >  1:661162/1‑76 (MQ=255)                       cTTTCGTCTAACGACTCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGCTTCTCTAACGc                                              >  2:956985/1‑76 (MQ=255)                        tttCGGATAAGGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCt                                             >  1:1319712/1‑76 (MQ=255)                        tttCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCt                                             >  2:843013/1‑76 (MQ=255)                             gATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCgcc                                        <  2:59324/76‑1 (MQ=255)                               tAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCgccag                                      <  2:885496/76‑1 (MQ=255)                               tAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCgccag                                      <  1:499257/76‑1 (MQ=255)                               tAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCgccag                                      >  1:295681/1‑76 (MQ=255)                               tAACGACGCTCCATACGATTGATCATTGTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCgccag                                      <  2:725742/76‑1 (MQ=255)                                aaCGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCgccagc                                     <  1:634618/76‑1 (MQ=255)                                 acgacgCTNCTTACGATTGNTCATTTTCTGCGTCGCGTTGGANTCAGATACCAGTTTCTCTAACGCTTCgccagcc                                    <  2:1188431/76‑1 (MQ=255)                                   gacgCTCCATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCaa                                  <  1:342897/76‑1 (MQ=255)                                         ccATACGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCgg                            >  1:425691/1‑76 (MQ=255)                                            tACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGc                         >  1:1364143/1‑76 (MQ=255)                                             aCGATTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCa                        >  1:1070313/1‑76 (MQ=255)                                                aTTGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTg                     <  1:899157/76‑1 (MQ=255)                                                 ttGATCATTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTgg                    >  2:542163/1‑76 (MQ=255)                                                       aTTTTCTGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTgg              >  1:574252/1‑76 (MQ=255)                                                         tttCTGCGTCGCGTTGGACTCAGATACCATTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGgtg            <  2:367736/76‑1 (MQ=255)                                                          ttCTGCGTCGCGTTGGACTCAGACCCCATTTTCTCAACCGCTTCCCCAGCCAATGCCGGTGCACTGGCGTTGgtgt           <  1:710147/76‑1 (MQ=255)                                                             tGCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCa        >  1:117019/1‑76 (MQ=255)                                                              gCGTCGCGTTGGACTTAGATACCAGTTTCTCTANCGCTTCGCCAGCCAATGCCGGTNCACTGGCGATNGTGTGCAg       >  2:1188978/1‑76 (MQ=255)                                                              gCGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAg       >  2:141073/1‑76 (MQ=255)                                                               cGTCGCGTTGGACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGc      >  2:56687/1‑76 (MQ=255)                                                                  cgctTTGGACTCGGNTCCCATTTTCCCTGACGCTNCGCGATCCNATGCCGGTGCACTGGCGTTAGTGTGCAGCGtt   <  2:1169707/76‑1 (MQ=255)                                                                   gcgTTGGCCTCAGATCCCAGTTTCTCAACCGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTg  <  2:569813/76‑1 (MQ=255)                                                                         |                                                                      TAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTG  >  REL606/3838649‑3838790

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑