Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,935,823 Δ1 bp 100% intergenic (+35/+1) ywbE → / ← rlmI conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,935,8230T.95.4% 224.6 / 4.7 65intergenic (+35/+1)ywbE/rlmIconserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base . (38/24);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935685‑3935953
                                                                                                                                          |                                                                                                                                  
cAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCGTGTTCat                                                                                                                                    >  5:267072/1‑139 (MQ=255)
  gTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCatta                                                                                                                                  >  2:298687/1‑138 (MQ=255)
     aaGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCattattt                                                                                                                               >  3:119489/1‑135 (MQ=255)
        gATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAtttttt                                                                                                                           >  5:45412/1‑139 (MQ=255)
        gATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAttttt                                                                                                                            >  3:318574/1‑138 (MQ=255)
           atTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCCGCTTTTTTATTATTTTTTCTg                                                                                                                        >  6:247100/1‑139 (MQ=255)
             tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGCCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTg                                                                                                                        >  3:87144/1‑137 (MQ=255)
             tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTg                                                                                                                        >  4:105929/1‑137 (MQ=255)
                    aaaaCAAACTTTAATCCGCATGGAATAAAAGACCGTCTTCAAGAGGCCCCCATCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                               <  3:424303/139‑1 (MQ=255)
                               tCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAgg                                                                                                     >  8:30727/1‑138 (MQ=255)
                                cATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAgg                                                                                                     >  6:312292/1‑137 (MQ=255)
                                      cATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGcc                                                                                             >  7:8773/1‑139 (MQ=255)
                                       aTGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGCTCGGTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGGTTTTTCATTATTTTTTATGCAGAAGGAACACTTTCAGGTAGTTgacg                                                                                            >  1:320719/1‑136 (MQ=255)
                                       aTGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAACAGCCCCCGGACCGGGCCAGGTTTTGATTTTTTTTTGTTGATGAGGGACCAGTTTCGGTTGGTTGctt                                                                                            >  1:197362/1‑137 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGCCCACTTTCAGGTATTTGCCTTCCGGAAATGCCg                                                                               >  7:352873/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  1:62980/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  1:62070/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  6:390466/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGACTTCCGGAAATGCCg                                                                               >  5:184056/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGa                                                                              >  4:207428/1‑138 (MQ=255)
                                                      gTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGa                                                                              >  5:188304/1‑138 (MQ=255)
                                                        cTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAAt                                                                           >  3:307284/1‑139 (MQ=255)
                                                                gggCCGCAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   <  3:145461/138‑1 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGCTCGTTTCTACATAAAAAAAAGCTGCCCCTGCACAGGGCAGCGTTTTGATTATTTTTTCTGCAGAAGGACCACGTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  3:290324/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTTCCCCTGAACAGGGCAGCTTTTTCATAATTTTTTATTCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  5:234614/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  6:151293/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  7:413515/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   <  8:8773/138‑1 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACACGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  1:163007/1‑138 (MQ=255)
                                                                      cATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcct                                                             >  6:422243/1‑139 (MQ=255)
                                                                         cGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCg                                                          >  6:217486/1‑139 (MQ=255)
                                                                         cGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCg                                                          >  2:400519/1‑139 (MQ=255)
                                                                          ggCCGCGACCAGCAGATCTTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCgg                                                         <  7:30727/139‑1 (MQ=255)
                                                                                      cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTTTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             >  6:166616/1‑139 (MQ=255)
                                                                                      cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGt                                                                                                 >  8:134054/1‑87 (MQ=255)
                                                                                      cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGt                                                                                                 <  7:134054/87‑1 (MQ=255)
                                                                                      cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             >  3:421324/1‑139 (MQ=255)
                                                                                      cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGGAAATTc                                             >  5:25995/1‑139 (MQ=255)
                                                                                            gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  5:138787/138‑1 (MQ=255)
                                                                                            gTTTATACATAAAAAAATGCTGCCCCTGAACAGGGAAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  4:290324/138‑1 (MQ=255)
                                                                                               tCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGa                                     <  4:307284/138‑1 (MQ=255)
                                                                                                    cataaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                <  3:207428/135‑1 (MQ=255)
                                                                                                     tAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTTATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                              <  4:242897/139‑1 (MQ=255)
                                                                                                       aaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGc                            >  4:155276/1‑139 (MQ=255)
                                                                                                        aaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCg                           <  8:98835/139‑1 (MQ=255)
                                                                                                             aGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCg                           >  8:113134/1‑134 (MQ=255)
                                                                                                             aGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCg                           <  7:113134/134‑1 (MQ=255)
                                                                                                                  ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  <  6:188304/138‑1 (MQ=255)
                                                                                                                  ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  <  3:105929/138‑1 (MQ=255)
                                                                                                                    ccTGAACAGGGCAGCTTTTTCATTTTTTTTTCTGCAGAAGGACCACTTTCAGGGAGTTGCCTTACGGAAAAGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCCCAAAGATCGTGTAGCGGTCATTTGTTTc               >  2:378863/1‑139 (MQ=255)
                                                                                                                            gggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAgg       >  7:22920/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  4:318435/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  4:427245/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  5:169740/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  1:311354/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  7:390130/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  5:376708/1‑139 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      <  2:163007/139‑1 (MQ=255)
                                                                                                                             ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACAACTTTCAGGTAGTTGCCTTCCGGAAATGCCTAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTAATCAATTATCGTGTAGAGTTCATTTGTTTATTTAAAtgc      >  6:411013/1‑139 (MQ=255)
                                                                                                                                aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  1:400519/139‑1 (MQ=255)
                                                                                                                                aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  5:390466/139‑1 (MQ=255)
                                                                                                                                aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  5:217486/139‑1 (MQ=255)
                                                                                                                                aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  1:413316/139‑1 (MQ=255)
                                                                                                                                 gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCTGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  6:169740/139‑1 (MQ=255)
                                                                                                                                 gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  6:376708/139‑1 (MQ=255)
                                                                                                                                 gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  1:309614/139‑1 (MQ=255)
                                                                                                                                 gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  8:390130/139‑1 (MQ=255)
                                                                                                                                 gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  8:413515/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
CAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935685‑3935953

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: