Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,468 A→T 18.2% Q385L (CAG→CTG)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4680AT18.2% 335.6 / 71.2 149Q385L (CAG→CTG) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base A (56/66);  new base T (13/14);  total (69/80)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-02

Only 100 of 155 total aligned reads displayed.

TTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666369‑1666567
                                                                                                   |                                                                                                   
ttttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCa                                                                                                     >  1:1057513/1‑100 (MQ=255)
  ttACTCGGATTGACCGGGATGGCGACGCTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAAAGACCGGGGAGTATCAga                                                                                                   <  2:122030/100‑2 (MQ=255)
  ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                   <  2:2605883/100‑1 (MQ=255)
   tACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCt                                                                                                  <  2:937744/100‑1 (MQ=255)
    aCTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                 >  1:2238625/1‑100 (MQ=255)
     cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                <  2:2801240/100‑1 (MQ=255)
     cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  1:2583673/1‑100 (MQ=255)
     cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  2:33502/1‑100 (MQ=255)
     cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  2:1935057/1‑100 (MQ=255)
        ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                             <  2:250415/100‑1 (MQ=255)
         gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTccc                                                                                            >  2:2480225/1‑100 (MQ=255)
          cgatctcGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                           <  1:431835/94‑1 (MQ=255)
            tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:2823723/100‑1 (MQ=39)
            tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:2445411/100‑1 (MQ=39)
            tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:1849460/100‑1 (MQ=39)
            tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:796684/100‑1 (MQ=39)
            tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:99022/100‑1 (MQ=39)
            tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:2846792/100‑1 (MQ=39)
             gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                        >  1:462168/1‑100 (MQ=255)
             gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                        <  1:2682082/100‑1 (MQ=255)
             gACCGGGATGGCGACACTGTTAGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGCCCGGGGAGTAGCTGATACTCCCCGGt                                                                                        <  2:2358916/100‑1 (MQ=38)
             gACCGGGATGGCGACACTGGTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                        <  2:1047948/100‑1 (MQ=255)
              aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                       >  2:1729373/1‑100 (MQ=255)
              aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                       <  1:1403635/100‑1 (MQ=255)
               ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                      >  1:3198078/1‑100 (MQ=255)
                cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTccccagatcg                                                                                     >  2:431833/1‑94 (MQ=255)
                  ggATCGCTACACTGTTGGCGGCGACCACGCACGTGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcattt                                                                                   <  2:3110683/100‑6 (MQ=25)
                   gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                  <  1:1542482/100‑1 (MQ=255)
                      ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                               >  1:2271147/1‑100 (MQ=255)
                       gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                              <  1:2532310/100‑1 (MQ=255)
                        cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                             >  2:1866101/1‑100 (MQ=255)
                        cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                             >  2:3031902/1‑100 (MQ=255)
                        cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTACTGATTGc                                                                             >  1:1753467/1‑100 (MQ=255)
                        cGACACTGTTAGCGGCGACCACGAACGCGCCGATTATCTCGACGTTGATGATATGGGAAATGACCGGGGAGTGTCAGCGACTCCCCGGTTGATTGATTGc                                                                             >  2:3539441/1‑100 (MQ=255)
                          acacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                           <  2:3461903/100‑1 (MQ=255)
                             cTGTTGGCGGCGACCACGCACGCGCCGAGTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCTCCGGTTTATTGATTGCCTGCg                                                                        <  2:1679164/100‑1 (MQ=255)
                             cTGTTGGCGGCGACCACGCACGCGCAGATTATGTAGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                        >  2:2505343/1‑100 (MQ=255)
                               gTTGGCGGCGACCACGCACGCGCCGCTTTTGTCGACGTTGATGAGATGTGCAATGACCGGGGCGTCTCATCTACTCGCCGGGGTATCGGTTGCCTGc                                                                         >  2:1129896/1‑97 (MQ=255)
                                ttGGCGGCGACAACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                      <  1:386152/99‑1 (MQ=255)
                                 tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                    <  1:1418384/100‑1 (MQ=255)
                                  ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                   >  1:1075787/1‑100 (MQ=255)
                                  ggcggcGACCACGCACGCGCCGAATATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                   >  1:2004592/1‑100 (MQ=255)
                                     ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGc                                                                 >  2:3580342/1‑99 (MQ=255)
                                       cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              >  1:3577275/1‑100 (MQ=255)
                                         aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGCTTTATTGATTGCCTGCGTAATTGCGTc                                                              <  1:740177/98‑1 (MQ=255)
                                            acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTaa                                                         >  2:2773346/1‑100 (MQ=255)
                                                  ctgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                   <  1:1093088/98‑1 (MQ=255)
                                                   cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                  >  2:723477/1‑100 (MQ=255)
                                                   cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                  <  1:628219/100‑1 (MQ=255)
                                                    gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcag                                                 >  2:1093080/1‑98 (MQ=255)
                                                    gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  1:1777926/100‑1 (MQ=255)
                                                    gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  2:1170345/100‑1 (MQ=255)
                                                    gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 >  2:1928474/1‑100 (MQ=255)
                                                       gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                              <  1:3513018/100‑1 (MQ=255)
                                                         ttATGTCGACGTTGATGATATGTAAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCCGACg                                              <  2:350917/98‑1 (MQ=255)
                                                           aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                          <  1:723480/100‑1 (MQ=255)
                                                            tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                         >  1:123174/1‑100 (MQ=255)
                                                             ntCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTacac                                        >  2:1883868/2‑100 (MQ=255)
                                                                 aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                    <  2:2524878/100‑1 (MQ=255)
                                                                   gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                     <  2:2446251/97‑1 (MQ=255)
                                                                   gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACAGTGAc                                  >  2:2423294/1‑100 (MQ=255)
                                                                    ttGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                   >  2:1631716/1‑98 (MQ=255)
                                                                       atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                              <  1:3100776/100‑1 (MQ=255)
                                                                        tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                             >  1:922013/1‑100 (MQ=255)
                                                                          atatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                           <  2:1306744/100‑1 (MQ=255)
                                                                           naggtgAAATGACCGGGGAGTATCAGCGACTCCCCGGGTTATGGAGTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGGTACACCGTGACtctacct                           >  2:729543/4‑99 (MQ=255)
                                                                             tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  2:2072899/100‑1 (MQ=255)
                                                                             tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  1:1938293/100‑1 (MQ=255)
                                                                             tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  1:576732/100‑1 (MQ=255)
                                                                              gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                       >  1:391872/1‑100 (MQ=255)
                                                                               tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                      >  1:1838276/1‑100 (MQ=255)
                                                                               tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                      <  1:3027030/100‑1 (MQ=255)
                                                                                gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     <  1:998029/100‑1 (MQ=255)
                                                                                gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  1:1510903/1‑100 (MQ=255)
                                                                                 aataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                    >  1:330804/6‑100 (MQ=37)
                                                                                 aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACATTACACCGTGACTCTATCTACCGCCa                    <  1:54165/100‑1 (MQ=255)
                                                                                   ataaCCGGGGAGTATCATCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                  <  2:3001909/97‑1 (MQ=255)
                                                                                    aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                 >  2:14476/3‑100 (MQ=38)
                                                                                     gACCGGGGAGTATCAGCTACTCCCCGCTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGGCGTTACACCGTGACTCTATCTACCGCCAGcac                >  2:2304932/1‑100 (MQ=255)
                                                                                      aCCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcaca               >  1:2163667/1‑100 (MQ=255)
                                                                                       ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:3136594/1‑100 (MQ=255)
                                                                                       ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              <  1:1659659/100‑1 (MQ=255)
                                                                                       ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              <  1:3528976/100‑1 (MQ=255)
                                                                                       ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACAGTGACTCTATCTACCGCCAGCacac              >  1:2674026/1‑100 (MQ=255)
                                                                                       ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTCATTGCGTCGGTAATTTCGCGGACGATCCACCGTGACTCTATCTACCGCCAGCacac              >  2:2544345/1‑100 (MQ=38)
                                                                                       ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:2821277/1‑100 (MQ=39)
                                                                                       ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:1057699/1‑100 (MQ=39)
                                                                                        cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt             >  2:2009518/1‑100 (MQ=255)
                                                                                        cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt             >  1:2283682/1‑100 (MQ=255)
                                                                                         ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  1:3337216/1‑100 (MQ=255)
                                                                                         ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  2:677419/1‑100 (MQ=255)
                                                                                             aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcgga        <  1:2462280/100‑3 (MQ=255)
                                                                                              gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       >  1:2263831/1‑100 (MQ=255)
                                                                                              gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       <  1:1176425/100‑1 (MQ=255)
                                                                                              gTATCAGCTACTCCCCGGTTTATTGATTGCCGGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       <  2:468619/100‑1 (MQ=255)
                                                                                               tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagc      <  2:391870/100‑1 (MQ=255)
                                                                                                 tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt    >  1:3123901/1‑100 (MQ=255)
                                                                                                   aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTACTTTCGCGGACGTTACACCGTGACGCTATCTACCGCCGGCACACTGCGCAGCAtcg  >  2:2225374/1‑98 (MQ=255)
                                                                                                   aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCCGCACACTCCGCAGCATAg  >  2:1777807/1‑100 (MQ=255)
                                                                                                   aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg  >  2:213452/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666369‑1666567

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