Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,468 A→T 15.5% 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,4680AT15.5% 367.8 / 54.2 147Q385L (CAG→CTG) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base A (58/66);  new base T (10/13);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 148 total aligned reads displayed.

TTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666370‑1666567
                                                                                                  |                                                                                                   
tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                    <  1:1492588/100‑1 (MQ=255)
 ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                   <  1:2617158/100‑1 (MQ=255)
 ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                   <  2:3018923/100‑1 (MQ=255)
  tACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCt                                                                                                  <  2:2007778/100‑1 (MQ=255)
    cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  1:1338809/1‑100 (MQ=255)
    cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  1:710693/1‑100 (MQ=255)
       ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                             <  1:3397676/100‑1 (MQ=255)
       ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                             <  1:2072500/100‑1 (MQ=255)
        gATTGACCGGGATGGCGGCAGTGTGGGAGGCGACCACGCACGCGCCGATTATGTCGGCGTTGATGTTATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                              <  2:927893/98‑1 (MQ=255)
        gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTccc                                                                                            <  1:2123536/100‑1 (MQ=255)
         aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                           <  1:410776/100‑1 (MQ=255)
           tgacatggATGGTTATACTGTTGGCGGCGAAAACGCACGGGCGGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATTCTCCCCgg                                                                                         <  1:813545/94‑1 (MQ=255)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATAAGCTACTCCCCgg                                                                                         >  1:870836/1‑100 (MQ=255)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:2638137/100‑1 (MQ=39)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:49475/100‑1 (MQ=39)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGACGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:1791397/100‑1 (MQ=38)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGCTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                        <  2:2503707/100‑1 (MQ=38)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                         >  1:2440462/1‑99 (MQ=255)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                        >  2:1250281/1‑100 (MQ=255)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                        <  2:1675572/100‑1 (MQ=255)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAATATCAGCTACTCCCCGGt                                                                                        >  1:2322261/1‑100 (MQ=255)
            cttccgatctgCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                        <  1:3737986/90‑1 (MQ=37)
               cGGTATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATAGGTGAAATCACTGGTGTGTAGCTGATACTCCCCGGtcat                                                                                     <  2:2418281/100‑4 (MQ=255)
               cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTa                                                                                     <  2:1720208/100‑1 (MQ=255)
                 ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                   >  2:1478528/1‑100 (MQ=255)
                   aTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattga                                                                                 <  1:2469197/100‑1 (MQ=255)
                    tGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                                <  2:424/100‑1 (MQ=255)
                     ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                               >  1:1881742/1‑100 (MQ=255)
                      gcgcCACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                               >  2:1210466/1‑99 (MQ=255)
                      gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGACATGACCGGGGAGTATCAGCTACTCCCCGGTTCattgatt                                                                               >  2:413198/1‑99 (MQ=255)
                      gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                              >  1:1755804/1‑100 (MQ=255)
                      gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                              <  2:2974354/100‑1 (MQ=255)
                      gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                              <  2:2611915/100‑1 (MQ=255)
                      gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtagatcggaag                                                                              >  2:3737947/1‑90 (MQ=38)
                       cGACACTGTTCGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                             <  2:2734964/100‑1 (MQ=255)
                         acacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                           >  2:3197988/1‑100 (MQ=255)
                          cacTGTTGGCGGCGACCACGCACGCGTCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                          <  1:447197/100‑1 (MQ=255)
                          cacTGTTGGCGGCGACCACGCACGAGCAAATTATGTCGACGTAGATGATATGTGAAATTACCGGGGAGTATCAGTTACTCCCCAGTTTattgaatgac                                                                            >  2:3750064/1‑96 (MQ=255)
                            cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                        <  1:397546/100‑1 (MQ=255)
                             tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGt                                                                       >  1:2898270/1‑100 (MQ=255)
                                 ggcggcGACCACGCACGCGCAGATTATGTCGACGTTGATGATATTTGAAATTACCGGGGAGTATCAGCTACTCCCTGGTTTATTGATTGCCTGCGTAAt                                                                    <  2:349708/99‑1 (MQ=255)
                                  gcggctACCACGCACGCGCCTATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTTCCTGCGTAATTg                                                                  <  1:1139798/100‑1 (MQ=255)
                                  gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAAAATAAGATACTCCACGGTTTATTGATTGCCTGCGTAATTg                                                                  >  1:2214535/1‑100 (MQ=255)
                                      cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              <  1:3427785/100‑1 (MQ=255)
                                      cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              <  1:2582791/100‑1 (MQ=255)
                                       gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCg                                                             <  2:3858676/100‑1 (MQ=255)
                                        aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                            <  1:494798/100‑1 (MQ=255)
                                        aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                            <  2:2044913/100‑1 (MQ=255)
                                        aCCACGCACGCGCCGAATATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                            >  2:1644880/1‑100 (MQ=255)
                                                cgcgCCGATTATGTTGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                    <  1:3302282/100‑1 (MQ=255)
                                                 gcgcCGTTTATGTCGACGATGATGAAATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTAGGTGATTTcgc                                                   >  2:1271382/1‑100 (MQ=255)
                                                  cgcCGATTATGTGGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                  <  1:3624670/100‑1 (MQ=255)
                                                  cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                  <  1:3871530/100‑1 (MQ=255)
                                                  cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                  <  1:181846/100‑1 (MQ=255)
                                                   gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  1:2493708/100‑1 (MQ=255)
                                                   gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 >  2:2218227/1‑100 (MQ=255)
                                                     cGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                               >  2:2917775/1‑100 (MQ=255)
                                                     cGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                               >  2:3847950/1‑100 (MQ=255)
                                                      gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                              >  2:2395813/1‑100 (MQ=255)
                                                       aTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGt                                             <  1:2426299/100‑1 (MQ=255)
                                                          aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                          >  1:146019/1‑100 (MQ=255)
                                                           tGTCGACGTTGATGATATTTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                         <  1:1250295/100‑1 (MQ=255)
                                                             tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAcc                                       <  2:1400923/100‑1 (MQ=255)
                                                              cGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                      >  1:2227329/1‑100 (MQ=255)
                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                  <  2:870825/100‑1 (MQ=255)
                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                  <  2:3555053/100‑1 (MQ=255)
                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                  <  1:283735/100‑1 (MQ=255)
                                                                  gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                  <  1:1381155/100‑1 (MQ=255)
                                                                     gatgatATGTGAAATGACCGGGGAGTATCAGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAAAGTGActct                               >  2:195812/1‑100 (MQ=255)
                                                                      atgatATGTGAAATGACCGGGGAGTATAAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                              <  2:3576951/100‑1 (MQ=255)
                                                                       tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                             >  2:1842569/1‑100 (MQ=255)
                                                                            tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        >  2:58049/1‑100 (MQ=255)
                                                                             gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                       >  2:3150639/1‑100 (MQ=255)
                                                                             gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                       >  1:1251018/1‑100 (MQ=255)
                                                                             gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCAAACGTTACACCGTGACTCTATCTAcc                        >  1:893951/1‑99 (MQ=255)
                                                                             gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGACGGTAATTTCGCAGACGTTACACCGTGACTCTATCTAcc                        >  1:710984/1‑99 (MQ=255)
                                                                             gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATAGCGTCGGTAATATCGCGGACGTTACACCGTGACTCTATCTACCg                       >  1:121773/1‑100 (MQ=255)
                                                                               gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  2:968237/1‑100 (MQ=255)
                                                                                 aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                   <  2:1881727/100‑1 (MQ=255)
                                                                                  aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                  <  2:1755789/100‑1 (MQ=255)
                                                                                    gACCGGGGAGTATCAGCTACTCCCTGGTTTATTGATTGTCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcac                <  1:3660833/100‑1 (MQ=255)
                                                                                      ncGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  2:107952/2‑100 (MQ=255)
                                                                                      ccGGGGAGTAGCTGATACTCCACAGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTAACACCGTGACTCTATCTACCGCCAGCagac              >  2:229814/1‑100 (MQ=255)
                                                                                       cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt             >  1:2557051/1‑100 (MQ=255)
                                                                                       cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt             <  2:1210781/100‑1 (MQ=255)
                                                                                       cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  2:1385326/1‑99 (MQ=255)
                                                                                        ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  2:3389448/1‑100 (MQ=255)
                                                                                        ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  1:979611/1‑100 (MQ=255)
                                                                                          ggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg          <  1:318822/100‑1 (MQ=255)
                                                                                            aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCTGACGTTACACCGTGAATCTATTTACCGACAGCACACTacg          >  1:789882/1‑98 (MQ=255)
                                                                                             gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       <  2:506599/100‑1 (MQ=255)
                                                                                              tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcgc         <  1:195817/97‑1 (MQ=255)
                                                                                               aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCCCACTGCgcagca     >  2:966601/1‑100 (MQ=255)
                                                                                               aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca     >  1:470943/1‑100 (MQ=255)
                                                                                               aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca     >  1:1019842/1‑100 (MQ=255)
                                                                                                tCAGCTACTCCCCGGTTTATTGATTGCTTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt    <  1:1328354/100‑1 (MQ=255)
                                                                                                tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt    <  1:2541388/100‑1 (MQ=255)
                                                                                                 cAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACAGTGACTCTATCTACAGCCAGCACACTGCGCAGCATa   >  1:944804/1‑100 (MQ=255)
                                                                                                 cAGCTACTCCCCGGTTTATTGATTGCCTACGTAATTGCGTCTGTAATTTCGCGGACGTTACACCGTGACTCTATATACCGCCAGCACACTGCgcagca     >  2:2433042/1‑98 (MQ=255)
                                                                                                  aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg  <  2:3016201/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666370‑1666567

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