Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,817,619 A→G 5.2% F294F (TTT→TTC)  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,6190AG 5.2% 238.8 / 2.4 76F294F (TTT→TTC) chbFphospho‑chitobiase; general 6‑phospho‑beta‑glucosidase activity
Reads supporting (aligned to +/- strand):  ref base A (39/33);  new base G (2/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCAACGAACGAATTAAGCC  >  NC_000913/1817475‑1817765
                                                                                                                                                |                                                                                                                                                  
atgCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCaaaaa                                                                                                                                                <  5:257450/149‑1 (MQ=255)
    gggATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAgttg                                                                                                                                            <  1:326938/149‑1 (MQ=255)
    gggATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAgttg                                                                                                                                            >  8:263940/1‑149 (MQ=255)
    gggATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAgttg                                                                                                                                            >  8:257845/1‑149 (MQ=255)
    gggATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAgttg                                                                                                                                            >  2:133857/1‑149 (MQ=255)
          ttAACGTAATGTTCTGCGTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCCTGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTGCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGAGTTTTATACAGCTCAAAAAGTTGTTtc                                                                                                                                        <  4:340802/147‑1 (MQ=255)
            aaCGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCg                                                                                                                                     <  1:271816/148‑1 (MQ=255)
            aaCGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGa                                                                                                                                    >  3:136083/1‑149 (MQ=255)
             aCGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGATGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGa                                                                                                                                    <  6:7791/148‑1 (MQ=255)
                  aTGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCt                                                                                                                              <  6:199269/149‑1 (MQ=255)
                     ttCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGta                                                                                                                           >  2:36666/1‑149 (MQ=255)
                          cttgcttgTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGtactac                                                                                                                       <  5:189692/148‑1 (MQ=255)
                          cttgcttgTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGtactac                                                                                                                       <  8:1346/148‑1 (MQ=255)
                             gcttgTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGATGTTCAAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCtcc                                                                                                                   >  4:174665/1‑147 (MQ=255)
                                 gTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTc                                                                                                               >  4:286107/1‑149 (MQ=255)
                                  tCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGGTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTGTCTGTACTACCTGCGCTCg                                                                                                              >  8:310308/1‑149 (MQ=255)
                                    gTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTgcgcgagtg                                                                                                            >  2:82370/1‑144 (MQ=255)
                                    gTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTg                                                                                                            >  4:141112/1‑149 (MQ=255)
                                         aGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGcgccg                                                                                                       >  8:378222/1‑149 (MQ=255)
                                          gATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACttt                                                                                                                                >  3:174788/1‑123 (MQ=255)
                                          gATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACttt                                                                                                                                <  4:174788/123‑1 (MQ=255)
                                          gATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGcgccgc                                                                                                      <  2:37063/149‑1 (MQ=255)
                                              gCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCGAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCttt                                                                                                  >  1:157585/1‑149 (MQ=255)
                                               cGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTg                                                                                                 <  8:320797/149‑1 (MQ=255)
                                                 ttGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCGAAAAGTTGTTTCTCGACTTTCt                                                                                                                              <  3:113374/118‑1 (MQ=255)
                                                 ttGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCGAAAAGTTGTTTCTCGACTTTCt                                                                                                                              >  4:113374/1‑118 (MQ=255)
                                                 ttGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgta                                                                                               <  8:395084/149‑1 (MQ=255)
                                                 ttGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgta                                                                                               <  1:133857/149‑1 (MQ=255)
                                                 ttGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTATTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTAATAACTGCGCTCGTGCTCCGCCGTTg                                                                                                 >  2:343771/1‑147 (MQ=255)
                                                   gATCACTTCGCATGCTGCATCAGAGTAATAATCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAAGTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTTCCTGCGCTCGTGGGGCGCCTTTgtagt                                                                                             >  3:303509/1‑149 (MQ=255)
                                                   gATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgtag                                                                                              >  8:171268/1‑148 (MQ=255)
                                                   gATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgtag                                                                                              >  3:308724/1‑148 (MQ=255)
                                                    aTCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgtagt                                                                                             <  4:16482/148‑1 (MQ=255)
                                                       aCTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGGTCCAGTGCTTTCGGCTTAACTTTCAGCTCAGGATGTTTATACCGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTTCGCTCGTGCGCCGCCTTTGTAGTAtt                                                                                          <  5:150790/148‑1 (MQ=255)
                                                             cATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCAt                                                                                   >  8:6836/1‑149 (MQ=255)
                                                             cATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCAt                                                                                   >  3:105024/1‑149 (MQ=255)
                                                               tgctgcATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCAttt                                                                                 <  1:36666/149‑1 (MQ=255)
                                                                 ctgcATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTc                                                                                <  7:378222/148‑1 (MQ=255)
                                                                    cATCAGAGTAATTACGCCCACCGCGCTGTTCCAGTTCTTTCGGCGTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTTTAGTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATa                                                                            <  7:310308/149‑1 (MQ=255)
                                                                    cATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgta                                                                                               <  6:278499/130‑1 (MQ=255)
                                                                    cATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTgta                                                                                               >  5:278499/1‑130 (MQ=255)
                                                                    cATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAAt                                                                             <  7:263940/148‑1 (MQ=255)
                                                                    cATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAAt                                                                             <  7:257845/148‑1 (MQ=255)
                                                                    cATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTGTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCaa                                                                              <  8:50230/147‑1 (MQ=255)
                                                                     aTCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGGTTTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAg                                                                           >  2:68929/1‑149 (MQ=255)
                                                                         gagTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACCGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAg                                                                       >  4:312249/1‑149 (MQ=255)
                                                                         gagTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAg                                                                       >  8:62706/1‑149 (MQ=255)
                                                                            taataaGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCAt                                                                    >  3:224344/1‑149 (MQ=255)
                                                                              ataaGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATa                                                                            >  4:75302/1‑139 (MQ=255)
                                                                              ataaGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATa                                                                            <  3:75302/139‑1 (MQ=255)
                                                                                     cccccGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTg                                                           >  7:252658/4‑149 (MQ=255)
                                                                                       aCCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGACGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCt                                                         >  3:298058/1‑149 (MQ=255)
                                                                                        ccGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCtt                                                        <  3:368598/149‑1 (MQ=255)
                                                                                           cgcTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGaa                                                     <  1:238105/149‑1 (MQ=255)
                                                                                           cgcTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTATATACAGCTCAAAAAGTTGTTTCTCGACTTACTGTACTACCTGCGCTCGTGCGCCGCCTATGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGaa                                                     >  1:201606/1‑149 (MQ=255)
                                                                                            gcTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAg                                                    <  7:171268/149‑1 (MQ=255)
                                                                                                    aGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATAcag                                            <  4:136083/149‑1 (MQ=255)
                                                                                                     gTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATAcag                                            >  8:294324/1‑148 (MQ=255)
                                                                                                      ttCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCGAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATAcagca                                          <  2:157585/149‑1 (MQ=255)
                                                                                                      ttCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATAcagca                                          <  4:209165/149‑1 (MQ=255)
                                                                                                      ttCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATAcag                                            >  5:218139/1‑147 (MQ=255)
                                                                                                            cGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATaa                                     <  6:218139/148‑1 (MQ=255)
                                                                                                               cTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGaa                                  >  7:188631/1‑148 (MQ=255)
                                                                                                               cTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGaa                                  >  6:246038/1‑148 (MQ=255)
                                                                                                                ttAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACa                                >  3:370076/1‑149 (MQ=255)
                                                                                                                   aCTTTCAGCTCTGGATTTTTATACAGCTCAAAAAGTTGTTTCTCTACTTTCTGTTCTACATGCGCTCGTGGGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCTGCTTCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATg                              <  3:122675/148‑1 (MQ=255)
                                                                                                                   aCTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATgg                             <  8:252658/149‑1 (MQ=255)
                                                                                                                     tttCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATgg                             >  1:166741/1‑147 (MQ=255)
                                                                                                                     tttCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATgg                             >  4:131496/1‑147 (MQ=255)
                                                                                                                         aGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGgcagca                        <  1:82370/148‑1 (MQ=255)
                                                                                                                               ggATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCaa                 >  3:82533/1‑149 (MQ=255)
                                                                                                                                gATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATa                                                                            <  4:1540/89‑1 (MQ=255)
                                                                                                                                gATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATa                                                                            >  3:1540/1‑89 (MQ=255)
                                                                                                                                gATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCaac                >  1:85289/1‑149 (MQ=255)
                                                                                                                                           aGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCAACGAACGAAtt       >  7:132377/1‑147 (MQ=255)
                                                                                                                                             cTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCAACGAACGAATTAAGc   <  3:141112/149‑1 (MQ=255)
                                                                                                                                              tCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCAACGAACGAATTAAGcc  >  2:65724/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ATGCGGGATATTAACGTAATGTTCTGCTTGCTTGTCGTTGTAGATAGCGTTGATCACTTCGCATGCTGCATCAGAGTAATAAGCCCCACCGCGCTGTTCCAGTTCTTTCGGCTTAACTTTCAGCTCAGGATTTTTATACAGCTCAAAAAGTTGTTTCTCGACTTTCTGTACTACCTGCGCTCGTGCGCCGCCTTTGTAGTATTCGCCCATTTCAATAGCCAGCATCTCTTTCTGCTTGAAGTAATACAGCAGATAAGAACATGGCAGCAGATTCAACGAACGAATTAAGCC  >  NC_000913/1817475‑1817765

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