Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 557,870 (T)5→6 intergenic (+421/‑3) yddT → / → ydzN putative exported phage protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome557,8651.T98.9% 331.7 / ‑3.4 93intergenic (+416/‑8)yddT/ydzNputative exported phage protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (44/48);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAA  >  BS168_reference_genome/557729‑557990
                                                                                                                                          |                                                                                                                             
tgAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGtt                                                                                                                             >  4:147647/1‑139 (MQ=255)
tgAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGtt                                                                                                                             >  2:308062/1‑139 (MQ=255)
  aaaGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGtttt                                                                                                                           >  1:358213/1‑139 (MQ=255)
      tAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGttttttaa                                                                                                                       <  1:367078/139‑4 (MQ=255)
        aTGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGttttttaatt                                                                                                                     >  4:144134/1‑134 (MQ=255)
        aTGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGttttttaatt                                                                                                                     >  7:456593/1‑134 (MQ=255)
         tGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTg                                                                                                                     >  5:4660/1‑139 (MQ=255)
                tAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTg                                                                                                                >  5:400061/1‑137 (MQ=255)
                     tCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTaaa                                                                                                         <  3:401137/139‑1 (MQ=255)
                     tCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTaaa                                                                                                         <  1:308062/139‑1 (MQ=255)
                     tCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTaaa                                                                                                         >  5:387885/1‑139 (MQ=255)
                            aGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTc                                                                                                   <  7:229993/138‑1 (MQ=255)
                              aaaGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAAt                                                                                                >  4:291946/1‑139 (MQ=255)
                                   tttAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAGAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttt                                                                                           <  8:126544/139‑1 (MQ=255)
                                   tttAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttt                                                                                           >  8:8977/1‑139 (MQ=255)
                                   tttAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttt                                                                                           >  8:411134/1‑139 (MQ=255)
                                        gATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGta                                                                                      <  1:76774/139‑1 (MQ=255)
                                            aGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                   <  6:333105/138‑1 (MQ=255)
                                            aGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                   <  1:2212/138‑1 (MQ=255)
                                             gATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtagta                                                                                 <  5:380950/139‑1 (MQ=255)
                                             gATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                   >  4:82374/1‑137 (MQ=255)
                                                   gAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCAtt                                                                            <  8:456593/138‑1 (MQ=255)
                                                   gAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCAtt                                                                            <  6:368912/138‑1 (MQ=255)
                                                    aaaaaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTgg                                                                          >  3:315735/1‑139 (MQ=255)
                                                     aaaaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATTCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTggg                                                                         <  3:412730/139‑1 (MQ=255)
                                                     aaaaTTTTATATATTTTAAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTggg                                                                         <  2:344418/139‑1 (MQ=255)
                                                      taaTTTTATTTATTTTCAATAAAATACTTGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCTTTGTAGTAAATTGTTCAATGTTTTTTGTATAGTATCATTGGGa                                                                        <  8:198914/138‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  4:315735/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  8:454607/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  3:368435/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  6:387885/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  6:4660/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  2:293191/139‑1 (MQ=255)
                                                        aTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATc                                                                      <  3:147647/139‑1 (MQ=255)
                                                             atatatTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAata                                                                  <  2:358213/138‑1 (MQ=255)
                                                              tatatTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatata                                                                >  2:217463/1‑139 (MQ=255)
                                                              tatatTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatat                                                                 <  8:215103/138‑1 (MQ=255)
                                                                   tttCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  5:125170/127‑1 (MQ=255)
                                                                   tttCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       >  6:125170/1‑127 (MQ=255)
                                                                      cAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTtat                                                         <  7:411134/138‑1 (MQ=255)
                                                                      cAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTtat                                                         <  5:148799/138‑1 (MQ=255)
                                                                      cAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTtat                                                         <  1:217463/138‑1 (MQ=255)
                                                                                 tGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTAc                                              >  8:339867/1‑138 (MQ=255)
                                                                                   gTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTaa                                           >  6:128084/1‑139 (MQ=255)
                                                                                   gTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTaa                                           >  1:332402/1‑139 (MQ=255)
                                                                                       tataATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCaa                                       <  2:261353/139‑1 (MQ=255)
                                                                                         taataaATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCaag                                      >  8:410334/1‑138 (MQ=255)
                                                                                            taaATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                  <  1:143640/139‑1 (MQ=255)
                                                                                            taaATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                  <  3:401718/139‑1 (MQ=255)
                                                                                            taaATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAAAAGt                                  <  5:355528/139‑1 (MQ=255)
                                                                                              aaTGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                  >  8:161565/1‑137 (MQ=255)
                                                                                                  tataCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGATCAATGTTTTTTGTATAGTAGCATTGTGATCAATATATTGTTATAAATTATTTACTAACCAAGAAGTCAGta                             >  8:287318/1‑138 (MQ=255)
                                                                                                   ataCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGTTTATAAATTATTTACTAAAATAGAAGTCATtag                            >  7:422127/1‑137 (MQ=255)
                                                                                                   ataCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGtac                            >  3:134784/1‑138 (MQ=255)
                                                                                                   ataCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACAAAGAAGTCAGtac                            >  7:278924/1‑138 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTTTAGTAGCATTGGGATCAAGAGATGGGTATAAGTTATTTACGACCGAAGACGGCAGtgc                            >  6:248606/1‑135 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTggg                                                                         <  8:333360/92‑1 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTggg                                                                         >  7:333360/1‑92 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       >  4:329090/1‑94 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  3:329090/94‑1 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGtacta                          <  6:400061/139‑1 (MQ=255)
                                                                                                     aCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGtactac                         <  2:382411/139‑1 (MQ=255)
                                                                                                      cAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTactact                        <  7:454399/139‑1 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGta                                                                                      <  3:159915/74‑1 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGta                                                                                      >  4:159915/1‑74 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTa                                            >  4:82778/1‑116 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTa                                            <  3:82778/116‑1 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGta                             <  5:1295/131‑1 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGt                     >  2:130168/1‑139 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGt                     >  4:375152/1‑139 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTAATAACCAAGAAGTCAGta                             >  6:1295/1‑131 (MQ=255)
                                                                                                          ttGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGtt                    <  3:144134/139‑1 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCATGAAGTCAGTACTACTCGtt                    >  7:282010/1‑138 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCg                      <  1:260734/136‑1 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCg                      >  2:260734/1‑136 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGtt                    <  8:404887/138‑1 (MQ=255)
                                                                                                            gTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  1:320590/86‑1 (MQ=255)
                                                                                                            gTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       >  2:320590/1‑86 (MQ=255)
                                                                                                            gTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGtt                    <  7:18807/137‑1 (MQ=255)
                                                                                                            gTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGtt                    >  8:18807/1‑137 (MQ=255)
                                                                                                             tATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                  >  6:51704/1‑122 (MQ=255)
                                                                                                             tATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                  <  5:51704/122‑1 (MQ=255)
                                                                                                                 gagaTGTAATTAATTAAAAAGGATTTTTTTAATTGCTTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTaataat             <  5:367654/139‑1 (MQ=255)
                                                                                                                  agaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTtaataa              >  6:274918/1‑137 (MQ=255)
                                                                                                                   gaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAg           >  3:157075/1‑139 (MQ=255)
                                                                                                                   gaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAg           >  1:300407/1‑139 (MQ=255)
                                                                                                                     tGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGc          <  3:291946/138‑1 (MQ=255)
                                                                                                                     tGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCt         >  6:391891/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCTATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTATCCAAGATGTCAGTACTACTCGTTTAATAATAGCTTc       >  2:336341/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       >  3:427367/1‑75 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                       <  4:427367/75‑1 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc       >  5:172845/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc       <  5:274918/139‑1 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc       >  8:42312/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc       >  5:395566/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAAGTATTTACTAACAAAGAAGTCAGTACTACTCTTTTAATAATAGCTTc       >  6:367654/1‑139 (MQ=255)
                                                                                                                            aattaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTaa  >  1:355910/1‑139 (MQ=255)
                                                                                                                              ttaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAGTATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTactact                        <  7:451950/115‑1 (MQ=255)
                                                                                                                              ttaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAGTATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTactact                        >  8:451950/1‑115 (MQ=255)
                                                                                                                                          |                                                                                                                             
TGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAA  >  BS168_reference_genome/557729‑557990

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: