Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 67,130 G→T 5.2% E242D (GAG→GAT yabN → putative fusion methylase and nucleotide pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome67,1300GT 5.2% 310.0 / 2.8 90E242D (GAG→GATyabNputative fusion methylase and nucleotide pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base T (3/2);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATCAGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGACACAGACCATATGATCGAAGAGCTTGG  >  BS168_reference_genome/66993‑67267
                                                                                                                                         |                                                                                                                                         
gagaTCAGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGc                                                                                                                                          >  6:306143/1‑139 (MQ=255)
    tCAGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcg                                                                                                                                      <  5:282737/139‑1 (MQ=255)
    tCAGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcg                                                                                                                                      <  5:306143/139‑1 (MQ=255)
      cggACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  2:237978/138‑1 (MQ=255)
      aGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACGACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  5:297526/139‑1 (MQ=255)
      aGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  3:222138/139‑1 (MQ=255)
      aGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  1:79177/139‑1 (MQ=255)
         aCAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCTTCATCAGAGAGCTTCGCGGAc                                                                                                                                 <  6:14764/139‑1 (MQ=255)
           aGTGCCATTATTTGAGCTGGACCGCAATGTAGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGACTTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGAcc                                                                                                                                <  8:361139/138‑1 (MQ=255)
            gTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGcccg                                                                                                                               >  1:342471/1‑137 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGCGCTTCGCGGACCTaa                                                                                                                             >  1:173918/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGATCTTCGCGGACCTaa                                                                                                                             >  7:253498/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGATCTTCGCGGACCTaa                                                                                                                             >  5:70336/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  3:321212/132‑1 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  3:326489/132‑1 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    >  4:326489/1‑132 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    >  4:321212/1‑132 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGa                                                                                                                                  >  2:50618/1‑134 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGa                                                                                                                                  <  1:50618/134‑1 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTaa                                                                                                                             >  3:340256/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTaa                                                                                                                             >  7:286585/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTaa                                                                                                                             >  5:54979/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTaa                                                                                                                             >  6:287013/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTaa                                                                                                                             >  1:115894/1‑139 (MQ=255)
             tGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGATTCGCGGACCTaa                                                                                                                             >  8:193904/1‑139 (MQ=255)
                  ttattTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGAt                                                                                                                        >  2:67238/1‑139 (MQ=255)
                       tGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  4:27007/122‑1 (MQ=255)
                       tGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    >  3:27007/1‑122 (MQ=255)
                       tGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGa                                                                                                                                  <  6:54049/124‑1 (MQ=255)
                       tGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGa                                                                                                                                  >  5:54049/1‑124 (MQ=255)
                       tGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGAcc                                                                                                                                >  3:200848/1‑126 (MQ=255)
                       tGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGAcc                                                                                                                                <  4:200848/126‑1 (MQ=255)
                         aGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTg                                                                                                                 <  6:54979/139‑1 (MQ=255)
                              gACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAACCATCATCAGACACCTTCGCGGACCTAATGGATGTCCGGGGGag                                                                                                             >  1:365090/1‑137 (MQ=255)
                                           cTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCa                                                                                               >  5:143924/1‑139 (MQ=255)
                                               gAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCa                                                                                            <  5:346604/138‑1 (MQ=255)
                                               gAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCa                                                                                            <  2:173918/138‑1 (MQ=255)
                                                     cTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGCGCGCTTCGCGGACCTAATGGATGTCATTGGGATAAAAAGCATCCGAATCAGTCAtt                                                                                      >  4:194505/1‑138 (MQ=255)
                                                     cTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCAtt                                                                                      >  1:239652/1‑138 (MQ=255)
                                                     cTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCAtt                                                                                      >  1:359677/1‑138 (MQ=255)
                                                        aaCAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAg                                                                                  >  6:133583/1‑139 (MQ=255)
                                                         tcAAGTGTGTATATTCCGCCGAGTAAAGAAGAGATGCGGGTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAg                                                                                  <  2:365090/137‑1 (MQ=255)
                                                         acatttgtgtTTATTCCGGCTTTTAAATATGTTAAGATGCTTTATCAGTAATTCTCCACTTTCTTAAGCTTAATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAg                                                                                  <  2:342471/133‑1 (MQ=255)
                                                         aCAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAg                                                                                  <  2:269625/138‑1 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGATCTTAGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGAGGCATCAGTGATTAAAGCAGTACAt                                                                          >  1:348762/1‑139 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    <  4:267514/81‑1 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgcg                                                                                                                                    >  3:267514/1‑81 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgc                                                                                                                                     <  4:194708/80‑1 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTcgc                                                                                                                                     >  3:194708/1‑80 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTa                                                                             >  4:372596/1‑136 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTa                                                                             <  3:372596/136‑1 (MQ=255)
                                                                tgtATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACAt                                                                          >  5:20834/1‑139 (MQ=255)
                                                                      ttCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTga                                                                    <  2:115894/139‑1 (MQ=255)
                                                                       tCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCGAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTg                                                                     >  1:303641/1‑137 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAGAGGCAGCCGCCTCAGTCCTTACAGCAGTGCAGGAGTgggg                                                                  >  1:319389/1‑136 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAg                                                                                                    <  4:157308/105‑1 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAg                                                                                                    >  3:157308/1‑105 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACa                                                                           <  7:2441/130‑1 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACa                                                                           >  8:2441/1‑130 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTg                                                                     >  1:189155/1‑136 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTg                                                                     >  3:115929/1‑136 (MQ=255)
                                                                        ccgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTg                                                                     >  8:175911/1‑136 (MQ=255)
                                                                         cgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTgaaga                                                                 <  1:106817/139‑1 (MQ=255)
                                                                         cgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTgaaga                                                                 <  7:233887/139‑1 (MQ=255)
                                                                         cgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTgaaga                                                                 <  7:250161/139‑1 (MQ=255)
                                                                         cgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTgaaga                                                                 <  8:213816/139‑1 (MQ=255)
                                                                         cgccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTgaaga                                                                 <  3:77866/139‑1 (MQ=255)
                                                                          gccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGAt                                                                                                                        <  6:238504/83‑1 (MQ=255)
                                                                          gccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGAt                                                                                                                        >  5:238504/1‑83 (MQ=255)
                                                                          gccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCAt                                                                                       <  8:42235/116‑1 (MQ=255)
                                                                          gccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCAt                                                                                       >  7:42235/1‑116 (MQ=255)
                                                                           ccgATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAACAGCAGCCGCATCAGTCATTAAAGCAGTACATGATTGAAGAAt                                                               >  4:340940/1‑139 (MQ=255)
                                                                            cgATTAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGCCCTAATGGATGTGCGTGGGATAAAAAGCAGACGCATCAGTCATTAATCAAGTACATGATTGAAGAAAGt                                                             >  5:84490/1‑139 (MQ=255)
                                                                              aTTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGtt                                                            >  5:6480/1‑139 (MQ=255)
                                                                                        agaagCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTtct                                                  >  7:77022/1‑139 (MQ=255)
                                                                                                ctTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCg                                                <  1:87090/133‑1 (MQ=255)
                                                                                                ctTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCg                                                >  2:87090/1‑133 (MQ=255)
                                                                                                         gAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACgaag                                 <  1:204677/139‑1 (MQ=255)
                                                                                                         gAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACgaag                                 <  5:12069/139‑1 (MQ=255)
                                                                                                          aaTTCTCCACGTTCCGAAGCATCATCAGAGATCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACgaaga                                <  6:70336/139‑1 (MQ=255)
                                                                                                          aaTTCTCCACGTTCCGAAGCATCATAAGAGATCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGACTGTTATGAGCTTCTCGAAGCTATTGACgaaga                                <  8:253498/139‑1 (MQ=255)
                                                                                                                 cACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGaca                           <  5:133583/137‑1 (MQ=255)
                                                                                                                 cACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGaca                           <  1:67238/137‑1 (MQ=255)
                                                                                                                 cACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGAcaca                         <  5:287013/139‑1 (MQ=255)
                                                                                                                   cGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGACACAGa                       <  8:286585/139‑1 (MQ=255)
                                                                                                                         gAAGCATCATAAGAGAGCTTCGCGGACCTATTGGATGTCCGTGGGAAAAAAAGCGGACGCATCAGTCTTTTATGCAGTACATGATTGATGAATGTTATGAGCTTCTCGACGGTATTGACGCAGAAGACACAGACCttat                 >  6:168865/1‑139 (MQ=255)
                                                                                                                         caaccatcatCACCGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGACACAGACCatat                 <  1:362539/135‑1 (MQ=255)
                                                                                                                               tcatcaGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGACACAGACCATATGATCg            >  8:224046/1‑138 (MQ=255)
                                                                                                                                catcaGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAg                                                                                  >  7:144434/1‑67 (MQ=255)
                                                                                                                                catcaGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAg                                                                                  <  8:144434/67‑1 (MQ=255)
                                                                                                                                      agagCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGaca                           >  5:153675/1‑116 (MQ=255)
                                                                                                                                      agagCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGaca                           >  5:145888/1‑116 (MQ=255)
                                                                                                                                      agagCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGaca                           <  6:145888/116‑1 (MQ=255)
                                                                                                                                      agagCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGaca                           <  6:153675/116‑1 (MQ=255)
                                                                                                                                        agCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGACACAGACCATATGATCGAAGAGCTTgg  >  1:247837/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                         
GAGATCAGGACAGTGCCATTATTTGAGCTGGACCGCAATGTCGCTTTGAACAACTTAACAAGTGTGTATATTCCGCCGATTAAAGAAGAGAAGCTGCTTTATCACGAATTCTCCACGTTCCGAAGCATCATCAGAGAGCTTCGCGGACCTAATGGATGTCCGTGGGATAAAAAGCAGACGCATCAGTCATTAAAGCAGTACATGATTGAAGAATGTTATGAGCTTCTCGAAGCTATTGACGAAGAAGACACAGACCATATGATCGAAGAGCTTGG  >  BS168_reference_genome/66993‑67267

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: