Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,073,367:2	+A	6.9%	coding (351/612 nt)	scoC ←	transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,073,367	2	.	A	6.9%	204.9 / 11.7	58	coding (351/612 nt)	scoC	transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
Reads supporting (aligned to +/- strand):  ref base . (23/31);  new base A (2/2);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

AATTTCCATAAAATCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCGTTCTTCCAGCTTTTTTGAAAAGTTGAATGCGGTTGATACGTGCATG  >  BS168_reference_genome/1073229‑1073505                                                                                                                                             |                                                                                                                                           aaTTTCCATAAAATCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc                                                                                                                                              <  7:233669/139‑1 (MQ=255)  aTTTCCATAAAATCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           <  3:252398/139‑1 (MQ=255)  aTTTCCATAAAATCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           <  8:65488/139‑1 (MQ=255)            aaTCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATATATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGaa                                                                                                                                 <  4:253478/139‑1 (MQ=255)            aaTCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGaa                                                                                                                                 <  5:211486/139‑1 (MQ=255)                 tcGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCt                                                                                                                            >  8:160504/1‑139 (MQ=255)                 tcGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCt                                                                                                                            >  2:159536/1‑139 (MQ=255)                 tcGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCt                                                                                                                            >  4:148174/1‑139 (MQ=255)                      gTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑CAGCTCCAGAATACCTCGGTa                                                                                                                       <  1:299013/139‑1 (MQ=255)                       tAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGCAAAAGCTCCAGAATACCTCgg                                                                                                                         >  1:31381/1‑138 (MQ=255)                       tAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCt                                                                                                                            <  4:193348/133‑1 (MQ=255)                       tAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCt                                                                                                                            >  3:193348/1‑133 (MQ=255)                       tAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTAc                                                                                                                      <  2:241415/139‑1 (MQ=255)                            atGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGCCCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACATCGGTACCTTcc                                                                                                                 >  6:6047/1‑139 (MQ=255)                               gCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCt                                                                                                                >  3:245598/1‑137 (MQ=255)                               gCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCa                                                                                                              <  5:125536/139‑1 (MQ=255)                               gCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCa                                                                                                              <  7:63655/139‑1 (MQ=255)                               gCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCa                                                                                                              <  7:72438/139‑1 (MQ=255)                                  tATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           <  2:159909/107‑1 (MQ=255)                                  tATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           >  1:159909/1‑107 (MQ=255)                                  tATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGt                                                                                                            >  8:92418/1‑138 (MQ=255)                                    tCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTc                                                                                                                                      <  6:286731/110‑1 (MQ=255)                                    tCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTc                                                                                                                                      >  5:286731/1‑110 (MQ=255)                                    tCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGaa                                                                                                                                 <  6:240999/115‑1 (MQ=255)                                    tCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGaa                                                                                                                                 >  5:240999/1‑115 (MQ=255)                                    tCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCaa                                                                                                         <  6:226127/139‑1 (MQ=255)                                    tCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCaa                                                                                                         <  3:148174/139‑1 (MQ=255)                                         cacaTCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGa                                                                                                    <  1:218190/139‑1 (MQ=255)                                         cacaTCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGa                                                                                                    <  4:280263/139‑1 (MQ=255)                                            atcatTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGaa                                                                                                                                 >  7:121272/1‑107 (MQ=255)                                            atcatTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGaa                                                                                                                                 <  8:121272/107‑1 (MQ=255)                                            atcatTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACat                                                                                                 <  6:273995/139‑1 (MQ=255)                                              catTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           <  7:78961/95‑1 (MQ=255)                                              catTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           >  8:78961/1‑95 (MQ=255)                                              catTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACatat                                                                                               >  7:42165/1‑139 (MQ=255)                                              catTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACatat                                                                                               >  2:26160/1‑139 (MQ=255)                                                      ccACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGCAAAAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTc                                                                                         >  5:165742/1‑139 (MQ=255)                                                      ccACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCa                                                                                                              >  8:108076/1‑116 (MQ=255)                                                      ccACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCa                                                                                                              <  7:108076/116‑1 (MQ=255)                                                      ccACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGa                                                                                                    <  8:120770/126‑1 (MQ=255)                                                      ccACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGa                                                                                                    >  7:120770/1‑126 (MQ=255)                                                      ccACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGCCAATTGGACATATGTGTTc                                                                                         <  2:223984/137‑1 (MQ=255)                                                       cACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGCAAAAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTcc                                                                                        <  2:31381/139‑1 (MQ=255)                                                          ttCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCATTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           <  3:215645/83‑1 (MQ=255)                                                          ttCAGGAAATTTTCAGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTcaagc‑‑a                                                                                                                                           >  4:215645/1‑83 (MQ=255)                                                            cAGGAAATTTTCCGAATAAATGATATAAAGGCATTGACCCTTTAAAAACAGCGTTGAGTGTCGTATCAAAATCCTCAAGC‑‑AAGCTCCAGTATACCTCGGTAACTTCCTCTGTCAATTGGACATATGTGTTCCGTTTa                                                                                   <  4:31312/137‑1 (MQ=255)                                                                       ccGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGt                                                                       <  4:245598/138‑1 (MQ=255)                                                                           aTAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  <  1:159536/139‑1 (MQ=255)                                                                           aTAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  <  6:24518/139‑1 (MQ=255)                                                                           aTAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  <  6:15549/139‑1 (MQ=255)                                                                            tAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGgag                                                                 <  7:92418/139‑1 (MQ=255)                                                                                    ataAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTa                                                          >  8:10340/1‑138 (MQ=255)                                                                                        aGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  >  6:214667/1‑126 (MQ=255)                                                                                        aGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  <  5:214667/126‑1 (MQ=255)                                                                                             gtgACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATcccc                                                 >  7:49530/1‑138 (MQ=255)                                                                                                    tttAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcc                                         >  3:193754‑M1/1‑131 (MQ=255)                                                                                                    tttAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcc                                         >  5:63733‑M1/1‑131 (MQ=255)                                                                                                        aaaaCAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgat                                     <  7:10340‑M1/139‑13 (MQ=255)                                                                                                                 ttGCGTGTCGGATCAAATTCCTCAAGCAAAAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCGTTCTTCCAGCTTTTTTGa                              <  6:165742/139‑1 (MQ=255)                                                                                                                           gATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgatgacgagtacattttgcacc                  >  2:52953‑M1/1‑108 (MQ=255)                                                                                                                           gATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgatgacgagtacattttgcacc                  >  3:236132‑M1/1‑108 (MQ=255)                                                                                                                              cAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgatgacgagtacattttgcacccgc               <  8:49530‑M1/139‑35 (MQ=255)                                                                                                                                aaTTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  >  2:240714/1‑86 (MQ=255)                                                                                                                                aaTTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  >  6:115577/1‑86 (MQ=255)                                                                                                                                aaTTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  <  5:115577/86‑1 (MQ=255)                                                                                                                                aaTTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGga                                                                  <  1:240714/86‑1 (MQ=255)                                                                                                                                    ccTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgatgacgagtacattttgcacccgcgggaa          <  4:193754‑M1/138‑40 (MQ=255)                                                                                                                                           c‑‑aagcTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgatgacgagtacattttgcacccgcgggaacagcagaa  <  5:6047‑M1/139‑48 (MQ=255)                                                                                                                                           c‑‑aagcTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCcacgatcccgatgacgagtacattttgcacccgcgggaacagcagaa  <  5:182503‑M1/139‑48 (MQ=255)                                                                                                                                             |                                                                                                                                           AATTTCCATAAAATCATCGCCGTAGATATGGCGTATCATACACATCATTTCTGCCACTTCAGGAAATTTTCCGAATAAATGATATAAAGGCTGTGACCCTTTAAAAACAGCGTTGCGTGTCGGATCAAATTCCTCAAGC‑‑AAGCTCCAGAATACCTCGGTACCTTCCTCAGTCAATTGGACATATGTGTTCCGTTTATCATTCAGCCGTTTGGAGAACCTTAAATATCCCCGTTCTTCCAGCTTTTTTGAAAAGTTGAATGCGGTTGATACGTGCATG  >  BS168_reference_genome/1073229‑1073505

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