Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 340,408 (GC)6→5 coding (682‑683/948 nt) mhpR ← DNA‑binding transcriptional activator, 3HPP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606340,3980G.98.3% 553.6 / 1.2 115coding (693/948 nt)mhpRDNA‑binding transcriptional activator, 3HPP‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (62/51);  total (63/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*REL606340,3990C.97.4% 543.6 / 3.4 115coding (692/948 nt)mhpRDNA‑binding transcriptional activator, 3HPP‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (62/50);  total (63/52)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGG  >  REL606/340323‑340467
                                                                           ||                                                                     
gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTg                                                                        <  2:4179160/76‑1 (MQ=255)
gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTg                                                                        <  1:1265246/76‑1 (MQ=255)
gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTg                                                                        <  2:4764806/76‑1 (MQ=255)
gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTg                                                                        <  2:196550/76‑1 (MQ=255)
gcgcAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTg                                                                        <  1:3024103/76‑1 (MQ=255)
    aGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcg                                                                    >  1:691069/1‑76 (MQ=255)
    aGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcg                                                                    >  2:5480707/1‑76 (MQ=255)
    aGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcg                                                                    >  1:3233132/1‑76 (MQ=255)
    aGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTccgcg                                                                    >  1:3391873/1‑76 (MQ=255)
       ggTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgc                                                                 >  1:4641428/1‑76 (MQ=255)
       ggTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgc                                                                 <  2:1959274/76‑1 (MQ=255)
          aCGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgcc                                                              >  1:117242/1‑75 (MQ=255)
          aCGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgcc                                                              >  1:2529494/1‑75 (MQ=255)
           cGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgcca                                                             <  1:705646/76‑3 (MQ=255)
           cGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgcca                                                             <  2:6139081/76‑3 (MQ=255)
            ggCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgccag                                                            >  2:584617/1‑73 (MQ=255)
            ggCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgccag                                                            <  2:234809/76‑4 (MQ=255)
            ggCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTgcgcgcgcgccag                                                            >  1:783744/1‑73 (MQ=255)
                 tAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATgg                                                     >  1:4052789/1‑76 (MQ=255)
                  agagGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGc                                                    >  2:2817119/1‑76 (MQ=255)
                  agagGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGc                                                    >  1:950932/1‑76 (MQ=255)
                    agGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCtt                                                  <  2:1648191/76‑1 (MQ=255)
                     ggCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTc                                                 <  1:5026026/76‑1 (MQ=255)
                      gCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTcc                                                >  1:4791937/1‑76 (MQ=255)
                        gATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTcccg                                              >  1:112392/1‑74 (MQ=255)
                        gATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAg                                              >  2:3799929/1‑76 (MQ=255)
                         aTCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGc                                             >  2:687890/1‑76 (MQ=255)
                         aTCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGc                                             >  2:4072306/1‑76 (MQ=255)
                           cttctCCTCTTGACCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCtt                                           <  2:2091898/76‑1 (MQ=255)
                           cttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTTCAGCtt                                           >  1:837982/1‑76 (MQ=255)
                           cttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCtt                                           <  1:3863407/76‑1 (MQ=255)
                           cttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCtt                                           <  2:2840898/76‑1 (MQ=255)
                           cttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCtt                                           >  1:5380946/1‑76 (MQ=255)
                            ttctCCTCCTGGTCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCttt                                          <  1:4404334/76‑1 (MQ=255)
                            ttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGGGCGCCAGAATGGCTTCCAGGttt                                          >  2:5135301/1‑76 (MQ=255)
                            ttctCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCttt                                          <  1:2272662/76‑1 (MQ=255)
                             tctcCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGGATGGCTTACAGCtttt                                         >  1:2261380/1‑75 (MQ=255)
                             tctcCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                         <  1:3633925/76‑1 (MQ=255)
                             tctcCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                         <  2:53781/76‑1 (MQ=255)
                             tctcCTCCTGACCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCACGCCAGAATGGCTTCCAGCTTTa                                         <  1:2478191/76‑1 (MQ=255)
                             tctcCACGTGATCCCAGCCGCGGTAGTTCTGTCCGTCACCCTC‑‑TGTGCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                         <  2:5500649/76‑1 (MQ=255)
                             tcccctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTa                                         <  1:2756377/73‑1 (MQ=255)
                              ctcctcCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTaa                                        <  2:690369/76‑1 (MQ=255)
                              ctcctcCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTaa                                        >  1:4608343/1‑76 (MQ=255)
                              ctccgcCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTaa                                        <  2:2056126/76‑1 (MQ=255)
                               tcctccTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAAc                                       <  1:3412608/76‑1 (MQ=255)
                               ccctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAAc                                       <  2:188021/75‑1 (MQ=255)
                                cctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCCGCTTTAACg                                      >  1:451749/1‑76 (MQ=255)
                                cctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACg                                      >  1:5503847/1‑76 (MQ=255)
                                cctcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACg                                      >  1:3978633/1‑76 (MQ=255)
                                  tcctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAAAGGt                                    >  1:5199736/1‑76 (MQ=255)
                                    ctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTTCAGGTTTAAAGGGTc                                  >  2:2485842/1‑76 (MQ=255)
                                    ctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTCACGGtgc                                  >  2:2325637/1‑74 (MQ=255)
                                    ctGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTc                                  >  1:1059140/1‑76 (MQ=255)
                                     tGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTcc                                 <  1:3763244/76‑1 (MQ=255)
                                     tGAGCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTcc                                 <  2:2423922/76‑1 (MQ=255)
                                       aTCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCg                               >  2:595730/1‑76 (MQ=255)
                                         cccAGCCGCGGTAGTTCTGTCCGTAACCTTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTg                             <  1:6229328/76‑1 (MQ=255)
                                          ccAGCCGCGGTATTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGc                            <  2:5809306/76‑1 (MQ=255)
                                           cAGCCGCGTTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGcc                           <  2:3702218/76‑1 (MQ=255)
                                            aGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCCGCTTTAACGGTTCCCATGCCa                          >  2:4017978/1‑76 (MQ=255)
                                            aGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCCGCGTTAACGGGTCCCCTGccc                          >  1:3571812/1‑75 (MQ=255)
                                            aGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCa                          <  2:37371/76‑1 (MQ=255)
                                            aGCCGCGGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCAAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCa                          <  2:5209803/76‑1 (MQ=255)
                                                gcgGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTg                      <  2:1166363/76‑1 (MQ=255)
                                                gcgGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTg                      <  2:6176402/76‑1 (MQ=255)
                                                gcgGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTg                      <  2:6003600/76‑1 (MQ=255)
                                                gcgGTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCCTGCCAGTTg                      >  2:5017229/1‑76 (MQ=255)
                                                  ggTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGAt                    >  2:2563303/1‑76 (MQ=255)
                                                   gTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGCTa                   >  1:1045167/1‑76 (MQ=255)
                                                   gTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATa                   <  2:3578091/76‑1 (MQ=255)
                                                   gTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATa                   >  1:6167122/1‑76 (MQ=255)
                                                   gTAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATa                   >  1:4321085/1‑76 (MQ=255)
                                                    tAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAg                  >  1:4611018/1‑76 (MQ=255)
                                                    tAGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAg                  <  1:1611016/76‑1 (MQ=255)
                                                     aGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCCGTTGATAGt                 >  2:4815735/1‑76 (MQ=255)
                                                     aGTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGGTGATAgg                 >  2:480062/1‑75 (MQ=255)
                                                      gTTCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtc                <  1:6263226/76‑1 (MQ=255)
                                                       ttCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtca               <  2:2995897/76‑1 (MQ=255)
                                                       ttCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtca               <  1:2316567/76‑1 (MQ=255)
                                                        tCTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtcat             >  1:4391185/1‑76 (MQ=255)
                                                         cTGTCCGTAGCCCTCTTTGCGCGCGCGCAC‑‑GAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGtcatc            <  2:2553322/76‑1 (MQ=255)
                                                         cTGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCACTTGATAGtaatt            >  1:3382951/1‑75 (MQ=255)
                                                          tGTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGGGCCAGTTGATAGgcagcc            >  1:450493/1‑73 (MQ=255)
                                                           gTCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCAc          >  2:1891068/1‑76 (MQ=255)
                                                            tCCGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCAcc         <  2:4052789/76‑1 (MQ=255)
                                                             ccGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCg        >  2:4709590/1‑76 (MQ=255)
                                                              cGTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCgg       <  1:4544313/76‑1 (MQ=255)
                                                               gTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCggg      <  1:4708967/76‑1 (MQ=255)
                                                               gTAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGCTGATAGTCATCACCggg      >  1:3279573/1‑76 (MQ=255)
                                                                tAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCgggg     <  1:1655260/76‑1 (MQ=255)
                                                                tAACCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGGGCCAGTTGATAGTCATCACCgggg     >  2:1491408/1‑76 (MQ=255)
                                                                 aaCCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCCGTTGATAGTCATCACCGGGGc    >  2:6050504/1‑76 (MQ=255)
                                                                 aaCCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCCGGTTATAGGCATTCACCgggg     >  2:4802298/1‑76 (MQ=255)
                                                                 aaCCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGc    <  2:3324990/76‑1 (MQ=255)
                                                                 aaCCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGc    >  1:4582667/1‑76 (MQ=255)
                                                                 aaCCCTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAAAGGTTCCCCTGCCAGTTGATAATTATCACCGGGGc    >  2:1401506/1‑76 (MQ=255)
                                                                   cccTCTTT‑‑GCGCGCGCGCCTGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCgg  <  2:506194/76‑1 (MQ=255)
                                                                   cccTCTTT‑‑GCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCgg  <  2:2116591/76‑1 (MQ=255)
                                                                   cccTCTTT‑‑GCGCGCGCGCAAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCgg  <  2:5073197/76‑1 (MQ=255)
                                                                           ||                                                                     
GCGCAGCGGTACGGCGATAGAGGCGATCTTCTCCTCCTGATCCCAGCCGCGGTAGTTCTGTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGG  >  REL606/340323‑340467

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