Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,648,641 A→G 5.6% F4S (TTT→TCT)  yqeF ← putative lipoprotein; putative esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,648,6410AG 5.6% 217.0 / 4.2 71F4S (TTT→TCT) yqeFputative lipoprotein; putative esterase
Reads supporting (aligned to +/- strand):  ref base A (30/37);  new base G (2/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

ACCCACTGCTGTGCTTTTGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCTCAGTGAATATCATTC  >  BS168_reference_genome/2648503‑2648773
                                                                                                                                          |                                                                                                                                    
aCCCACTGCTGTGCTTTTGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATaa                                                                                                                                      <  7:79930/139‑1 (MQ=255)
     ctgctgTGCTTTTGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGt                                                                                                                                 <  1:252858/139‑1 (MQ=255)
          gtgCTTTTGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCAt                                                                                                                            >  6:312890/1‑139 (MQ=255)
                ttGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATaa                                                                                                                                      <  2:9600/123‑1 (MQ=255)
                ttGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATaa                                                                                                                                      >  1:9600/1‑123 (MQ=255)
                 tGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTac                                                                                                                      <  7:339551/138‑1 (MQ=255)
                        tCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCa                                                                                                               >  6:301909/1‑138 (MQ=255)
                              tAATCCCATTTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  2:102194/139‑1 (MQ=255)
                              tAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  5:301909/139‑1 (MQ=255)
                              tAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  2:294736/139‑1 (MQ=255)
                              tAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  2:1666/139‑1 (MQ=255)
                                aTCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTc                                                                                                       <  1:351684/138‑1 (MQ=255)
                                   ccAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACtt                                                                                                   <  5:312890/139‑1 (MQ=255)
                                   ccAATTAGAGACACGCGTATTACTGTCACCAAATGCTACAACGGCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACtt                                                                                                   <  2:113748/139‑1 (MQ=255)
                                    cAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACttt                                                                                                  <  3:314029/139‑1 (MQ=255)
                                    cAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACttt                                                                                                  <  3:156961/139‑1 (MQ=255)
                                      aatAGAGCCACGCGAATTACAGTCACCAAATGCTACAACGTCTTCGTATCCACTTCCCTCGCATCATGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGATTCATCTTTACACCCCCAAAAGCCTCACTGTaa                                                                                                >  3:59541/3‑139 (MQ=255)
                                      aTTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTaa                                                                                                <  8:46842/139‑1 (MQ=255)
                                      aTTAGAGCCAAGAGTATTACTTTCAACAAATTCTACAACGTCTTCGTATCCACTTCCCTAGCATCCTGCTGTAACTAACAGCATAAGGATAATTATTATAAAATTTTTCATCTTTACACCCCCAAAAGCCTCACTTTaa                                                                                                <  8:49043/139‑1 (MQ=255)
                                         agagCCACGAGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATg                                                                                              <  4:262360/138‑1 (MQ=255)
                                            gCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTtat                                                                                          <  8:143744/139‑1 (MQ=255)
                                              cACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATAc                                                                                        >  7:95201/1‑139 (MQ=255)
                                               aCGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATAc                                                                                        >  8:117590/1‑138 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATaaa                                                                                                                                     >  1:60381/1‑92 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATaaa                                                                                                                                     <  2:60381/92‑1 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCa                                                                                                                             <  8:256089/100‑1 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCa                                                                                                                             >  7:256089/1‑100 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        >  4:71830/1‑121 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        >  5:46996/1‑121 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  3:71830/121‑1 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  6:46996/121‑1 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACct                                                                                      <  8:95201/139‑1 (MQ=255)
                                                cgcgTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACct                                                                                      <  7:117590/139‑1 (MQ=255)
                                                          tGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTctct                                                                            >  1:277001/1‑139 (MQ=255)
                                                            tCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTc                                                                          >  4:340608/1‑139 (MQ=255)
                                                            tCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTc                                                                          >  3:350398/1‑139 (MQ=255)
                                                             cACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCg                                                                         >  8:312532/1‑139 (MQ=255)
                                                              aCCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        >  1:116562/1‑107 (MQ=255)
                                                              aCCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  2:116562/107‑1 (MQ=255)
                                                              aCCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTaa                                                                                                <  4:263727/115‑1 (MQ=255)
                                                              aCCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTaa                                                                                                >  3:263727/1‑115 (MQ=255)
                                                              aCCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGt                                                                        >  4:65600/1‑139 (MQ=255)
                                                              aCCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGt                                                                        >  8:101085/1‑139 (MQ=255)
                                                               ccAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTaca                                                                                                                     <  1:161266/93‑1 (MQ=255)
                                                               ccAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTaca                                                                                                                     >  2:161266/1‑93 (MQ=255)
                                                               ccAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGtt                                                                       <  2:277001/139‑1 (MQ=255)
                                                                      cTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGaa                                                                >  8:21750/1‑139 (MQ=255)
                                                                      cTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGaa                                                                <  7:312532/139‑1 (MQ=255)
                                                                          accGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAACAGCCTCACTTTAATCTTATACATCTGTTTCTCTGCGTTCAGTGAAGAGt                                                            >  8:274803/3‑139 (MQ=255)
                                                                              tCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGGTATACCTCTGTTTCTCTTCGTTCAATGAAGAGttttt                                                        >  4:245590/1‑139 (MQ=255)
                                                                                  cGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAATTGATTATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCg                                                      <  4:59541/137‑1 (MQ=255)
                                                                                      tCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCa                                                <  1:69245/139‑1 (MQ=255)
                                                                                      tCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCa                                                <  3:65600/139‑1 (MQ=255)
                                                                                         ccTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCa                                                <  4:294114/136‑1 (MQ=255)
                                                                                         ccTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCa                                                >  3:294114/1‑136 (MQ=255)
                                                                                            tCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        >  5:77947/1‑77 (MQ=255)
                                                                                            tCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCt                                                                                                        <  6:77947/77‑1 (MQ=255)
                                                                                              ccTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAGAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGttttt                                        >  6:73493/1‑139 (MQ=255)
                                                                                              ccTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAGAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGttttt                                        >  6:68129/1‑139 (MQ=255)
                                                                                                  gCATCCTGCTGTATCGAACATCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTcc                                     <  7:274803/138‑1 (MQ=255)
                                                                                                  gCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAGAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTcc                                     <  5:68129/138‑1 (MQ=255)
                                                                                                  gCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAGAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTcc                                     <  5:73493/138‑1 (MQ=255)
                                                                                                  gCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTcc                                     <  3:324763/138‑1 (MQ=255)
                                                                                                          ctgTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTGTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGATATCACTGTTTTTTCCTATAGTTaa                            >  4:68982/1‑139 (MQ=255)
                                                                                                             tAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTATACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAAc                           <  1:287249/137‑1 (MQ=255)
                                                                                                              aaCGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACgtg                        >  4:216838/1‑139 (MQ=255)
                                                                                                                  aaCAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGATGAGTTTTTCGATATCACTGTTTTTTCCTATAGTTAACGTGTccc                    >  8:265884/1‑139 (MQ=255)
                                                                                                                  aaCAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTccc                    >  6:281595/1‑139 (MQ=255)
                                                                                                                     agcagAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCt                 <  3:340608/139‑1 (MQ=255)
                                                                                                                     agcagAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCt                 <  7:101085/139‑1 (MQ=255)
                                                                                                                       cagAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGttt                                                          <  3:251011/96‑1 (MQ=255)
                                                                                                                       cagAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGttt                                                          >  4:251011/1‑96 (MQ=255)
                                                                                                                       cagAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCt                 >  1:187233/1‑137 (MQ=255)
                                                                                                                             gAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCTCAGTGaa          >  2:362757/1‑138 (MQ=255)
                                                                                                                               aaaGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCATATAGTTAACGTGTCACCCTCAGTGAatat       >  5:300848/1‑139 (MQ=255)
                                                                                                                                    atgatAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCTCAGTGAATATCATTc  <  5:281595/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
ACCCACTGCTGTGCTTTTGGATAATCGCGATAATCCCAATTAGAGCCACGCGTATTACTGTCACCAAATGCTACAACGTCTTCGTATCCCCTTCCCTCGCATCCTGCTGTAACGAACAGCAGAAGGAAAAGGATGATAAAATGTTTCATCTTTACACCCCCAAAAGCCTCACTTTAATGTTATACCTCTGTTTCTCTTCGTTCAATGAAGAGTTTTTCGCTATCACTGTTTTTTCCTATAGTTAACGTGTCCCCCTCAGTGAATATCATTC  >  BS168_reference_genome/2648503‑2648773

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