Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	409,328	A→G	7.7%	D208D (GAT→GAC)	tcyC ←	cystine ABC transporter (ATP‑binding protein)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	409,328	0	A	G	7.7%	152.3 / 4.3	52	D208D (GAT→GAC)	tcyC	cystine ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base A (17/31);  new base G (2/2);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CTCTGTTTTTCTTATTACAGCGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTGAGGTCGGTTCGTCAAA  >  BS168_reference_genome/409194‑409462                                                                                                                                       |                                                                                                                                       ctctGTTTTTCTTATTACAGCGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCGTCCg                                                                                                                                    >  3:47162/1‑139 (MQ=255) ctctGTTTTTCTTATTACAGCGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCGTCCg                                                                                                                                    >  3:41545/1‑139 (MQ=255)      tttttCTTATTACAGCGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCa                                                                                                                                  <  4:184438/136‑1 (MQ=255)              attaCAGCGGGTTCAAAATCCGGTTTAAGCACCGCTGTGTCCGATCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCGGCTCCCCGAGAACGCCGCAGTCGAGGAAGATGACTTCATCCGCACCCTCCTgcgc                                                                                                                       >  6:264839/1‑139 (MQ=255)              attaCAGCGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATACCGCCGCCGTCGATGAAGCTGACTTCATCCGCAACCTCCTgcgc                                                                                                                       >  8:249468/1‑139 (MQ=255)                     cGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACtt                                                                                                                 >  7:73844/1‑138 (MQ=255)                     cGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTg                                                                                                                >  8:339204/1‑139 (MQ=255)                        gTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGCACTTGAtt                                                                                                             >  6:86396/1‑139 (MQ=255)                              actCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAACATTGGCTCCGGCGGGCCCTGCTCCACGAAAACGCCGCCGTCGATGAAGCTGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTgg                                                                                                       <  3:109523/137‑1 (MQ=255)                                   ggTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCac                                                                                                  >  2:94646/1‑139 (MQ=255)                                    gTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacg                                                                                                 <  4:334198/139‑1 (MQ=255)                                     tttAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacga                                                                                                <  6:124631/139‑1 (MQ=255)                                     tttAAGAACCGCTGTGTCCGTTCTTCTTTTGGGGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacga                                                                                                <  4:52827/139‑1 (MQ=255)                                       tAAGAACCGCTGTGTCGGTTCTTCTTTTGGTGCTGAGAAAATTTGCTGCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacgac                                                                                               <  3:311536/138‑1 (MQ=255)                                       tAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacgac                                                                                               <  5:173312/138‑1 (MQ=255)                                           accggCTGTGTCCGTTCTTCTTTTGGTGCGTAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGCCTTCATCCGCAACCTCCTGCGCTAACTTGATTTCGTGGGGCACGACGACCa                                                                                          <  6:43117/135‑1 (MQ=255)                                           accctgtgttTCGGTTTTTTTTTTGGGGCGTTGAAATTTTGGCCGGGCGGGCCCTGCCCCACTTACCCGCCGCGGTGTATTACGATCACTTCATCCGCAACCTCCTGCGCGAACTTTATTTCGTGGGTCCCGACGACCa                                                                                          <  6:270409/133‑1 (MQ=255)                                           acccGCTGTGTCGTTTCTTCTTTTGGTGGGGAGAAAATTGGCGCCGGCGGGCCCTGAACCACGAAAACGCCGCCGTCGATTAAGTTGACTTCATCCGCAACCTCCTGCTCGAACTTGATTTCGTGGGTCACGACGACCa                                                                                          <  8:365747/137‑1 (MQ=255)                                              cGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacgac                                                                                               >  2:54892/1‑131 (MQ=255)                                              cGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacgac                                                                                               <  1:54892/131‑1 (MQ=255)                                                   tgtCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACATCATTTCGTGGGTCACGACGACAATGGTCCAg                                                                                  >  8:18474/1‑139 (MQ=255)                                                     tCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTAATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCacc                                                                                 >  5:242254/1‑136 (MQ=255)                                                         ttttttttttGGTGCGGATAAATGTTGCTCGGGGGGCCCGTGCTCCACGATAACGCGGCCGTCGATGAAGATGAATTCATCCGCAACCTCCTGCGCTAACTTGTTTTCGTGGGTCACGACGACCATGGTCCAGCCTTCa                                                                            <  6:242254/133‑1 (MQ=255)                                                         ttcttcTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCa                                                                            <  7:10139/139‑1 (MQ=255)                                                           cttcttTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCGTCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCAt                                                                           <  4:41545/138‑1 (MQ=255)                                                           cttcttTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCGTCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCAt                                                                           <  4:47162/138‑1 (MQ=255)                                                               ttttGGGGGGGACAAAATTTGCGCCGGCGGGCCCTGCTCCACGATAACGCCGCCGTCGATGAAGAGGACTTCATCCGCAACCTCCTGCGCGAACTGGATTTCGTGGGGCACGACGACCATGGTCCAGCCTTCATTGGcc                                                                      <  4:428166/139‑1 (MQ=255)                                                                 ttGGGGCGGAGAAAATTTGCTCCGGCGGACCCTGCCCCACGATAACGCCGCCGGCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGGGGGTCACGACGACCATGGTCCAGCCTTCATTGGcc                                                                      <  5:399260/137‑1 (MQ=255)                                                                         gagaAAATTCGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCTAACTTGATTTCGTGGGGCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGa                                                            <  8:225886/139‑1 (MQ=255)                                                                          agaAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGAt                                                           <  8:73844/139‑1 (MQ=255)                                                                           gaAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATa                                                          <  1:175973/139‑1 (MQ=255)                                                                           gaAAATGTGCTCCGGCGGTCCCTGCTCCACTAAAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCACGTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATa                                                          <  1:88354/139‑1 (MQ=255)                                                                                     gccGGCGGTCCCGGCTCCACTATAAGGCCGCCGTCGATGAATATGACTTCATCCGCAACCTCCTGCGCTAACTTTATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCa                                                <  2:194761/138‑1 (MQ=255)                                                                                       cggcggGCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCa                                                <  7:96179/137‑1 (MQ=255)                                                                                               ccTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGAc                                      <  3:315493/139‑1 (MQ=255)                                                                                               ccTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGAAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGAc                                      <  7:60944/139‑1 (MQ=255)                                                                                               ccTGCGCCACGATAAGGCCGCCGTCGATGAAGATGTCTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGGCACGACGAACATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGAc                                      <  7:18474/139‑1 (MQ=255)                                                                                                      cACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCt                                 <  5:219123/137‑1 (MQ=255)                                                                                                        cGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGTCCAAGTCCTTGATAAACTTCAGCACCTCTCCGACAAGCTCggg                             >  3:184998/1‑139 (MQ=255)                                                                                                             aCGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCaa                        >  2:113947/1‑139 (MQ=255)                                                                                                                ccgccgTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacgac                                                                                               <  4:182584/65‑1 (MQ=255)                                                                                                                ccgccgTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCacgac                                                                                               >  3:182584/1‑65 (MQ=255)                                                                                                                   ccgTCGATGAAGATGTCTTCATACGCAACCTCCTGCGCGAACTTGATTTAGTGGGTCACGAAGACCATGGTCCAGCCTTCATTGGGCAATTCCTTGATAACCTTAATCACCTCTACGACACGCAAGGGATCAAgcg                     >  6:56974/1‑136 (MQ=255)                                                                                                                   ccgTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGCCCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTg                  >  4:402881/1‑139 (MQ=255)                                                                                                                   ccgTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGAcc                                                                                           <  7:209336/66‑1 (MQ=255)                                                                                                                   ccgTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGAcc                                                                                           >  8:209336/1‑66 (MQ=255)                                                                                                                   ccgTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTg                  >  7:320645/1‑139 (MQ=255)                                                                                                                            aaGAGGACTTCATCCGCACCCTCCTGCGCGAACTTGATTTCGGGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTGAGGTCGGt          <  4:300520/138‑1 (MQ=255)                                                                                                                                 gACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATc                          <  3:380338/117‑1 (MQ=255)                                                                                                                                 gACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATc                          >  4:380338/1‑117 (MQ=255)                                                                                                                                  aCTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTGAGGTCGGTTCGTCa    <  1:169574/138‑1 (MQ=255)                                                                                                                                   cTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTGAGGTCGGTTCGTCaaa  <  5:86396/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       CTCTGTTTTTCTTATTACAGCGGGTTCAAAATCCGGTTTAAGAACCGCTGTGTCCGTTCTTCTTTTGGTGCGGAGAAAATTTGCTCCGGCGGTCCCTGCTCCACGATAACGCCGCCGTCGATGAAGATGACTTCATCCGCAACCTCCTGCGCGAACTTGATTTCGTGGGTCACGACGACCATGGTCCAGCCTTCATTGGCCAAGTCCTTGATAACCTTCAGCACCTCTCCGACAAGCTCGGGATCAAGCGCTGAGGTCGGTTCGTCAAA  >  BS168_reference_genome/409194‑409462

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: ‑