Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 165,751 Δ1 bp intergenic (+44/‑3) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome165,7510C.97.8% 234.1 / 4.5 92intergenic (+44/‑3)rrnI‑5S/trnI‑Asnribosomal RNA‑5S/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (35/54);  minor base T (0/2);  total (35/57)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATCTCATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGC  >  BS168_reference_genome/165693‑165889
                                                          |                                                                                                                                           
ggACGCCGGCAAGCAGATTTATTAATTCTATCATGTTTGTTTATGAAGGTTTATATTATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTGACCGATCGGACGTAGGTTCGAATCCTAGCTGCGGAGCCATTGCTTCCATAGcaca                                                <  1:3350701/150‑1 (MQ=18)
     ccgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAgg                                           <  2:1589837/150‑1 (MQ=21)
     ccgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAgg                                           >  2:3645031/1‑150 (MQ=21)
     ccgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAgg                                           <  1:3220390/150‑1 (MQ=21)
     ccgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAgg                                           <  2:867123/150‑1 (MQ=21)
     ccgccAAGCAGATTTATTAATTCTACCATGTATGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAgg                                           <  2:1709992/150‑1 (MQ=21)
      cgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAgga                                          >  2:2004387/1‑149 (MQ=21)
      cgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGt                                          >  2:3884995/1‑150 (MQ=21)
      cgccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGt                                          <  1:3420709/150‑1 (MQ=21)
       gccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTa                                         >  1:3612940/1‑150 (MQ=21)
       gccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTa                                         <  2:260530/150‑1 (MQ=21)
        ccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTag                                        <  2:2345110/150‑1 (MQ=21)
         cAAGCAGATTTTTTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTCCCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTaga                                       <  2:1063510/150‑1 (MQ=21)
         cAAGCAGATTTATTAATCCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTaga                                       <  1:1927361/150‑1 (MQ=21)
          aaGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTagag                                      <  1:1398927/150‑1 (MQ=21)
           aGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGc                                     >  2:3311605/1‑150 (MQ=21)
            gCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTGTATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCa                                    <  2:2981565/150‑1 (MQ=21)
            gCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCa                                    <  1:3597302/150‑1 (MQ=21)
            gCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCa                                    <  2:3882219/150‑1 (MQ=21)
             cAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCAc                                   <  2:2397008/150‑1 (MQ=21)
               gATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCCCtt                                 >  2:3040631/1‑150 (MQ=21)
               gATTTATTAATTCTACCATGTTTGTTTATGAAGCTT‑‑TATCTTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACtt                                 <  2:2974160/150‑1 (MQ=21)
                aTTTATTACTTCTACCATGTCTGTTTATGAAGCTTTATATTACTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTACCCGATCGGTCGAAGGTTCGCATCCCACCTGCGGAGCCCTTGCTCCCATAGCTCAGCAGGTAGAGCACTTc                                <  2:4033190/150‑1 (MQ=18)
                aTTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTc                                <  1:2892004/150‑1 (MQ=21)
                aTTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTc                                >  2:734793/1‑150 (MQ=21)
                aTTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTc                                >  2:52319/1‑150 (MQ=21)
                aTTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTc                                <  2:800981/150‑1 (MQ=21)
                 tttATTAATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTcc                               >  2:2628607/1‑150 (MQ=21)
                  ttattaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCa                              <  1:3408341/150‑1 (MQ=21)
                  ttattaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCa                              >  1:1302857/1‑150 (MQ=21)
                   tattaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTACCCGATCGGCCGTAGGTTCGAATCCTGCCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCAt                             <  2:10056/150‑1 (MQ=21)
                   tattaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTGGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCAt                             <  1:2641758/150‑1 (MQ=21)
                   tattaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCAt                             >  1:937201/1‑150 (MQ=21)
                   tattaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCAt                             >  2:3614990/1‑150 (MQ=21)
                    attaATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGGGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATg                            >  1:1704390/1‑150 (MQ=21)
                      taATTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGt                          >  1:2943756/1‑150 (MQ=21)
                       aaTTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTa                         >  1:1965252/1‑150 (MQ=21)
                       aaTTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTa                         >  2:2411630/1‑150 (MQ=21)
                       aaTTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCACAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTa                         <  1:2100337/150‑1 (MQ=21)
                        aTTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTaa                        <  2:3964979/150‑1 (MQ=21)
                        aTTCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGAAGGTTCGAATCCTACCTGCGGAGCCATTGCTACCATAGCTCAGCAGGTAGAGCACTTCCTAGGTaa                        >  1:2849838/1‑150 (MQ=21)
                         ttCTACCATTTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAg                       <  1:3618025/150‑1 (MQ=21)
                         ttCTACCATGTTTGTTTATGAAGCTTTATATTCTTCCGCAGTACCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAg                       <  1:3960128/150‑1 (MQ=21)
                         ttCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAg                       <  2:2269947/150‑1 (MQ=21)
                          tCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAgg                      >  2:1674043/1‑150 (MQ=21)
                          tCTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCAAAGCTCACCAGGGAGGGCACTTCCATGGTAAgg                      >  1:2998180/1‑150 (MQ=18)
                           cTACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGa                     <  1:529182/150‑1 (MQ=21)
                           cTACCATGTGTGTTTATGAAGCTTTATATTATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGCCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGa                     <  1:2120702/150‑1 (MQ=21)
                            tACCGTGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTTACCGATCGGTGGTAGGTTCGAATCCTCCCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAggta                    <  2:2657071/150‑3 (MQ=21)
                            tACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGaa                    <  2:3381869/150‑1 (MQ=21)
                            tACCATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGaa                    <  2:3331931/150‑1 (MQ=21)
                              ccATGTTTGTTTATGAAGCTTTATATTGTTCCGCAGTAGCTCAGTGGTAGAGCTATGGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAaga                  <  1:3821161/150‑1 (MQ=21)
                              ccATGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCACTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAaga                  >  1:3155021/1‑150 (MQ=21)
                                aTGTTTGTTTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAgg                <  2:1350847/150‑1 (MQ=21)
                                aTGTTTGTTTATGAAGCTTAATATTATTCCGCAGTAGCACAGTGGTAGAGCTATGGGCTGTTAACCGATGGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAgg                <  1:1790261/150‑1 (MQ=21)
                                 tgtttgttTATGAAGCTTTATATTATTCCGCATTAGCCCAGTGGTAGAGTTATTGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGt               <  1:1283638/150‑1 (MQ=21)
                                 tgtttgttTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGt               <  2:2224317/150‑1 (MQ=21)
                                 tgtttgttTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGt               <  1:201611/150‑1 (MQ=21)
                                 tgtttgttTATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGACGTAGGTTCGAGTCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGt               >  1:3967546/1‑150 (MQ=21)
                                  gtttgtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGGGCCATTGCTTCCATAGCTCAGCAGGGAGAGCACTTCCATGGTAAGGAAGAGGTc              >  1:2689404/1‑150 (MQ=21)
                                  gtttgtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTc              <  2:3746638/150‑1 (MQ=21)
                                  gtttgtttAGGAAGCTTTATATTATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATGCTACCTGGGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACGTCCATGGTAAGGAAGAGGTc              <  1:2196140/150‑1 (MQ=21)
                                   tttgtttATGAAGCTTTATATTATCCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCa             <  2:668056/150‑1 (MQ=21)
                                    ttgtttgtGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAg            <  1:2628607/150‑1 (MQ=21)
                                    ttgtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAg            <  2:216719/150‑1 (MQ=21)
                                    ttgtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGAAAGGAAGAGGTCAg            >  2:1575602/1‑150 (MQ=21)
                                    ttgtttATGAAGCTTCATATTATTCCGCAGAAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAg            <  2:877127/150‑1 (MQ=21)
                                    ttgtttATGAAGCTT‑‑TATCTTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTACCCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAg            <  1:3822204/150‑1 (MQ=21)
                                      gtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGGGGTCAGCg          >  2:1315505/1‑150 (MQ=21)
                                      gtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg          <  1:123391/150‑1 (MQ=21)
                                      gtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg          <  1:879350/150‑1 (MQ=21)
                                      gtttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGCAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg          >  1:2484009/1‑150 (MQ=21)
                                      gtttATGAAGCTTTATATTATCCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg          <  1:3724039/150‑1 (MQ=21)
                                       tttATGAAGCTTTATATTATTCCGCAGTAGCTCGCTGGTAGAGCTATCGGCTGTTACCGGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCgg         <  2:2191075/150‑1 (MQ=21)
                                        ttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGt        >  1:3312829/1‑150 (MQ=21)
                                        ttATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCGGTTACCCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGt        <  2:2500121/150‑1 (MQ=21)
                                         tATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGtt       <  1:2677742/150‑1 (MQ=21)
                                         tATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGtt       >  1:2163278/1‑150 (MQ=21)
                                         tATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGtt       >  2:1268731/1‑150 (MQ=21)
                                         tATGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGtt       >  2:1799155/1‑150 (MQ=21)
                                          aTGAAGCTTTATATTATTCCGCAGTAGCTCGGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc      <  2:3179211/150‑1 (MQ=21)
                                          aTGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCCTAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc      <  2:173256/150‑1 (MQ=21)
                                          aTGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc      <  2:2121966/150‑1 (MQ=21)
                                          aTGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCAGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc      <  1:2735651/150‑1 (MQ=21)
                                           tGAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCg     <  2:1649073/150‑1 (MQ=21)
                                            gAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGGGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGa    >  1:852458/1‑150 (MQ=21)
                                            gAAGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGa    <  1:3072927/150‑1 (MQ=21)
                                             aaGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCTTAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAg   >  1:1100360/1‑150 (MQ=21)
                                             aaGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAg   >  2:2389028/1‑150 (MQ=21)
                                             aaGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAg   >  2:34152/1‑150 (MQ=21)
                                              aGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGGGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGc  >  1:71426/1‑150 (MQ=21)
                                              aGCTTTATATTATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGc  >  2:3791791/1‑150 (MQ=21)
                                                          |                                                                                                                                           
GGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATCTCATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTGCGGAGCATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGC  >  BS168_reference_genome/165693‑165889

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: