Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,364,921	A→G	7.5%	L58L (TTG→CTG)	qcrA ←	menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,364,921	0	A	G	7.5%	194.6 / 7.2	67	L58L (TTG→CTG)	qcrA	menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base G (2/3);  total (29/38)
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.50e-02

CTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAGCCGCCATAAAACCTCCTACGCCTGTGAGCGTATAATTCAAA  >  BS168_reference_genome/2364783‑2365057                                                                                                                                           |                                                                                                                                         cTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCaa                                                                                                                                          >  7:301245/1‑139 (MQ=255)   aaaTGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACt                                                                                                                                        >  6:379840/1‑139 (MQ=255)   aaaTGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACt                                                                                                                                        >  3:10919/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGCCCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAAc                                                                                                                                 >  8:255708/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTTACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAAc                                                                                                                                 >  7:139065/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAAc                                                                                                                                 >  6:278050/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAAc                                                                                                                                 >  7:371463/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAAc                                                                                                                                 >  5:216059/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCGCTTTTGTCAACTCATCAAc                                                                                                                                 >  4:9758/1‑139 (MQ=255)          tACAAATTGGCGAAAGTGCTACAATATCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATAAATATTAAAATCAAAGAGCTGAGGCTCTTTTGTCAACTCATCAAc                                                                                                                                 >  4:107244/1‑139 (MQ=255)           aCAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCaa                                                                                                                                  <  2:358824/137‑1 (MQ=255)              aaTTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTg                                                                                                                             <  5:379840/139‑1 (MQ=255)               aTTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGa                                                                                                                            <  3:270118/139‑1 (MQ=255)                   gCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAAc                                                                                                                        <  8:301245/139‑1 (MQ=255)                     gAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCt                                                                                                                      <  5:5542/139‑1 (MQ=255)                     gAAAGTGCTACAATTTCATCCCCGTTGTTGAAGACCCAGGCTGATCTTGCCTCATCTGACTCATACCAGGCATCGACATGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACACCCt                                                                                                                      <  4:37689/139‑1 (MQ=255)                        aGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGaa                                                                                                                   <  8:380731/139‑1 (MQ=255)                                            gTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTc                                                                                               >  4:213495/1‑139 (MQ=255)                                            gTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTc                                                                                               <  2:178027/139‑1 (MQ=255)                                            gTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTc                                                                                               <  6:216059/139‑1 (MQ=255)                                            gTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTc                                                                                               <  3:310128/139‑1 (MQ=255)                                            gTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTc                                                                                               <  4:104902/139‑1 (MQ=255)                                            gTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCAGAGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTc                                                                                               <  7:186490/139‑1 (MQ=255)                                                            gCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGt                                                                               <  8:301123/139‑1 (MQ=255)                                                             cTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTACCCTGTCGATTTTAATACCGGGTc                                                                              <  3:9758/139‑1 (MQ=255)                                                                aTCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGACCCATCTCCTGCTTCCCTGTCGATTTTTATAACGGGTCGAg                                                                           >  7:108957/1‑139 (MQ=255)                                                                          ttctGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  >  3:58425/1‑138 (MQ=255)                                                                            ctGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  7:243668/136‑1 (MQ=255)                                                                            ctGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  7:133909/136‑1 (MQ=255)                                                                            ctGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  7:83763/136‑1 (MQ=255)                                                                            ctGACTCATACCAGGCATCGACTTGATTAAATTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGGCGAGTgcgaagcg                                                                  <  5:310102/136‑1 (MQ=255)                                                                            ctGAATCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  4:10919/136‑1 (MQ=255)                                                                               aCTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGTGTCGAGTGCGAAGCGAACCa                                                             <  5:265711/138‑1 (MQ=255)                                                                               aCTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCa                                                             <  1:218638/138‑1 (MQ=255)                                                                                   aTACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGc                                                        <  5:278050/139‑1 (MQ=255)                                                                                    tACCAGGCATCGCCTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGCCAACCTGAACCATATCCTGCTTCCCTGTCGAGTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCa                                                       >  6:410917/1‑139 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACTCTGACAACCTGAACCAGATCCTGCTTCCCTGTCGATTTTAAAAACGGGTCGGGGGCGAAGCGAACCaatggg                                                        >  7:166794/1‑133 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGACCCATATCCTGCTTCCCTGTCGATTTTAAAACCGGGTCGAGTGCGAAGCGAACCATAGGCa                                                       >  4:394968/1‑139 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCata                                                                                                              >  8:274152/1‑84 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCata                                                                                                              <  7:274152/84‑1 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATAcc                                                                                   >  6:209816/1‑111 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATAcc                                                                                   <  5:209816/111‑1 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  6:209321/128‑1 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  >  6:198334/1‑128 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  >  5:209321/1‑128 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  5:198334/128‑1 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCa                                                       >  6:217050/1‑139 (MQ=255)                                                                                    tACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATAACGGGTCGATTTCGCAGCGAACCCTAGGc                                                        >  8:353366/1‑138 (MQ=255)                                                                                      ccAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCAt                                                      <  7:255708/138‑1 (MQ=255)                                                                                         ggCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  <  1:79008/123‑1 (MQ=255)                                                                                         ggCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaagcg                                                                  >  2:79008/1‑123 (MQ=255)                                                                                         ggCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTgcgaag                                                                    <  8:363383/121‑1 (MQ=255)                                                                                         ggCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACAGGGTCGAGTgcgaag                                                                    >  7:363383/1‑121 (MQ=255)                                                                                           cATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCACCGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGAGTTTAAGACCGGGGCGAGTGGGAAGCGCACCATACGCATGAGCa                                                 >  4:362648/1‑138 (MQ=255)                                                                                           cATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCa                                                             <  8:176764/126‑1 (MQ=255)                                                                                           cATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCa                                                             >  7:176764/1‑126 (MQ=255)                                                                                              cGACTTGGTTAATTTTGAAATCAAAGCGCTGCGGCTCTTTTGTCAGCTCATCTACGCTGACAACCTGAACCATGTCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCAAAGCGAAcc                                                              >  7:252226/1‑122 (MQ=255)                                                                                              cGACTTGGTTAATTTTGAAATCAAAGCGCTGCGGCTCTTTTGTCAGCTCATCTACGCTGACAACCTGAACCATGTCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTCCAAAGCGAAcc                                                              <  8:252226/122‑1 (MQ=255)                                                                                                 cTTGGTTAATTTTGAAATCAAAGCTCTGCTGCTCTTTTGTCAGCTCATCTACGATGACAACCTGAACCATGTCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTCCAAAGCGAACCATAGGCATGAGCATACTAgc                                          <  8:323484/139‑1 (MQ=255)                                                                                                 cTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAgc                                          <  1:399939/139‑1 (MQ=255)                                                                                                            ttAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAGCCGCCATaaaa                                <  4:41218/138‑1 (MQ=255)                                                                                                            ttAAAATCAAAGCGCGGCGGCTCTTTTGGCAACTCATCAACGCTGACAACCTGAACCATATCCTGCGTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAGCCGCCATaaaa                                <  6:147681/138‑1 (MQ=255)                                                                                                               aaaTCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGAttt                                                                                          >  4:210031/1‑77 (MQ=255)                                                                                                               aaaTCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGAttt                                                                                          <  3:210031/77‑1 (MQ=255)                                                                                                                  tCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAGCCGCCATAAAACCTCCTa                         <  8:139065/139‑1 (MQ=255)                                                                                                                    aaaGCGCTGCGGCTCTTTTGTCAGCTCATCTACGCTGACAACCTGAACCATGTCCTGCTTCCCTGTc                                                                                               <  8:210565/67‑1 (MQ=255)                                                                                                                    aaaGCGCTGCGGCTCTTTTGTCAGCTCATCTACGCTGACAACCTGAACCATGTCCTGCTTCCCTGTc                                                                                               >  7:210565/1‑67 (MQ=255)                                                                                                                           tGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAGCCGCCATAAAACCTCCTACGCCtgt                  >  7:226561/1‑137 (MQ=255)                                                                                                                                          aaCTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCAAGAGCATACTAGCCGCCATAAAACCTCCTACGCCTGTGAGCGTATAATTCaaa  >  8:181032/1‑138 (MQ=255)                                                                                                                                           |                                                                                                                                         CTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAATCAAAGCGCTGCGGCTCTTTTGTCAACTCATCAACGCTGACAACCTGAACCATATCCTGCTTCCCTGTCGATTTTAATACCGGGTCGAGTGCGAAGCGAACCATAGGCATGAGCATACTAGCCGCCATAAAACCTCCTACGCCTGTGAGCGTATAATTCAAA  >  BS168_reference_genome/2364783‑2365057

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