Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,535,528 (GCTG)3→4 coding (448/597 nt) yhhF → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,535,5151.G84.4% 223.2 / 39.0 64coding (435/597 nt)yhhFpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/10);  new base G (0/54);  total (0/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-03
*REL6063,535,5152.G84.4% 220.3 / 39.0 64coding (435/597 nt)yhhFpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/10);  new base G (0/54);  total (0/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-03
*REL6063,535,5153.C84.4% 212.9 / 39.9 64coding (435/597 nt)yhhFpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/10);  new base C (0/54);  total (0/64)
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.50e-04
*REL6063,535,5154.T84.4% 226.4 / 34.9 64coding (435/597 nt)yhhFpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/10);  new base T (0/54);  total (0/64)
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.07e-02

Only 100 of 108 total aligned reads displayed.

GCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑GGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTG  >  REL606/3535440‑3535572
                                                                            ||||                                                             
gCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAc                                                                   <  2:9188575/76‑1 (MQ=255)
gCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAc                                                                   >  2:10281449/1‑76 (MQ=255)
gCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAc                                                                   <  1:7370545/76‑1 (MQ=255)
gCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAc                                                                   <  1:4489621/76‑1 (MQ=255)
 cATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTTGGAGAGACCATAAATTAACCGCAAGATAAC‑‑‑‑g                                                              >  1:2745170/1‑76 (MQ=255)
 cATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑g                                                              <  2:643387/76‑1 (MQ=255)
 cATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑g                                                              <  2:5886279/76‑1 (MQ=255)
 cATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑g                                                              <  1:4544257/76‑1 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTTGAAGAGACGATAAATTTACTGGAAGATAAA‑‑‑‑gg                                                             >  1:4391107/1‑76 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTTCTGGAAGATAAA‑‑‑‑gg                                                             >  2:771671/1‑76 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑gg                                                             >  1:2484745/1‑76 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑gg                                                             <  2:7496585/76‑1 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑gg                                                             >  1:2421106/1‑76 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑gg                                                             >  2:3509373/1‑76 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑gg                                                             <  1:7205247/76‑1 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAA‑‑‑‑gg                                                             >  2:4095784/1‑76 (MQ=255)
  ataataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAAAGACGATAAAATTACTGGAAGATAAA‑‑‑‑gg                                                             >  1:2039821/1‑76 (MQ=255)
   taataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTTGAAGAGACGATAAATTTACTGGAAGATAAA‑‑‑‑ggc                                                            >  2:8044740/1‑76 (MQ=255)
   taataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggc                                                            <  1:5915535/76‑1 (MQ=255)
   taataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggc                                                            <  1:1122490/76‑1 (MQ=255)
   taataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑agc                                                            >  2:2977257/1‑76 (MQ=255)
    aataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggct                                                           >  2:10194909/1‑76 (MQ=255)
    aataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggct                                                           >  1:5499862/1‑76 (MQ=255)
    aataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggct                                                           >  1:4519573/1‑76 (MQ=255)
    aataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGAAGATAAATTTACTGGAAGATAAC‑‑‑‑ggct                                                           >  1:2696815/1‑76 (MQ=255)
    aataTCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAAACGATAAATTTACTGGGAGATAAC‑‑‑‑gggt                                                           >  1:1844021/1‑74 (MQ=255)
     atatCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctg                                                          <  1:7566825/76‑1 (MQ=255)
     atatCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctg                                                          <  1:9728169/76‑1 (MQ=255)
       atCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTTGAAGAGACGATAAAATTACTGGAAGAGAAA‑‑‑‑ggcgggg                                                        >  1:1522635/1‑75 (MQ=255)
       atCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctggg                                                        >  1:10280665/1‑75 (MQ=255)
       atCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑c                                                    <  1:9684819/76‑1 (MQ=255)
       atCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑c                                                    <  1:8144648/76‑1 (MQ=255)
       atCGTGTTTGTCGATCCACCGGTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑c                                                    <  2:2528746/76‑1 (MQ=255)
        tCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑ct                                                   <  2:7411347/76‑1 (MQ=255)
        tCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑ct                                                   <  1:174206/76‑1 (MQ=255)
         cgggtTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑ctg                                                  <  1:9192339/73‑1 (MQ=255)
         cGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctggggg                                                      >  1:4219669/1‑73 (MQ=255)
         cGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctggggg                                                      >  1:9093288/1‑73 (MQ=255)
         cGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑ctg                                                  <  1:6165137/76‑1 (MQ=255)
         cGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAAAGACGATAAATTTACTGGAAGATAAA‑‑‑‑ggctggggg                                                      >  1:9388026/1‑73 (MQ=255)
            gtTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑ggctgg‑‑‑‑ctggct                                               <  2:8063797/76‑1 (MQ=255)
             tttGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctg                                              <  1:9354754/76‑1 (MQ=255)
             tttGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctg                                              <  1:7155633/76‑1 (MQ=255)
             tttGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctg                                              <  1:2866737/76‑1 (MQ=255)
             tttGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctg                                              <  1:5409581/76‑1 (MQ=255)
              ttGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctgg                                             <  1:5010270/76‑2 (MQ=255)
               tGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctggc                                            <  2:5531884/76‑3 (MQ=255)
                gTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctggct                                           <  2:3820240/76‑4 (MQ=255)
                gTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑Ggctgg‑‑‑‑ctggctggct                                           <  2:2539676/76‑4 (MQ=255)
                               cgccgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGAtt                                <  1:563084/76‑1 (MQ=255)
                                gccgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGAttt                               <  1:9430896/76‑1 (MQ=255)
                                gccgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGAttt                               <  1:6763137/76‑1 (MQ=255)
                                gccgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGAttt                               <  1:1789077/76‑1 (MQ=255)
                                 ccgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTa                              <  2:2627680/76‑1 (MQ=255)
                                  cgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTAt                             <  1:1204323/76‑1 (MQ=255)
                                  cgTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTAt                             <  1:2160810/76‑1 (MQ=255)
                                     ggCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTc                          <  2:676284/76‑1 (MQ=255)
                                     ggCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTc                          <  2:8781346/76‑1 (MQ=255)
                                     ggCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTc                          <  2:418783/76‑1 (MQ=255)
                                      gCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCg                         <  1:8546458/76‑1 (MQ=255)
                                      gCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCg                         <  1:1200578/76‑1 (MQ=255)
                                       cTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGa                        <  1:3421148/76‑1 (MQ=255)
                                       cTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGa                        <  1:3308498/76‑1 (MQ=255)
                                       cTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGa                        <  1:5197108/76‑1 (MQ=255)
                                         tGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGaaa                      <  1:7309646/76‑1 (MQ=255)
                                         tGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGaaa                      <  2:9534272/76‑1 (MQ=255)
                                         tGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGaaa                      <  2:4405548/76‑1 (MQ=255)
                                          gTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAg                     <  2:6351561/76‑1 (MQ=255)
                                          gTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAg                     <  1:6419679/76‑1 (MQ=255)
                                          gTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAg                     <  2:406553/76‑1 (MQ=255)
                                          gTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAg                     <  1:6050657/76‑1 (MQ=255)
                                          gTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAg                     <  2:5256648/76‑1 (MQ=255)
                                          gTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑GGCTGGCTGTCTGGCTGACGAAGCCCTGATTTATGTCGAAAg                     <  1:6831454/76‑1 (MQ=255)
                                           ttAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGc                    <  2:5899365/76‑1 (MQ=255)
                                            tAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCg                   <  1:3448933/76‑1 (MQ=255)
                                            tAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCg                   <  1:4267474/76‑1 (MQ=255)
                                             agaagaGACGATAAATTTACTGGAAGTTAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGa                  <  1:8666168/76‑1 (MQ=255)
                                             agaagaGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGa                  <  1:2851965/76‑1 (MQ=255)
                                             agaagaGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGa                  <  2:850604/76‑1 (MQ=255)
                                              gaagaGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGaa                 <  1:4251997/76‑1 (MQ=255)
                                               aagaGACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACCAAGCCCTGATTTATGTCGAAAGCGAAg                <  1:2917772/76‑1 (MQ=255)
                                                agagACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGt               <  2:5266168/76‑1 (MQ=255)
                                                agagACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGt               <  1:6165440/76‑1 (MQ=255)
                                                agagACGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGt               <  1:8248725/76‑1 (MQ=255)
                                                   gaCGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGa            <  1:9949265/76‑1 (MQ=255)
                                                    aCGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaa           <  2:8278839/76‑1 (MQ=255)
                                                    aCGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaa           <  2:8829308/76‑1 (MQ=255)
                                                    aCGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaa           <  1:9778457/76‑1 (MQ=255)
                                                     cGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaa          <  2:7857832/76‑1 (MQ=255)
                                                     cGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaa          <  2:5430237/76‑1 (MQ=255)
                                                     cGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaa          <  1:6399461/76‑1 (MQ=255)
                                                     cGATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaa          <  1:6748399/76‑1 (MQ=255)
                                                      gATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaaa         <  1:3312542/76‑1 (MQ=255)
                                                      gATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaaa         <  2:5962454/76‑1 (MQ=255)
                                                      gATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaaa         <  1:1868220/76‑1 (MQ=255)
                                                      gATAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGaaaa         <  2:2629721/76‑1 (MQ=255)
                                                        tAAATTTACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACg       <  2:672538/76‑1 (MQ=255)
                                                            tttACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCt   <  2:5499862/76‑1 (MQ=255)
                                                            tttACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCt   <  2:2738935/76‑1 (MQ=255)
                                                             ttACTGGAAGATAACGGCTGGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGCCGAAAACGGTCTg  <  2:5744177/76‑1 (MQ=255)
                                                                            ||||                                                             
GCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAAC‑‑‑‑GGCTGG‑‑‑‑CTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTG  >  REL606/3535440‑3535572

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