Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,602,723	C→T	5.5%	R1147R (CGC→CGT)	bpr →	bacillopeptidase F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,602,723	0	C	T	5.5%	248.1 / 2.8	73	R1147R (CGC→CGT)	bpr	bacillopeptidase F
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base T (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

AGTCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGGAGAGAAAG  >  BS168_reference_genome/1602593‑1602859                                                                                                                                   |                                                                                                                                         aGTCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGt                                                                                                                                  >  7:273072/1‑139 (MQ=255)   tCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGTAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGtt                                                                                                                                 >  8:320485/1‑138 (MQ=255)   tCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGtt                                                                                                                                 >  3:254536/1‑138 (MQ=255)   tCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGtt                                                                                                                                 >  2:349313/1‑138 (MQ=255)   tCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCAGGCCGCAACTGGCAGTATATAGATTGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATTTTATTCGCTCATTAGtt                                                                                                                                 >  4:126408/1‑138 (MQ=255)       aaCCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTa                                                                                                                            <  1:215265/139‑1 (MQ=255)             aTCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGTGCATTAGTTGATTATGAAgc                                                                                                                      >  5:314986/1‑139 (MQ=255)             aTCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGTGCATTAGTTGATTATGAAgc                                                                                                                      >  5:309703/1‑139 (MQ=255)             aTCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAgc                                                                                                                      >  2:62882/1‑139 (MQ=255)                cGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAgctg                                                                                                                    <  8:273072/138‑1 (MQ=255)                 gTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAgctgct                                                                                                                  <  1:349313/139‑1 (MQ=255)                  tATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGctgctg                                                                                                                 >  7:129301/1‑139 (MQ=255)                       aacTGGATTGGCAATTGTAGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTAcc                                                                                                             <  3:350089/136‑1 (MQ=255)                       aCCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTAcc                                                                                                             <  3:282938/138‑1 (MQ=255)                           ggATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCt                                                                                                            >  7:24468/1‑135 (MQ=255)                                       gATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCga                                                                                            <  1:87495/139‑1 (MQ=255)                                            aaCCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       <  1:165403/139‑1 (MQ=255)                                            aaCCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       <  3:17603/139‑1 (MQ=255)                                            aaCCGGCCAGAAATCCGGCCGCAACTGGCAGTATATAGATCGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCAATAGTTGATTATGAAGCTGCTGTACCAGAGCTTACTTCACCgacagac                                                                                       <  3:126408/139‑1 (MQ=255)                                               cGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACaaa                                                                                    >  2:324369/1‑139 (MQ=255)                                                 gCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       >  3:63304/1‑134 (MQ=255)                                                 gCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       <  4:63304/134‑1 (MQ=255)                                                 gCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAAt                                                                                   >  8:250339/1‑138 (MQ=255)                                                 gCCAGAATTCCGGCCGAAACTGGCAGTATATAGATGGAAAGTGGCAGCGAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAAt                                                                                   >  3:345336/1‑138 (MQ=255)                                                    aGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGTGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               <  6:309703/139‑1 (MQ=255)                                                    aGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGTGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               <  6:314986/139‑1 (MQ=255)                                                    aGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               <  3:266805/139‑1 (MQ=255)                                                    aGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               <  8:16259/139‑1 (MQ=255)                                                      aaTTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  7:126366/139‑1 (MQ=255)                                                      aaTTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  1:324369/139‑1 (MQ=255)                                                      aaTTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  5:326830/139‑1 (MQ=255)                                                      aaTTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  7:43937/139‑1 (MQ=255)                                                      aaTTCCGCCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  1:62882/139‑1 (MQ=255)                                                       aTTCCGGCCTCAACTGGAAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  4:3976/138‑1 (MQ=255)                                                           cGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACaaa                                                                         <  5:68906/138‑1 (MQ=255)                                                           cGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAAt                                                                        <  2:119095/139‑1 (MQ=255)                                                           cGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAAt                                                                        <  1:318359/139‑1 (MQ=255)                                                              ccGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAg                                                                     <  4:62346/139‑1 (MQ=255)                                                                       gCAGTATATAGATGGAAAGTGGTAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGt                                                             <  7:320485/138‑1 (MQ=255)                                                                       gCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGt                                                             <  1:999/138‑1 (MQ=255)                                                                       gCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGt                                                             <  7:250339/138‑1 (MQ=255)                                                                            atataGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACATACAAATCATACACAAATAAGGATAGCGTCACTg                                                        >  5:356127/1‑138 (MQ=255)                                                                            atataGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               >  8:40742/1‑115 (MQ=255)                                                                            atataGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               <  7:40742/115‑1 (MQ=255)                                                                                 gATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAg                                                  >  4:35226/1‑139 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGtt                                                                                                                                 <  3:97423/56‑1 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGtt                                                                                                                                 >  4:97423/1‑56 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       <  5:303866/98‑1 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       >  6:303866/1‑98 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGaaa                                              >  1:345528/1‑139 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGaaa                                              >  4:325282/1‑139 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGaaa                                              >  4:259226/1‑139 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGAAAAAGGaaa                                              >  6:367862/1‑139 (MQ=255)                                                                                     gAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGAGTCACTGTAAAAGGaaa                                              >  5:301793/1‑139 (MQ=255)                                                                                       aaGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       >  3:218298/1‑96 (MQ=255)                                                                                       aaGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCgacagac                                                                                       <  4:218298/96‑1 (MQ=255)                                                                                       aaGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATCCACAAATAAGGATAGCGTCACTGTAAAAGGAAACg                                            >  8:354432/1‑139 (MQ=255)                                                                                        aggcgCACCCAGGTGACAATGCGTAGGGCAACTATCTGATACGCGCATTCGTGGATTATGACGCTGCTGTCCCTTACTTAACTTCACCGACAGACAAAACATACACAAATACCGATAGCGTCACTGAAAAAGGAAACGc                                           <  2:331358/135‑1 (MQ=255)                                                                                                   ggTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGca                                <  3:325282/139‑1 (MQ=255)                                                                                                    gTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGcac                               >  4:245524/1‑139 (MQ=255)                                                                                                           agggtaTGGAAAATATATTATTCGCTCATTAGTTGATTATGAAGCTGCTTTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGCAAACGCTTCTCCTGGCACAACGGTa                        <  4:182715/134‑1 (MQ=255)                                                                                                           aGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTa                        <  3:245524/139‑1 (MQ=255)                                                                                                           aGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTa                        >  3:69345/1‑139 (MQ=255)                                                                                                           aGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACACACAAATCATAAACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTa                        <  5:367862/139‑1 (MQ=255)                                                                                                             cGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACATACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTaca                      <  6:356127/139‑1 (MQ=255)                                                                                                               gATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             <  6:353195/82‑1 (MQ=255)                                                                                                               gATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATaca                                                                             >  5:353195/1‑82 (MQ=255)                                                                                                                  ggCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACAtt                  >  3:182715/1‑138 (MQ=255)                                                                                                                          tatGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACg                                            >  4:320094/1‑104 (MQ=255)                                                                                                                          tatGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACg                                            <  3:320094/104‑1 (MQ=255)                                                                                                                           atGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTAATTCACCTACAGAAAAATCATACACAAATAAGGATAGCGCCACTGTAAACGTAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGga         <  8:24468/138‑1 (MQ=255)                                                                                                                               ttCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGgagaga     >  3:302351/1‑138 (MQ=255)                                                                                                                               ttCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGgagaga     >  1:4445/1‑138 (MQ=255)                                                                                                                                 cgcgCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               >  7:206808/1‑62 (MQ=255)                                                                                                                                 cgcgCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATa                                                                               <  8:206808/62‑1 (MQ=255)                                                                                                                                 cgcgCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGt                                                             <  6:301177/80‑1 (MQ=255)                                                                                                                                 cgcgCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGt                                                             >  5:301177/1‑80 (MQ=255)                                                                                                                                  gcgcATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGTATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGGAGagaaag  >  3:220295/1‑138 (MQ=255)                                                                                                                                  gcgcATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGGAGagaaag  >  3:369714/1‑138 (MQ=255)                                                                                                                                  gcgcATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAAGAAGGATAGCTTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAAAGCAGagaaag  >  5:4408/1‑138 (MQ=255)                                                                                                                                   |                                                                                                                                         AGTCGAAACCTGATCCGTATTCACCTGGATTGGCAATGGATGAAACCGGCCAGAATTCCGGCCGCAACTGGCAGTATATAGATGGAAAGTGGCAGCCAGGTGACAAAGCGGATGGCAACTATATGATTCGCGCATTAGTTGATTATGAAGCTGCTGTACCTGAGATTACTTCACCGACAGACAAATCATACACAAATAAGGATAGCGTCACTGTAAAAGGAAACGCTTCTCCTGGCACAACGGTACACATTTATAATGGAGAGAAAG  >  BS168_reference_genome/1602593‑1602859

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