Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,020,611 A→C 16.5% S99S (TCT→TCG yxzG ← putative nucleic acid binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,020,6110AC16.5% 136.5 / 6.5 67S99S (TCT→TCGyxzGputative nucleic acid binding protein
Reads supporting (aligned to +/- strand):  ref base A (41/14);  new base C (2/10);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 3.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGAAGAGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACACCCGTCTCACGACCGTTTCATC  >  BS168_reference_genome/4020479‑4020748
                                                                                                                                    |                                                                                                                                         
acCAGAAGAGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAgacgac                                                                                                                                     >  2:28260/1‑139 (MQ=255)
 ccAGAAGAGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCGTGACTCACATCGGAAGAAATCAAATAAAGAGACGACt                                                                                                                                    >  1:64909/1‑139 (MQ=255)
  ctgacgagGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACAGCTAGGCCTTCCCCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGCCAGAAATCAAATAAAGAGACGACtt                                                                                                                                   <  2:102082/134‑1 (MQ=255)
  cAGAAGAGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACt                                                                                                                                    >  8:245457/1‑138 (MQ=255)
  cAGAAGAGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCACAAATAAAATAAAGCCACGACtt                                                                                                                                   <  2:289439/139‑1 (MQ=255)
  cAGAAGAGCCCACTGCGTAACCCGCGAAACTTGTCAGTTTCAGCTGTTCAACATCTCGGCCTCCACCGTTCGAAGCCCAGCATGTTATCTCATCCATCTCTTGACCCACCTCGGCAGAAATCAAAAAACGCCACGACtt                                                                                                                                   <  2:253338/139‑1 (MQ=255)
     aagaGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACTGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGCCGACTTAAg                                                                                                                                >  1:168221/1‑139 (MQ=255)
       gagGCCACTGAGTAAACAGCTTAAATTGTCATTTTCAGCTGTTTAACATCTTGGCATTCACCGTTCCAATCACAGCATTTTATCTCTTCCATCTCTTGACTAACATCGGAATAAAAAAAATAAACACACGACTTAAGGt                                                                                                                              <  6:333255/139‑1 (MQ=255)
          gCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTCATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTACTGTCtc                                                                                                                           >  5:129011/1‑138 (MQ=255)
          gCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGGTATCTCTTCCCTCTCTTGACTCACATCGGCAGAAATCAAATAAAGCGACGACTTAAGGTCTg                                                                                                                           >  4:98918/1‑139 (MQ=255)
           ccACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATAAAATAACGACACGACTTAAGGTCTGc                                                                                                                          <  1:166443/139‑1 (MQ=255)
            cACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGCCCTCTGAAGGTGTGCa                                                                                                                         >  7:303968/1‑139 (MQ=255)
            cACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCa                                                                                                                         >  8:315892/1‑139 (MQ=255)
             acggCGAAAACAGATTAAATGGGCAGTTTCAGCTGATGAAAATATAGCCCTTCACGGTTCGAATCACAGCTTGTTAACTCATCCAGCTCTGCACTCCACTCGGCAGAAACAAAATAACGCCACGACTTAAGGTCTGCAc                                                                                                                        <  4:17768/136‑1 (MQ=255)
               tGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGCCTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTg                                                                                                                      >  4:346577/1‑139 (MQ=255)
               tGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTg                                                                                                                      >  5:333255/1‑139 (MQ=255)
                    aaaCAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCTAATAAAGAGACGACTTAAGGTCTGCCCTGTAAc                                                                                                                  >  2:186393/1‑138 (MQ=255)
                    aaaCAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAAc                                                                                                                  >  8:3844/1‑138 (MQ=255)
                         gCGTAACTTGTCTGTTTAAGCTGTTCAACATCTAGGCATTCACCTTTCGAATCAAAGATTTTTATCTCTACCATTTCTTCACTCAAATCGGCACAAATCAAATAAAGACACGACTTAAGGTCTGCACTGTAACAATTGa                                                                                                            <  8:9963/139‑1 (MQ=255)
                         gCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACCGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTg                                                                                                             >  1:195693/1‑138 (MQ=255)
                         gCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTg                                                                                                             >  4:227436/1‑138 (MQ=255)
                           gTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTGTGCACTGTAACAATTGCAc                                                                                                          >  7:208248/1‑139 (MQ=255)
                           gTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAAc                                                                                                          >  3:429422/1‑139 (MQ=255)
                           gTAACATGTCAGATTCAGCAGTTGAACATCTAGGCCTTCACCGTACGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGCCAGAATTCAAATAAAGAGACGACTTAAGCTCTGCACAGTAACAATTCAAc                                                                                                          >  6:117441/1‑139 (MQ=255)
                            tAACTTGTCAGTTTCAGCTGTTGCACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGCCTCACATCGGCAGCAATCAAATAAAGCGACGACTTAAGGTCTGCACTGTAACAAGTGCACg                                                                                                         >  6:51699/1‑139 (MQ=255)
                                  gTCAGTTTGAGCTGTTGAACAGCTAGTCCTGCACCGGGCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGTCTTAAGGTCTGCACTGAAACAATTGAACGTCcca                                                                                                    <  8:303968/138‑1 (MQ=255)
                                    cAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCCCATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCccaac                                                                                                  >  3:350868/1‑138 (MQ=255)
                                    cAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCccaacc                                                                                                 >  7:9963/1‑139 (MQ=255)
                                    cAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCccaac                                                                                                  >  7:110273/1‑138 (MQ=255)
                                                 tGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCGTCCTTCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCAAACCAAAACACcatctt                                                                                    >  5:355769/1‑139 (MQ=255)
                                                  gAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATACAGAGGCGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACcttctt                                                                                    <  7:306654/138‑1 (MQ=255)
                                                    ccATCTAGCCCTTCACGGGTCGAATCCCACCTTTATATCGCTACCATCTCTTCACTCACCTCCGCAGTAAACAAATAAACCCACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCgg                                                                                 <  4:2434/138‑1 (MQ=255)
                                                    aCATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAAAAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCgg                                                                                 <  5:156806/139‑1 (MQ=255)
                                                          aGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATcg                                                                           >  4:54083/1‑139 (MQ=255)
                                                          aGGCCTTCACCGTTCAAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATcg                                                                           >  6:107122/1‑139 (MQ=255)
                                                              cTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAg                                                                       <  7:245457/139‑1 (MQ=255)
                                                              cTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATAAAATAACGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAg                                                                       <  3:411050/139‑1 (MQ=255)
                                                                      ttCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAAc                                                                                          >  1:6111/1‑112 (MQ=255)
                                                                      ttCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAAc                                                                                          <  2:6111/112‑1 (MQ=255)
                                                                      tgcGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATAAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCaa                                                               <  8:77858/137‑1 (MQ=255)
                                                                       cccaaacacaCCATGTAATCTATCCCATCTCTTCACACACATCGCAAAAAAAAAAAAAAACCCACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCaa                                                               <  3:227436/132‑1 (MQ=255)
                                                                              acaGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGCACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCt                                                       >  6:679/1‑139 (MQ=255)
                                                                              acaGCATGTTATCTCATCCATCTATTGACACACATCGGCATAAAACAAATAAAGCCACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCt                                                       <  3:123368/139‑1 (MQ=255)
                                                                                aGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCt                                                       >  8:293856/1‑137 (MQ=255)
                                                                                    tGTTATCTCTTCCATCTCTTGACTCCCATCGGCAGAAATCAAAAAAAGAGACGCCTTAAGGTCTGCACTGTAACAAATGAACGTCCCATCCAAAACACCTACTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTcc                                                 >  1:60459/1‑139 (MQ=255)
                                                                                       tATCTCTTCCATCTCATGACTCACATCGGAAGAAATCAAATAAAGAGACGACTTAAGGGCTGCACTGTAACAATTGAACTTCCCATCCAAAAACCCTTCTTCTGATATCCCTAGAAGAACAATTCCTTCTGAGTCCGGc                                              >  4:2963/1‑139 (MQ=255)
                                                                                             ttCCACCTCTTCACACACATCGGCACAAATAAAAAACCGCGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGctct                                           <  1:321125/136‑1 (MQ=255)
                                                                                              tCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGGTATCGCGAGCTGATCTAATGCTTCTGAGTCCGGCGCTTTAt                                       >  2:386334/1‑139 (MQ=255)
                                                                                              tCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTAt                                       >  7:180311/1‑139 (MQ=255)
                                                                                              tCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTAt                                       >  4:115189/1‑139 (MQ=255)
                                                                                              tCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTATGCACTGTAACAATTGAACGTCACAACCAAAAAAACTTCTTCTGATATCGCGAGCTGATCAAATGATTCTGATTACGGCTCTTTa                                        >  7:188921/1‑138 (MQ=255)
                                                                                               ccATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATc                                      <  4:429422/139‑1 (MQ=255)
                                                                                                  tctcTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGa                                   >  3:154141/1‑139 (MQ=255)
                                                                                                     ctTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCGGTATCCGa                                   >  6:194032/1‑136 (MQ=255)
                                                                                                      ttGACTCACATCGGCAGAAATCAAATAAAGCGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACa                                 <  2:14667/137‑1 (MQ=255)
                                                                                                       tGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGCTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGa                               >  3:396423/1‑138 (MQ=255)
                                                                                                           taacaTCGGCAGAAATCAAATAAAGAGACGACTTAATGTCTGCACTGTAAAAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTTATCAAATGCTTCTGAGTCCGGCTCTTTATACGAAAGAATCca                          >  4:178047/3‑139 (MQ=255)
                                                                                                             acaTCGCCAGAAATCAAATAAAGACACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCcaca                        <  5:107122/139‑1 (MQ=255)
                                                                                                              caTCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCacac                       >  8:389349/1‑139 (MQ=255)
                                                                                                              caTCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGGAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCcaccc                       >  4:91581/1‑137 (MQ=255)
                                                                                                               aTCGGCAGAAATCAAATAAACACACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACAcc                      <  3:178047/139‑1 (MQ=255)
                                                                                                               aTCGCCACAAATAAAAAAACGCGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACAcc                      <  1:386334/139‑1 (MQ=255)
                                                                                                                 cgccAGAATTAAAATAAAGCCACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACAcc                      <  4:154141/134‑1 (MQ=255)
                                                                                                                       aaaTCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAGCCAACAAACCTTCTTCGGATTTCGCGAGCTGAGCAAATGCTACTGAGTCCGGCTCTATAACCGACAGAATCCAGACCCGtctccc              >  2:202100/1‑137 (MQ=255)
                                                                                                                         aTCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACACCCGTCTCACGa            >  4:126921/1‑139 (MQ=255)
                                                                                                                         aTCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACACCCGTCTCACGa            >  4:12737/1‑139 (MQ=255)
                                                                                                                             aaTAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCACACCCAAAAAACCTTCTTCGGATATCGCGAGCAGATCAAATGCTTCTGAGTCCGGCTATTTATCCGACAGCATCCACACCCGTCTCACGACCgg        >  4:351548/1‑138 (MQ=255)
                                                                                                                                    agaCGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACACCCGTCTCACGACCGTTtcatc  >  3:100332/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                         
ACCAGAAGAGGCCACTGCGTAAACAGCGTAACTTGTCAGTTTCAGCTGTTGAACATCTAGGCCTTCACCGTTCGAATCACAGCATGTTATCTCTTCCATCTCTTGACTCACATCGGCAGAAATCAAATAAAGAGACGACTTAAGGTCTGCACTGTAACAATTGAACGTCCCAACCAAAAAACCTTCTTCGGATATCGCGAGCTGATCAAATGCTTCTGAGTCCGGCTCTTTATCCGACAGAATCCACACCCGTCTCACGACCGTTTCATC  >  BS168_reference_genome/4020479‑4020748

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