Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,817,424 A→G C359C (TGT→TGC)  chbF ← phospho‑chitobiase; general 6‑phospho‑beta‑glucosidase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,817,4240AG98.5% 255.8 / ‑4.3 66C359C (TGT→TGC) chbFphospho‑chitobiase; general 6‑phospho‑beta‑glucosidase activity
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (36/29);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

CAGTAACACATCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTACAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCT  >  NC_000913/1817280‑1817569
                                                                                                                                                |                                                                                                                                                 
cAGTAACACATCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTc                                                                                                                                             <  1:963204/151‑1 (MQ=255)
   tAACACATCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCAtt                                                                                                                                          >  1:576111/1‑151 (MQ=255)
    aaCACATCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCAtt                                                                                                                                          <  1:855124/150‑1 (MQ=255)
    aaCACATCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCAt                                                                                                                                           >  1:702682/1‑149 (MQ=255)
       acaTCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCt                                                                                                                                       <  1:490513/150‑1 (MQ=255)
         aTCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACt                                                                                                                                    >  2:531618/1‑151 (MQ=255)
               aaaTTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGccc                                                                                                                                >  1:597651/1‑149 (MQ=255)
                 aTTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTc                                                                                                                            <  1:642745/151‑1 (MQ=255)
                      ccGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCg                                                                                                                       <  2:393685/151‑1 (MQ=255)
                          tAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAat                                                                                                                   >  1:990381/1‑151 (MQ=255)
                           aaGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAat                                                                                                                   >  1:110763/1‑150 (MQ=255)
                             gTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAAtatt                                                                                                                <  1:771652/151‑1 (MQ=255)
                             gTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAAtatt                                                                                                                <  1:775245/151‑1 (MQ=255)
                              tGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTACAGGTCATTTCTACTGCCCAGTCTGCCGGAAtatt                                                                                                                >  2:1097810/1‑150 (MQ=255)
                                 ggcgTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATc                                                                                                             >  2:1020150/1‑150 (MQ=255)
                                  gcgTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCaa                                                                                                           <  2:492233/151‑1 (MQ=255)
                                   cgTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTAGCGTCGCGCCATCGCGCCGCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATc                                                                                                             <  1:209829/148‑1 (MQ=255)
                                     ttACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCt                                                                                                        >  1:472961/1‑151 (MQ=255)
                                     ttACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCt                                                                                                        >  2:1165891/1‑151 (MQ=255)
                                          ggCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGccc                                                                                                    >  2:1037917/1‑150 (MQ=255)
                                            cagcaATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCatg                                                                                                 <  1:1041115/151‑1 (MQ=255)
                                            cagcaATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCatg                                                                                                 <  2:1094761/151‑1 (MQ=255)
                                               caATCTCGAAGCCTTTAATGGTGTGAATCAACCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCatgatg                                                                                              <  2:668803/151‑1 (MQ=255)
                                                aaTCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGc                                                                                             <  2:1153174/151‑1 (MQ=255)
                                                   ctcGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCggg                                                                                          >  2:650136/1‑151 (MQ=255)
                                                    tcGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCggg                                                                                          >  2:218584/1‑150 (MQ=255)
                                                     cGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGa                                                                                         <  2:617586/150‑1 (MQ=255)
                                                      gAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGat                                                                                        >  2:994194/1‑150 (MQ=255)
                                                          ccTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTa                                                                                    >  2:675710/1‑150 (MQ=255)
                                                              tAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGt                                                                                <  2:1154388/150‑1 (MQ=255)
                                                                 tGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATg                                                                            >  1:386013/1‑151 (MQ=255)
                                                                  ggTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATg                                                                            <  1:987104/150‑1 (MQ=255)
                                                                     gtgAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTc                                                                         <  2:116681/150‑1 (MQ=255)
                                                                       gAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCtgc                                                                      >  1:263727/1‑151 (MQ=255)
                                                                        aaTCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGTTGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAGTCTGCCCATGATGCGGGATATTAACGTAACGTTCtgc                                                                      >  1:239622/1‑150 (MQ=255)
                                                                         aTCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCtgctt                                                                    >  1:437955/1‑151 (MQ=255)
                                                                           cAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCtgcttg                                                                   <  1:1169291/150‑1 (MQ=255)
                                                                                  aTTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGt                                                           >  2:759570/1‑151 (MQ=255)
                                                                                  aTTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGt                                                           <  2:334693/151‑1 (MQ=255)
                                                                                        ttATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGTCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGt                                                        >  2:115481/1‑148 (MQ=255)
                                                                                        ttATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGt                                                        >  2:81624/1‑148 (MQ=255)
                                                                                        ttATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTa                                                       <  1:856257/149‑1 (MQ=255)
                                                                                         tATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGt                                                        >  1:777276/1‑147 (MQ=255)
                                                                                            catcGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGc                                                 <  1:675687/151‑1 (MQ=255)
                                                                                               cGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGtt                                              >  2:219310/1‑151 (MQ=255)
                                                                                                gAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTg                                             >  2:596984/1‑151 (MQ=255)
                                                                                                  aaTGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGAt                                           >  2:934611/1‑151 (MQ=255)
                                                                                                  aaTGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGAt                                           <  2:1053024/151‑1 (MQ=255)
                                                                                                   aTGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATc                                          <  2:935996/151‑1 (MQ=255)
                                                                                                   aTGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATc                                          >  1:817309/1‑151 (MQ=255)
                                                                                                    tGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATc                                          >  2:733771/1‑150 (MQ=255)
                                                                                                      cGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACt                                       >  2:939526/1‑151 (MQ=255)
                                                                                                       gTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACtt                                      >  2:1148365/1‑151 (MQ=255)
                                                                                                         aaTGCGCGGATGTGGCGTCGCGCCATCGCTCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTc                                     <  2:660736/150‑1 (MQ=255)
                                                                                                         aaTGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTc                                     <  2:990441/150‑1 (MQ=255)
                                                                                                             cgcgGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCAtg                                <  2:642643/151‑1 (MQ=255)
                                                                                                               cgGATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCAtgct                              >  1:939534/1‑151 (MQ=255)
                                                                                                                 gATGTGGCGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCAtgctgc                            >  1:245265/1‑151 (MQ=255)
                                                                                                                        cGTCGCGCCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCag                       >  2:568563/1‑149 (MQ=255)
                                                                                                                             cgcCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGtaat                 >  2:169054/1‑150 (MQ=255)
                                                                                                                              gcCATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCCCATGCTGCATCAGAGtaataa               <  1:220011/151‑1 (MQ=255)
                                                                                                                               ccATCGCGCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGtaata                >  1:525290/1‑149 (MQ=255)
                                                                                                                                   cgcgCCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGcccc          <  2:1072898/151‑1 (MQ=255)
                                                                                                                                     cgcCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCAc        <  1:162285/151‑1 (MQ=255)
                                                                                                                                     cgcCCCAGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCAc        >  2:188748/1‑151 (MQ=255)
                                                                                                                                           aGCTTGCAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCt  >  2:120789/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CAGTAACACATCGTTAAATTCTCCGCTAAGTGCGGCGTTACTGGCAGCAATCTCGAAGCCTTTAATGGTGTGAATCAGCCCCATTACTTTATCATCGAAATGCGTAATGCGCGGATGTGGCGTCGCGCCATCGCGCCCCAGCTTACAGGTCATTTCTACTGCCCAGTCTGCCGGAATATTATCAATCTGCCCATGATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCT  >  NC_000913/1817280‑1817569

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