Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,556,359	G→T	5.6%	E725D (GAG→GAT)	pycA →	pyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,556,359	0	G	T	5.6%	241.7 / 4.3	71	E725D (GAG→GAT)	pycA	pyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base T (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GATATCCTTGACAAGAACCGGACGAAGTACGACCTTGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGACGTGGCGGTCAG  >  BS168_reference_genome/1556222‑1556490                                                                                                                                          |                                                                                                                                    gaTATCCTTGACAAGAACCGGACGAAGTACGACCTTGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGc                                                                                                                                    >  2:33858/1‑139 (MQ=255)    atCCTTGACAAGAACCGGACGAAGTACGACCTTGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCg                                                                                                                                 >  4:246346/1‑139 (MQ=255)         tGACAAGAACCGGACGAAGTACGACCTTCCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCAGATTCACGCGATAAAAGATATTGAAGCGCTGTTAAAACCGCAGGCTGCATATCAGCTCGTTGCt                                                                                                                            >  5:25527/1‑139 (MQ=255)                          aGTACGACCTTGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGAt                                                                                                           <  3:246346/139‑1 (MQ=255)                               gACCTTGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACa                                                                                                      >  2:271519/1‑139 (MQ=255)                                    tGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCg                                                                                                 >  1:23210/1‑139 (MQ=255)                                    tGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCg                                                                                                 >  4:267401/1‑139 (MQ=255)                                    tGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCg                                                                                                 >  6:123863/1‑139 (MQ=255)                                    tGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCg                                                                                                 >  7:268599/1‑139 (MQ=255)                                    tGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCg                                                                                                 >  8:282851/1‑139 (MQ=255)                                    tGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGCCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCg                                                                                                 >  5:11321/1‑139 (MQ=255)                                      cATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGt                                                                                               <  2:23210/139‑1 (MQ=255)                                       atatTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGtt                                                                                              <  5:123863/139‑1 (MQ=255)                                       atatTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGtt                                                                                              <  3:333234/139‑1 (MQ=255)                                       atatTATACAACGATTGCGAAGTAGCTAGAGGCGGCCGGAGCCCATATTCTCTGGATTAAAGATAAGGCAGGGCTGTTAAAACCGCAGGCTTCATATGAGCACGTTTCTTCCTTGAAAGAAACGATCGACATTCCGGtt                                                                                              <  3:248182/139‑1 (MQ=255)                                           tataCATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCAcc                                                                                          >  5:330672/1‑139 (MQ=255)                                           tataCATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATACAGGTTCAcc                                                                                          >  4:357421/1‑139 (MQ=255)                                            ataCATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCAcc                                                                                          <  2:53937/138‑1 (MQ=255)                                               cATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGATTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCa                                                                                      <  3:34582/139‑1 (MQ=255)                                                aTCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCa                                                                                      >  8:54950/1‑138 (MQ=255)                                                aTCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCa                                                                                      >  7:2244/1‑138 (MQ=255)                                                aTCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTGCACCTTCa                                                                                      >  5:240238/1‑138 (MQ=255)                                                 tCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATa                                                                                    <  8:268599/139‑1 (MQ=255)                                                  cGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATAc                                                                                   >  5:387370/1‑139 (MQ=255)                                                   gATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCTGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  <  8:380721/139‑1 (MQ=255)                                                   gATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCAc                                                                                           >  8:295144/1‑130 (MQ=255)                                                   gATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCAc                                                                                           <  7:295144/130‑1 (MQ=255)                                                    aTGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTTAAACCGCAGGCAGCATATGAGCTCGTTTCTGCGTTGATAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  >  2:336777/1‑138 (MQ=255)                                                    aTGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAACCCGCAGGCTGCCTATGAGCTCGTTGCTGCGTTGATAGAAACGATCGACATTCCGGTTCCCCTTCAGACGc                                                                                 >  8:384092/1‑139 (MQ=255)                                                    aTGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGATCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  >  2:9315/1‑138 (MQ=255)                                                    aTGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGATCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  >  2:226595/1‑138 (MQ=255)                                                    aTGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  >  8:309506/1‑138 (MQ=255)                                                          aaGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATa                                                                           <  3:411889/139‑1 (MQ=255)                                                            ggAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCTTTGAAAGAAACGAt                                                                                                           <  8:297671/105‑1 (MQ=255)                                                            ggAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGAt                                                                                                           >  7:297671/1‑105 (MQ=255)                                                             gtgtttGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGTTATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTTAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACg                                                                         <  6:240238/134‑1 (MQ=255)                                                              aGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAg                                                                       <  5:391977/139‑1 (MQ=255)                                                                    cggcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGaaa                                                                 <  3:305992/138‑1 (MQ=255)                                                                      gcggcCTGATCCCATATTCTCGTGAGTAAAGATCTGGCACGGATGTTAAAACCGCAGGCTGCATAGTAGCGCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  4:229725/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATTAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  3:362148/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGATCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  1:9315/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGATCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  1:226595/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  1:271519/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  5:355577/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  8:248542/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  7:54950/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  6:330672/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  6:387370/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  3:267401/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATATTCTCGGGATAAAAGATATGGCAGGGCTGGTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  5:132645/139‑1 (MQ=255)                                                                      gcggcCGGAGCCCATAATCTCGGGATTTAAGATATGGCAGGGCTGTCAACACCGCAGGCGGCATATGAGCTCGGTTCTGCGGGGAACGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACg                                                               <  4:210836/139‑1 (MQ=255)                                                                        ggcCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGt                                                             <  1:33858/139‑1 (MQ=255)                                                                          ccGGAGCTCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  <  1:272802/116‑1 (MQ=255)                                                                          ccGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACg                                                                                  >  2:272802/1‑116 (MQ=255)                                                                          ccGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTAt                                                           >  2:64617/1‑139 (MQ=255)                                                                                    tatTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATGCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTAtatgtat                                                 >  6:361649/1‑139 (MQ=255)                                                                                           gggATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAg                                          >  2:277432/1‑139 (MQ=255)                                                                                                      tatGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGcc                               <  8:2244/139‑1 (MQ=255)                                                                                                       atGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGt                                                                                               >  5:155521/1‑74 (MQ=255)                                                                                                       atGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGt                                                                                               <  6:155521/74‑1 (MQ=255)                                                                                                          gCAGGGCTGTTAATACCGCAGGCTGCATATTAGCTCGTTTCTGCGTTGACAGAATCTATCGACATTCCGGTTCACGTTCATTCGCTTGATACGAGAGGAAAGGGTATTTATATGTATGCGAAAGCTGTTGTAGCCGGCg                           <  6:25527/139‑1 (MQ=255)                                                                                                             gggCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTg                        <  7:282851/139‑1 (MQ=255)                                                                                                              ggCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGa                       >  3:409757/1‑139 (MQ=255)                                                                                                              ggCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAAAGATCGACATTCCGGTTCAAATTCATACGCATTATACGAGCGGAAACGGTATTTATATTTATGCGAAAGCTGTTGAAGACGGAGTTGa                       >  7:248776/1‑139 (MQ=255)                                                                                                                 tGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGatat                    >  4:330991/1‑139 (MQ=255)                                                                                                                        aCCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGAc             >  3:403205/1‑139 (MQ=255)                                                                                                                           gCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGACGTg          <  1:277432/139‑1 (MQ=255)                                                                                                                             aGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAACCGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGACGTGGc        >  1:361927/1‑139 (MQ=255)                                                                                                                             aGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGACGTGGc        >  3:66782/1‑139 (MQ=255)                                                                                                                                cTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGACGTGGCGGt     <  5:69508/139‑1 (MQ=255)                                                                                                                                   cATATGAGCTCGTTTCTGCGTGGAAAGAAGCGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGATGCCGGCGTTGATATCATAGCCGTGGCGGTccg  >  2:252563/1‑137 (MQ=255)                                                                                                                                          |                                                                                                                                    GATATCCTTGACAAGAACCGGACGAAGTACGACCTTGCATATTATACATCGATGGCGAAGGAGCTTGAGGCGGCCGGAGCCCATATTCTCGGGATTAAAGATATGGCAGGGCTGTTAAAACCGCAGGCTGCATATGAGCTCGTTTCTGCGTTGAAAGAAACGATCGACATTCCGGTTCACCTTCATACGCATGATACGAGCGGAAACGGTATTTATATGTATGCGAAAGCTGTTGAAGCCGGCGTTGATATCATAGACGTGGCGGTCAG  >  BS168_reference_genome/1556222‑1556490

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