Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	4,087,723	G→C	14.2%	G90R (GGC→CGC)	yxcD →	conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	4,087,723	0	G	C	14.2%	114.8 / 5.6	68	G90R (GGC→CGC)	yxcD	conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base C (2/6);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAACTCAGCGGAT  >  BS168_reference_genome/4087587‑4087859                                                                                                                                         |                                                                                                                                          tGAAGCTACATTTGGCCTGCCCCACAAGCAGCTTATGTCAGATGACCTTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTATCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGc                                                                                                                                         >  4:134318/1‑139 (MQ=255) tGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGc                                                                                                                                         >  8:48849/1‑139 (MQ=255) tGAAGCTACATTTGGCCTGAACCACAAGCAGCTTATGTCAGATGAAATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGTACGTAACCGGCGTGGAATTTCAATTTTATAAGAAAAAAgcc                                                                                                                                         >  2:54388/1‑137 (MQ=255)   aaGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCtt                                                                                                                                       >  1:160786/1‑139 (MQ=255)   aaGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCtt                                                                                                                                       >  6:61734/1‑139 (MQ=255)            ttGGCCTGCACCACAAGCAGCTTATGTCATATGACATTTCAGAAGGCATTATTCAGTTCTTAGTGTAATACCACAAGTTCATTCCGGAGTTTACCGTCTTTTATTTGCATTTTGATAAGAAAAAAGGCTTCTCCGCC‑t                                                                                                                              <  1:40402/138‑1 (MQ=255)               gCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGTATTGCAATTTTATTAGAAAAAAGGCATCTCccccggggg                                                                                                                           >  5:43538/1‑134 (MQ=255)                   gCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCATAGAGCAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTTATAAGAAAAAAGGCATATCCGCC‑TTAGtattt                                                                                                                      >  8:13414/1‑139 (MQ=255)                    caccacAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGGAATTTGATAATAAAAAAGGCTTCTCCGCC‑TTGGTTTTTg                                                                                                                     >  3:35003/1‑139 (MQ=255)                        acaAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCTTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCaa                                                                                                                 <  6:207750/139‑1 (MQ=255)                        acaAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTTAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCaa                                                                                                                 <  1:165874/139‑1 (MQ=255)                        acaAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCaa                                                                                                                 <  3:199059/139‑1 (MQ=255)                              agcTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAAGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAg                                                                                                            >  6:129506/1‑138 (MQ=255)                                cTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCag                                                                                                         <  5:204519/139‑1 (MQ=255)                                    tGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTCCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCagaaga                                                                                                     <  8:30887/139‑1 (MQ=255)                                          attaCATTTCAGAAGTCATTATTCAGTTCTTACAGTAATACTACAAATTCAATCCGGACGTTACCGTCGTTTAATTGCAATTTTAAACCAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATaaaa                                                                                               <  5:31909/136‑1 (MQ=255)                                          aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGTATTTGAATTTTATTATAAAAAAGGCTTGTCCGCCGTTGGGTTTTG‑GAAAGAAGGCGAAGAATaaa                                                                                                >  1:241993/1‑138 (MQ=255)                                          aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGCATTTCAAATTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGGTTTTTTCAACCAAGCAGCAGAATaaa                                                                                                >  1:168475/1‑138 (MQ=255)                                          aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTGCCCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATaaa                                                                                                >  2:257268/1‑138 (MQ=255)                                          aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTGGTATTGCAAATTGTTAAGAAAAAAGGCTTCTCCGCC‑GTGGTTTTTGTCAACGCAGCAGAAGAATaaa                                                                                                >  3:206865/1‑138 (MQ=255)                                          aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCCTCGTGGAATTGCAATTTTCAAAGTAAAAAGGGTTCGCCGCC‑ATGGGTTTTGTCAAAGAAGGAGCAGAATtaa                                                                                                >  1:65734/1‑138 (MQ=255)                                           ttaCATTTAATATGGCATTATTCAGTTCTTCGAGGAATAATACAACTTCAATCCGGACGTTACCGTCGTTGAATTCCATTTTGATAACAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATaaaa                                                                                               <  5:35904/136‑1 (MQ=255)                                              cATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTTCAAATTGATAAGAAAAAAGGCGTGTCCGCC‑ATGGTTTTTTTCTAAGAAACAGATAAATAAAAAtt                                                                                            >  2:251118/1‑137 (MQ=255)                                                     taaGCCATTAATCAGTTCTATCCGGAATACCACAACATCACACCGCACGTTACCGCGGTTTAATTCCAATTTCTAAATAAAAAAGGCTTCTCCGCC‑TTGGATTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  4:262517/138‑1 (MQ=255)                                                     gAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  7:11301/139‑1 (MQ=255)                                                     gAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAACGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  2:160786/139‑1 (MQ=255)                                                     gAAGGCATTATTAAGTTCTTAGAGGAATCCCACAACTTCAATCCGGCCGTAACCGTCGTTGAATTGCAATTTGATAACAAAAAACGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  2:107039/139‑1 (MQ=255)                                                     gAAGCCATTAATCAGTTATTACACTAAAACCACAACTTCAATCCGCACGTTACCGGCGTTTAATTGAAATTTAATAAAAAAAACCGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  3:15848/139‑1 (MQ=255)                                                     catgtcATTATTCTTTTATTAGTGAAACACCACAACTTAAACCGGGCCGTTCCGGTTGTTTATTTGAATTTTAAAAACAAAAACCGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  6:43538/134‑1 (MQ=255)                                                      aaGGCATTATTCAGTTCTAAGAGCAATTCAACAACTTCAACACGGACGTAACCGTCGTTTAATTCAAAATTCATAAGAAAACCCGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCa                                                                                   <  7:104742/139‑1 (MQ=255)                                                           attattCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCGCCGCC‑GTGGTTTTTGTAAACGAAGCAGAACAATAAAAATGTCTTAATCAAACCt                                                                              >  2:186047/1‑139 (MQ=255)                                                           attattCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAAATTGATAAAAAAAAAGGCTTCGCCGCC‑ATGGTTTTTGTCAACGAAGCACAATAAAAAAAAAGTATTAATCaaaac                                                                               >  4:79033/1‑136 (MQ=255)                                                              attCAGTTCTTTGAGGAATACCACAAGTACAAGGCGGCCTTGACCGTCTTTTAATGGATATTTGATAAGAATAAAGGCTTCTCCGCC‑TTGGTTTTTGTCCACGAAGCAGAAGAATAAAAATGTCTTAATCAAACATTa                                                                            <  8:69166/138‑1 (MQ=255)                                                              attCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTa                                                                            <  2:12878/138‑1 (MQ=255)                                                                  aGTTCTTAGAGTAATACCACAACTTCAATCCGCACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcg                                                                       <  3:109245/139‑1 (MQ=255)                                                                  aGTTATTAGAGGAATACCACAACTTCAATCCGGACGTTACCGGCTTTTAATTGCAATTTGAAAAGAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcg                                                                       <  4:35003/139‑1 (MQ=255)                                                                   gTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAATGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    >  4:41629/1‑125 (MQ=255)                                                                   gTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                                                                                    <  3:41629/125‑1 (MQ=255)                                                                    ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGGATTTCAATTTGATAAGAAAAAACGCTTATCCGCC‑CTGGGTTTTGTCAACGAAGCACAAGAATAAAAATTTCTTAAAAAAACCTTACTCcccg                                                                     >  1:128947/1‑139 (MQ=255)                                                                    ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTTCAATTTGATAAGAAAAAAGGCTTCTCCCCC‑TTGGTTTTTTTGAACGAAGCAGAAGAATAAAAATGTATTAATAAAACCTTACTAcgcg                                                                     >  8:124283/1‑139 (MQ=255)                                                                    ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcgcg                                                                     >  2:41367/1‑139 (MQ=255)                                                                    ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcgcg                                                                     >  5:109001/1‑139 (MQ=255)                                                                    ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcgcg                                                                     >  8:198790/1‑139 (MQ=255)                                                                           agGATTACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATACGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTaa                                                              <  3:79033/139‑1 (MQ=255)                                                                              agtACCACAACTTCAATCCGGAGGTTACGGTCGTTCAATTCCAATTTGATAAGAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGt                                                           <  3:229813/137‑1 (MQ=255)                                                                              aaTACCACACCTACAATCCCGACGTTACCGTCGTTTAATTGCAATTTGATAACAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCCACGACGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGt                                                           <  5:129506/139‑1 (MQ=255)                                                                              aaTACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTTATAAGAAAAAACGCATCTCCGCC‑TTGGTTTTTGTCAACGAAGCATAATAATAAAAATGTCTTAATAAAACCAAACACCGCGCGGGTAATGt                                                           >  2:27994/1‑139 (MQ=255)                                                                                tACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGttt                                                         >  1:91558/1‑139 (MQ=255)                                                                                      aaCTACAACCCGCACGTTACCGTCTTTTAATTGAAATTTGATAAGAAAAAACGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTaa                                                   <  1:257268/139‑1 (MQ=255)                                                                                        cTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTTTCAACGAAGCAGAAGAATAAAAATGTCTTAATAAAACCTTACACCGCGCGGGGAAGGTTTTTTTAAtt                                                 >  1:225862/1‑138 (MQ=255)                                                                                           cAATCCGGACGTTACCGTCGTGGAATTGCAATTTGATTATAAAAAAAGCTTCACCCCC‑TTGGTTTTTTTAAACGAAGCAGAATAAAAAAAATTTCTTAATAAAAACATACACCGCCCCGGGAATGTTTTTTTAAttggt                                              >  4:229813/1‑135 (MQ=255)                                                                                             aTCCGGACGTTACCGTCGTTGAATTGGAAATTGATAAGAAAAAAAGGTTCTCCGCC‑GGGGTTTTTGTCAACGCAGCAGAATAATAAAAATTTGTTAATCAAACCTTACTCCGCGCGGGGAAGCTTGTTTTAATGGttgg                                            >  4:175143/1‑137 (MQ=255)                                                                                                cGGACGTTACCGTCGTTGAATTCCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTctct                                         <  6:138805/139‑1 (MQ=255)                                                                                                cGGAAGTTACCGTCGTTTAATTGAAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTctct                                         <  1:27994/139‑1 (MQ=255)                                                                                                  gACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTctct                                         >  1:50794/1‑137 (MQ=255)                                                                                                  gACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTctct                                         >  8:183291/1‑137 (MQ=255)                                                                                                        aCCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTATACGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATa                                  >  8:72676/1‑138 (MQ=255)                                                                                                            tcgtTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGtt                             >  2:9725/1‑139 (MQ=255)                                                                                                            tcgtTGAATTGCAATTTGATAAGAAAAAAGGCGTCTACGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGCGCATAAAGtt                             >  1:206275/1‑139 (MQ=255)                                                                                                               ttGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCGTAAAGTTAtt                          >  8:86418/1‑139 (MQ=255)                                                                                                                   aTTGCATTTTGATAAGAAAATAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGt                      <  7:13414/139‑1 (MQ=255)                                                                                                                          tttCATAAAAAAAACGGATTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAAAAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCt               <  4:203761/139‑1 (MQ=255)                                                                                                                              aTAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCCACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAAc           >  2:79824/1‑139 (MQ=255)                                                                                                                               tAAGAAAAACGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAACt          <  1:251118/139‑1 (MQ=255)                                                                                                                                 aGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACt                                                                          >  8:219791/1‑73 (MQ=255)                                                                                                                                 aGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACt                                                                          <  7:219791/73‑1 (MQ=255)                                                                                                                                       aaGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAACTCAGCGGAt  <  6:109001/139‑1 (MQ=255)                                                                                                                                        aGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAACTCAg       <  3:244464/133‑1 (MQ=255)                                                                                                                                        aGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAACTCAg       >  4:244464/1‑133 (MQ=255)                                                                                                                                         |                                                                                                                                          TGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCC‑TTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTTTTTTAATGGTTTCTCTGAGCATAAAGTTATTTCGTTTTTCCTAAACTCAGCGGAT  >  BS168_reference_genome/4087587‑4087859

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