Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,220,375	A→C	16.7%	N176T (AAC→ACC)	oppA →	oligopeptide ABC transporter (binding lipoprotein)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,220,375	0	A	C	16.7%	63.2 / 2.5	39	N176T (AAC→ACC)	oppA	oligopeptide ABC transporter (binding lipoprotein)
Reads supporting (aligned to +/- strand):  ref base A (22/10);  new base C (2/5);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 8.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTAATAATGAATCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAGATGATTATGTATCAAACGGGCCGTTCAAAATGACGG  >  BS168_reference_genome/1220239‑1220512                                                                                                                                         |                                                                                                                                          ccTAAAAATGAAACAAAATACGCTTACCTGCTCAACAACATAAAAGGGGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTTAATGACAAACCGCTGAAGGTTCAATTAAAAAAcc                                                                                                                                         <  8:46798/139‑1 (MQ=255)    aataatTAAACACAATCCGCTCACCACCTAAAATACAAAAAAGGTGCTGAAGCGGCCAAACCCGCAAACGCCAGCCTACACTCTGGGGCAGTAAAAGCTGTCAATCACAAAACGCTGAAGGTTCAATAAAAAACCCCGa                                                                                                                                      <  8:145943/139‑1 (MQ=255)           aaTCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAATGCTGTGAATGACAAAACGCGGAAGGTTGAATTAACTAACCCGACTCCata                                                                                                                               >  3:204376/1‑139 (MQ=255)           aaTCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGCCTCCata                                                                                                                               >  4:26531/1‑139 (MQ=255)           aaTCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCata                                                                                                                               >  6:134692/1‑139 (MQ=255)           aaTCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCata                                                                                                                               >  1:32441/1‑139 (MQ=255)             tCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATAtt                                                                                                                             >  5:196166/1‑139 (MQ=255)             acacaATACGCTTACCAGCAATACTACATAACCGGTGCTGAAGCGGCTAAACCCGTAAACGCCCGCCTACCCGAGGGGGCATAAAATGCTTTGAATCAAAAAACGCTGATGGTTTAAAAAAAAACCCCGACTCCATAtt                                                                                                                             <  4:120717/138‑1 (MQ=255)               acaATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCTACTCCATAttt                                                                                                                            >  2:88337/1‑138 (MQ=255)                          aCCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAAc                                                                                                                >  3:235094/1‑139 (MQ=255)                            cAGCTCTACTAAAAAACAGGTGCTGACGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCATTAAAAGCTGTGAATGACAAACCGCTGAAGGTTTAATTAAATAACCCGACTCCATATTTCACTGAATTAACTg                                                                                                              <  3:242655/139‑1 (MQ=255)                              gCTCTACAACAAAAACGGTGCTGACGCGGCTAATACCGCAAAAGGCAGCCTAGCCGATGGGCCAGTAAAAGGTGTCAATCAAAAACCGCTGACGGTTGAATAAAAAAACCCGACTCCATATTTCACTGAATTAACTGCg                                                                                                            <  8:55022/139‑1 (MQ=255)                               ctctACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGt                                                                                                           >  5:201507/1‑139 (MQ=255)                                   acaaaaTAAATGGTGCTGAAGGGGCTAATACCGGAAAAGGCAGCCTAGACGATGGGGCAGTAAAAGCTGTGAATGAAAAAACGCTGAAGGTTTAATTAAATACCCCGACTCCATATTTCACTGAATTAACTGCGTTCta                                                                                                       <  1:92699/134‑1 (MQ=255)                                   aacacAAAAACGGTGCGGAAGCGGCTAAACCCGAAACCGGCCGCCTCGCCGTGGTGGCATTAACAGCGGTTAATCAAAAACCGCTGACGGTTTAATAAAAAACCCCGACTCCATATTTCACTGAATTAACTGCGTTCta                                                                                                       <  6:196166/136‑1 (MQ=255)                                      acATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATAc                                                                                                    >  7:80057/1‑139 (MQ=255)                                      acATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATAc                                                                                                    >  4:49764/1‑139 (MQ=255)                                         tAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGta                                                                                                 >  7:184367/1‑139 (MQ=255)                                            aaGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATACCCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTatat                                                                                              >  4:152600/1‑139 (MQ=255)                                                 gCTGAAGCGGCGAGTACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGa                                                                                         >  1:210031/1‑139 (MQ=255)                                                 gCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGtatactgtct                                                                                         >  3:226942/1‑133 (MQ=255)                                                 gCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGtatactgtct                                                                                         >  3:176979/1‑133 (MQ=255)                                                          gcgAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAg                                                                                 >  1:191385/1‑138 (MQ=255)                                                             aaTACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATACCCCGCCTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGaaaa                                                                             >  8:87756/1‑139 (MQ=255)                                                                            aGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAAt                                                               >  1:225207/1‑138 (MQ=255)                                                                                    ctttgtgGCATTAAAAGGTGTGAAACAAAAAACGCTGAAGGTTTAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTgg                                                      <  5:140198/136‑1 (MQ=255)                                                                                     gATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGa                                                     <  8:20100/139‑1 (MQ=255)                                                                                        gtgGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAAt                                                   >  8:257243/1‑138 (MQ=255)                                                                                        gtgGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATa                                                  >  4:27649/1‑139 (MQ=255)                                                                                        gtgGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGAATATGCAGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGCAATa                                                  >  8:139706/1‑139 (MQ=255)                                                                                               tAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAACTGc                                           >  6:82804/1‑139 (MQ=255)                                                                                                 aaaGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCg                                         <  3:49764/139‑1 (MQ=255)                                                                                                  aaGCTGTCAATCACAAAACGCTGACGGGTCAATAAAAAACCCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCgg                                        <  5:54987/139‑1 (MQ=255)                                                                                                    gatttcaaTCACAAAACGCTGATGTTTTAATTAAAAAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGgag                                      <  4:156444/133‑1 (MQ=255)                                                                                                        tgAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAgatga                                  >  4:139105/1‑139 (MQ=255)                                                                                                           aTGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAGATGAtt                                <  3:139105/138‑1 (MQ=255)                                                                                                           aTGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAGATGAtt                                <  5:134692/138‑1 (MQ=255)                                                                                                                         aaGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAAATGCAGAGAAAAATAAAAAGTGGAATACAAATGACGGAGATGATTATGTATCAAACGGGc                 >  1:204382/1‑139 (MQ=255)                                                                                                                                      aTAACCCGACTCCATAATTCACTGAATTAACTGCTGTCTATACTTATATGCCGATCAATACGAAAATTTCAGTGAAAAATAAAACGGGCAATACAACTGCCGGAGATGATTATGTATCAAACGGGCCGTTCAAAATGAc    <  8:239317/139‑1 (MQ=255)                                                                                                                                        aaCCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAGATGATTATGTATCAAACGGGCCGTTCAAAATGACgg  <  8:184367/139‑1 (MQ=255)                                                                                                                                        aaCCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAGATGATTATGTATCAAACGGGCCGTTCAAAATGACgg  <  7:180466/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                          CCTAATAATGAATCACAATACGCTTACCAGCTCTACTACATAAAAGGTGCTGAAGCGGCGAATACCGGAAAAGGCAGCCTAGACGATGTGGCAGTAAAAGCTGTGAATGACAAAACGCTGAAGGTTGAATTAAATAACCCGACTCCATATTTCACTGAATTAACTGCGTTCTATACGTATATGCCGATCAATAAGAAAATTGCAGAGAAAAATAAAAAGTGGAATACAAATGCCGGAGATGATTATGTATCAAACGGGCCGTTCAAAATGACGG  >  BS168_reference_genome/1220239‑1220512

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