Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,059,485 T→G 9.6% intergenic (‑20/+3) yxeN ← / ← yxeM putative ABC transporter (permease)/putative ABC transporter (binding lipoprotein)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,059,4850TG 9.6% 171.2 / 2.0 65intergenic (‑20/+3)yxeN/yxeMputative ABC transporter (permease)/putative ABC transporter (binding lipoprotein)
Reads supporting (aligned to +/- strand):  ref base T (12/46);  new base G (5/2);  total (17/48)
Fisher's exact test for biased strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTCCTCCGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATCTAGGGCCTTATTGACTTTTTTGCGGAGCTTGTCGTGCGCATCGTCCT  >  BS168_reference_genome/4059348‑4059618
                                                                                                                                         |                                                                                                                                      
tAAGTCCTCCGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGttt                                                                                                                                       <  8:130035/139‑1 (MQ=255)
tAAGTCCTCCGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGttt                                                                                                                                       >  2:71827/1‑139 (MQ=255)
     cctccGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTaa                                                                                                                                  <  1:297027/139‑1 (MQ=255)
     cctccGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTaa                                                                                                                                  <  4:391768/139‑1 (MQ=255)
     cctccGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTaa                                                                                                                                  <  6:49987/139‑1 (MQ=255)
      ctccGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAAt                                                                                                                                 <  6:360966/139‑1 (MQ=255)
      ctccGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTGTCTTAAt                                                                                                                                 >  4:19286/1‑139 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  3:382441/139‑1 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  5:76160/139‑1 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  2:41188/139‑1 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  1:71827/139‑1 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  1:336787/139‑1 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  7:390359/139‑1 (MQ=255)
          gATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  8:234620/139‑1 (MQ=255)
           aTAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             >  2:312731/1‑138 (MQ=255)
               tcgCAAAAATCATAGCAGCTATTGCCATAAACAACGTGATCGGAAGGGCCTGAATTAAAGTCGTGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGtt                                                                                                                        <  1:373705/137‑1 (MQ=255)
               tGGCAAAAATCATAGGAGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGtt                                                                                                                        <  6:398298/139‑1 (MQ=255)
               tGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGtt                                                                                                                        <  6:98616/139‑1 (MQ=255)
               tGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGtt                                                                                                                        <  4:97432/139‑1 (MQ=255)
                    aaaaTCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACt                                                                                                                   <  4:312232/139‑1 (MQ=255)
                          aTAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGat                                                                                                              <  1:69507/138‑1 (MQ=255)
                           tAGCTGCTATTGCCATAATGAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGatat                                                                                                            <  7:263503/139‑1 (MQ=255)
                           tAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTAatat                                                                                                            <  3:19286/139‑1 (MQ=255)
                           tAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTGCTACTGTGatat                                                                                                            <  8:58072/139‑1 (MQ=255)
                              ctgctATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGATTCTGTTCTACTGTGATATCtt                                                                                                         <  2:57533/139‑1 (MQ=255)
                                 ctATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCATCTGTTCTACTGTGATATCTTCAt                                                                                                      >  7:11791/1‑139 (MQ=255)
                                    ttGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTa                                                                                                    <  5:348634/138‑1 (MQ=255)
                                    ttGCCATAAACAAGGTGAGCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTa                                                                                                    <  5:342948/138‑1 (MQ=255)
                                                 ggTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGGTTGTGTTATACTGTGATATATTAATTATAGTAtttttaa                                                                                        >  1:308990/1‑135 (MQ=255)
                                                 ggTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             >  7:194942/1‑100 (MQ=255)
                                                 ggTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCtt                                                                                                                             <  8:194942/100‑1 (MQ=255)
                                                  gTGCTCGGAAGGGCCGGCATTAAACTCGGGCACGCTGATATCATGAATGCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGGGGCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCagag                                                                                     <  8:416276/139‑1 (MQ=255)
                                                  gTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCagag                                                                                     <  8:194186/139‑1 (MQ=255)
                                                  gTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCagag                                                                                     <  3:103979/139‑1 (MQ=255)
                                                  gTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTAATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCagag                                                                                     <  5:176618/139‑1 (MQ=255)
                                                          aaGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGtttttt                                                                             <  6:332414/139‑1 (MQ=255)
                                                            gggCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCa                                                                           <  7:58792/139‑1 (MQ=255)
                                                            gggCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTACCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCa                                                                           <  4:203642/139‑1 (MQ=255)
                                                             ggCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATTCTTATGTTAAACGGTGATATCTTCATTAAAGTTTTTTTATGAGAGTTTTTTGAg                                                                          >  3:373399/1‑139 (MQ=255)
                                                               ccTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAgtg                                                                        >  6:137176/1‑139 (MQ=255)
                                                                   aaTTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTc                                                                     <  1:396473/138‑1 (MQ=255)
                                                                   aaTTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTc                                                                     <  8:11791/138‑1 (MQ=255)
                                                                   aaTTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTc                                                                     <  3:399445/138‑1 (MQ=255)
                                                                        aaGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTGTGTTTATTGTTCTGTTGTAATGTGTTATAGTGATTAAATTATTTTTAAGAGTATTTTTTGATTGTTGCGTCtt                                                               >  4:38166/1‑139 (MQ=255)
                                                                        aaGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTGTCTTAATGCGTATGTTCTACGTTGGATTCTTGATTTAAGTATTTTTCAGAGAGTGTTTTGAGTTTTACGTCtt                                                               >  1:120030/1‑139 (MQ=255)
                                                                        aaGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGGGTCTTAATGCTTGGGGTCTTCTGTGATTTAGTCATTAAATTATTTTTCAGAGAGTTTTTTAAGTTTTGCGTCtg                                                               >  2:370193/1‑138 (MQ=255)
                                                                             gggAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTATTAATGCTTCTGTTCTACTGTGATATATTAATTAAATTATTTTTAATAGATTTTTTTAAGTGTTACGTATTTTCg                                                           >  6:111472/1‑138 (MQ=255)
                                                                                aaCGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAAt                                                       <  1:165492/139‑1 (MQ=255)
                                                                                aaCGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAAt                                                       <  1:312731/139‑1 (MQ=255)
                                                                                aaCGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAAt                                                       <  1:380781/139‑1 (MQ=255)
                                                                                aaCGCCGATATCATGAATTCCCAATCAATCGTGTTCAATCTGGTGAGTCACCGCCGTGTCTTAATGCTGCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAAt                                                       <  7:334372/139‑1 (MQ=255)
                                                                                 aCGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAAtt                                                      <  4:149277/139‑1 (MQ=255)
                                                                                 aCGCTGATATCATGAATTCCCAATCAATCGTGGTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAAtt                                                      <  6:14986/139‑1 (MQ=255)
                                                                                     tGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTTATGGTTCTGTTGTACTGTGTTATGTTAATTAAATTATTTTTAAGAGAGGTTTTTGAGTGTTGCGTGGTTTCGCAATTGAtg                                                  >  7:421730/1‑138 (MQ=255)
                                                                                       atatCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATTCTTCTGTTATACTGTGATATATTAATTAAAGTATTTTTCAGAGAGTTATTTAAGTGTTACGTCTTTAAGCAAttaaac                                                  >  3:167018/1‑133 (MQ=255)
                                                                                       atatCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATc                                                  <  5:167943/137‑1 (MQ=255)
                                                                                       atatCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATc                                                  <  5:44862/137‑1 (MQ=255)
                                                                                       atatCATGAATTACCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGATCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATc                                                  <  3:255912/137‑1 (MQ=255)
                                                                                                 ttCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATCTAGGGCCTTAtt                                     <  5:137176/139‑1 (MQ=255)
                                                                                                              gtgtTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGAAATTCATCTAGGGCCTTATTGACTTTTTTGCgt                        <  8:175750/139‑2 (MQ=255)
                                                                                                                 ttCATTCTGGTGAGTCACCGCCGTTTATTAATGCTTGTGTTCTACTGTGATAAATTCAATAAAGTATTTTTAAGAGAGTTTTTTAATTGTTCCGGATTTAAGGAATTCATCTTGGGCCGTTTTTGATTTTTTTCGGAGc                     >  3:101976/1‑139 (MQ=255)
                                                                                                                      tCTGGTGAGTCACCGCCGTTTCTTAATGCTTATGTTCTACTGTGATATATTCATTAAATTATTTTTAAGAGAGTTTTTTCATTTTTGAGTCTTTACGCAATTAATCTAGTGCCTTATTTACTTTTTTTGGGAGCTTGTc                >  3:84691/1‑139 (MQ=255)
                                                                                                                         ggTGAGTCACCGCCGTTTCTTAATGCTTCTGTTGTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTAAGTGTTGCGTGTTTTCGCAATTCAAATAGGGCCTTTTTTACTTTTTTTCGGCGGTTgtggtg             >  6:300886/1‑139 (MQ=255)
                                                                                                                          gTGAGTCACCGCCGTGTAGTAATGCTTCTGTTCTAATTTGATATGTTCATTAAAGTATTTTTAAGAGAGTTTTTTCATTGTTGCGTCTTTACGCAATTAATCTAGGGCCTTATTGACTTTTTTGCGGAGCTTGTAGTgc            >  3:201123/1‑139 (MQ=255)
                                                                                                                             aGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATCTAGGGCCTTATTGACTTTTTTGCGGAGCTTGTCGTGCGCa         <  7:366464/139‑1 (MQ=255)
                                                                                                                             aGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATCTAGGGCCTTATTGACTTTTTTGCGGAGCTTGTCGTGCGCa         <  4:373399/139‑1 (MQ=255)
                                                                                                                                    gccgTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATCTAGGGCCTTATTGACTTTTTTGCGGAGCTTGTCGTGCGCATCGTCCt  <  5:354776/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TAAGTCCTCCGATAATGGCAAAAATCATAGCTGCTATTGCCATAAACAAGGTGATCGGAAGGGCCTGAATTAAAGTCGGGAACGCTGATATCATGAATTCCCAATCAATCGTGTTCATTCTGGTGAGTCACCGCCGTTTCTTAATGCTTCTGTTCTACTGTGATATCTTCATTAAAGTATTTTTCAGAGAGTTTTTTCAGTGTTCCGTCTTTACGCAATTCATCTAGGGCCTTATTGACTTTTTTGCGGAGCTTGTCGTGCGCATCGTCCT  >  BS168_reference_genome/4059348‑4059618

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: