Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 3,770,066 (A)8→9 intergenic (‑258/+38) ywzF ← / ← csbD conserved protein of unknown function/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,770,0581.A93.8% 190.1 / 6.4 64intergenic (‑250/+46)ywzF/csbDconserved protein of unknown function/stress response protein
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base A (36/24);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

ATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGC  >  BS168_reference_genome/3769929‑3770187
                                                                                                                                   |                                                                                                                                    
ataAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATAAATTCAAAGAATATTTAAAAAATTTAGTTACATTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑‑aaaaaaa                                                                                                                              >  6:326928/1‑137 (MQ=255)
ataAAAGATTACAACATATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑‑aaaaaaaa                                                                                                                            <  7:291616/138‑1 (MQ=255)
      gTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAATAGCG‑‑aaaaaaaaagac                                                                                                                        >  3:64420/1‑132 (MQ=255)
      gTTTACAACTTATGTCATCATTATATAGAAGAGGCCGATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGAACAAATCAGCTTCTTTTTTTGTTATTTTTATTCACATCAGATAAAGAG‑‑aaaaaaaaagac                                                                                                                        >  7:227199/1‑132 (MQ=255)
        ttACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGa                                                                                                                    >  7:47685/1‑138 (MQ=255)
         tACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGaaa                                                                                                                  >  4:21974/1‑139 (MQ=255)
         tACAACTTATGACATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGATACCTTGACCAAATCAGCATCTTTTTTTGTTATATTTATTAACATCAGAAAAAGCGAAAAAAAAAGACAAGaaa                                                                                                                  >  8:63164/1‑139 (MQ=255)
           cAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGaaaaa                                                                                                                >  4:23550/1‑137 (MQ=255)
                 aTGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAGTTCGt                                                                                                          >  5:295424/1‑139 (MQ=255)
                    tgatgatTATATTGAAGAGGCATATTTATTCTAAGAATATTTAAAAAATTTTGTTAACTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGTAAAATCGAAAAAAAAAGACAAGAAAGATTTCGTTTc                                                                                                       <  8:78544/134‑1 (MQ=255)
                    tcatcaTTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTc                                                                                                       <  1:3001/139‑1 (MQ=255)
                      atcatTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAAGACAAGAAACATTTCGTTTcc                                                                                                      <  6:152658/139‑1 (MQ=255)
                      atcatTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTc                                                                                                       >  6:343022/1‑137 (MQ=255)
                      atcatTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTc                                                                                                       <  5:343022/137‑1 (MQ=255)
                        catTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGc                                                                                                   >  4:166262/1‑139 (MQ=255)
                        catTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTAACATAAGAAAAATCTA‑‑AAAAAAAAGACAAGAAACATTACGTTTCCTGc                                                                                                   >  2:276610/1‑139 (MQ=255)
                         aatatctaGAAGAGGCTTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAATACAAGAAACATTTCGTTtgctgc                                                                                                   <  6:112415/136‑1 (MQ=255)
                                  aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAACACGCGAAAAAAAAAGCCAAGAAGCATTTCGTTGCCGGCGTCTTCCGTa                                                                                         >  2:142916/1‑139 (MQ=255)
                                  aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCg                                                                                           >  6:224329/1‑137 (MQ=255)
                                  aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCg                                                                                           <  5:224329/137‑1 (MQ=255)
                                  aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                         >  2:189638/1‑139 (MQ=255)
                                  aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                         >  7:175644/1‑139 (MQ=255)
                                  aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                         >  2:89192/1‑139 (MQ=255)
                                    gagGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCTA‑‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTAc                                                                                        >  6:53439/1‑138 (MQ=255)
                                    gagGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTAc                                                                                        >  8:33660/1‑138 (MQ=255)
                                     atgacTATTAATTCAAAGAATATTTAAAAAATTTAGTTAAATTTAAAAAAAAAGATTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTAGTTTCCTGACTCTTCCTTACGa                                                                                      >  8:89063/5‑139 (MQ=255)
                                      ggCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCGGCGTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  1:246491/1‑139 (MQ=255)
                                       gcatATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  7:210098/4‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAACCATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  4:16654/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGag                                                                                     >  5:19540/1‑137 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  1:185750/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  6:79142/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  2:17863/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  7:29120/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  3:223815/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  8:124598/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  3:313149/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  4:12112/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                     >  5:285398/1‑138 (MQ=255)
                                       gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAtt                                                                                    >  5:101533/1‑138 (MQ=255)
                                            ttAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCt                                                                               >  5:329096/1‑139 (MQ=255)
                                            ttAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑‑AAAAAAACAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTATTCt                                                                               >  4:36007/1‑139 (MQ=255)
                                              aaTTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCtt                                                                              <  1:343235/138‑1 (MQ=255)
                                                  cAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCTCATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGc                                                                                                   <  1:112355/113‑1 (MQ=255)
                                                  cAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCTCATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGc                                                                                                   >  2:112355/1‑113 (MQ=255)
                                                          atTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGttt                                                                                                        <  6:78870/100‑1 (MQ=255)
                                                          atTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGttt                                                                                                        >  5:78870/1‑100 (MQ=255)
                                                                    tttAGTTACCTTGACCAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCtt                                                      <  6:101533/139‑1 (MQ=255)
                                                                     ttAGTTACCTTGAGCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCATAAAAAGCGAAAAAAAAATACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCt                                                       <  8:79019/138‑1 (MQ=255)
                                                                        gTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTcc                                                   <  6:286042/139‑1 (MQ=255)
                                                                             cTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGata                                              <  6:329096/139‑1 (MQ=255)
                                                                             cTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGata                                              <  7:89063/139‑1 (MQ=255)
                                                                             cTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGata                                              <  7:124598/139‑1 (MQ=255)
                                                                                 aCCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCt                                          >  8:62260/1‑139 (MQ=255)
                                                                                 aCCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCATCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTTGATTTCCCGATATcat                                          >  8:50572/1‑137 (MQ=255)
                                                                                          gCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGAttt                                 <  3:166262/139‑1 (MQ=255)
                                                                                          gCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGAttt                                 <  8:29120/139‑1 (MQ=255)
                                                                                                 tttttGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCtt                           <  7:7182/138‑1 (MQ=255)
                                                                                                 tttttGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCt                            <  8:210098/137‑1 (MQ=255)
                                                                                                                ttCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCttttt            <  2:33257/138‑1 (MQ=255)
                                                                                                                    caTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCt       >  2:274220/1‑139 (MQ=255)
                                                                                                                    caTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCt       >  7:158330/1‑139 (MQ=255)
                                                                                                                     aTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTAt                  <  2:215156/127‑1 (MQ=255)
                                                                                                                     aTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTAt                  >  1:215156/1‑127 (MQ=255)
                                                                                                                     aTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCtt      <  4:223815/139‑1 (MQ=255)
                                                                                                                      tCAGAAAAAGCG‑‑AAAAAAAATGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCtt      >  7:237561/1‑138 (MQ=255)
                                                                                                                        aGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCa    <  2:185750/138‑1 (MQ=255)
                                                                                                                         gAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGc  <  3:23550/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
ATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGC  >  BS168_reference_genome/3769929‑3770187

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: