Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,470 C→A 17.2% L386I (CTA→ATA)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4700CA17.2% 368.8 / 54.3 146L386I (CTA→ATA) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (55/64);  new base A (12/14);  total (68/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 149 total aligned reads displayed.

TTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGC  >  NC_000913/1666371‑1666568
                                                                                                   |                                                                                                  
ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                    <  2:3018923/100‑1 (MQ=255)
ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                    <  1:2617158/100‑1 (MQ=255)
   cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                 >  1:1338809/1‑100 (MQ=255)
   cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                 >  1:710693/1‑100 (MQ=255)
     cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                               >  1:2238937/1‑100 (MQ=255)
      ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                              <  1:3397676/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)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                          <  2:1223809/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATAAGCTACTCCCCgg                                                                                          >  1:870836/1‑100 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  1:49475/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  2:2638137/100‑1 (MQ=39)
           gACCGGGATGGTGACACTGTTGGCGGCGACCACGTACGTGCCGATTATGTCGATGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                         <  1:221435/100‑1 (MQ=37)
           gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGCTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                         <  2:2503707/100‑1 (MQ=38)
           gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                         >  2:1250281/1‑100 (MQ=255)
           gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                         <  2:1675572/100‑1 (MQ=255)
           cttccgatctgCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                         <  1:3737986/90‑1 (MQ=37)
             ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                          <  2:166459/97‑1 (MQ=39)
             ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtca                                                                                       <  2:1468730/100‑3 (MQ=38)
              cGGTATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATAGGTGAAATCACTGGTGTGTAGCTGATACTCCCCGGtcat                                                                                      <  2:2418281/100‑4 (MQ=255)
                ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                    >  2:1478528/1‑100 (MQ=255)
                  aTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattga                                                                                  <  1:2469197/100‑1 (MQ=255)
                   tGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                                 <  2:424/100‑1 (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)
                     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)
                         cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                           >  1:2429656/1‑100 (MQ=255)
                         cacTGTTGGCGGCGACCACGCACGAGCAAATTATGTCGACGTAGATGATATGTGAAATTACCGGGGAGTATCAGTTACTCCCCAGTTTattgaatgac                                                                             >  2:3750064/1‑96 (MQ=255)
                           cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                         <  1:397546/100‑1 (MQ=255)
                            tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGt                                                                        >  1:2898270/1‑100 (MQ=255)
                             gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAAAAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                       >  1:17701/1‑100 (MQ=255)
                                ggcggcGACCACGCACGCGCAGATTATGTCGACGTTGATGATATTTGAAATTACCGGGGAGTATCAGCTACTCCCTGGTTTATTGATTGCCTGCGTAAt                                                                     <  2:349708/99‑1 (MQ=255)
                                 gcggctACCACGCACGCGCCTATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTTCCTGCGTAATTg                                                                   <  1:1139798/100‑1 (MQ=255)
                                   ggcgACCACGCACGCGCTGTTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCg                                                                 <  1:2471787/100‑1 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                               <  1:2582791/100‑1 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                               <  1:3427785/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)
                                        ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                            >  2:1308976/1‑100 (MQ=255)
                                               cgcgCCGATTATGTTGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                     <  1:3302282/100‑1 (MQ=255)
                                                 cgcCGATTATGTGGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                   <  1:3624670/100‑1 (MQ=255)
                                                 cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                   <  1:181846/100‑1 (MQ=255)
                                                 cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                   <  1:3871530/100‑1 (MQ=255)
                                                    cGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                                >  2:3847950/1‑100 (MQ=255)
                                                    cGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                                >  2:2917775/1‑100 (MQ=255)
                                                      aTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGt                                              <  1:2426299/100‑1 (MQ=255)
                                                         aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                           >  1:146019/1‑100 (MQ=255)
                                                         aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                           <  2:2426502/100‑1 (MQ=255)
                                                         aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                           <  2:2121052/100‑1 (MQ=255)
                                                          tGTCGACGTTGATGATATTTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                          <  1:1250295/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:1381155/100‑1 (MQ=255)
                                                                 gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                   <  1:283735/100‑1 (MQ=255)
                                                                    gatgatATGTGAAATGACCGGGGAGTATCAGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAAAGTGActct                                >  2:195812/1‑100 (MQ=255)
                                                                     atgatATGTGAAATGACCGTGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                               <  1:2061134/100‑1 (MQ=255)
                                                                     atgatATGTGAAATGACCGGGGAGTATAAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                               <  2:3576951/100‑1 (MQ=255)
                                                                          atGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAc                          >  1:2545424/1‑100 (MQ=255)
                                                                           tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                         >  2:58049/1‑100 (MQ=255)
                                                                            gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                        >  2:1177826/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)
                                                                              gAAATGACCGGGGAGTATCAGATACTCCCCGGTTTATTGATTGCCTGCTTAATTGCGTCGGTACTTTCTCGGACGTTACACCGTTACTCTATCTACCGcc                      <  2:1862553/100‑1 (MQ=255)
                                                                                aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                    <  2:1991486/100‑1 (MQ=255)
                                                                                aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                    <  2:1881727/100‑1 (MQ=255)
                                                                                 aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                   <  2:2322239/100‑1 (MQ=255)
                                                                                 aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                   <  2:1755789/100‑1 (MQ=255)
                                                                                  tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                  <  2:1735057/100‑1 (MQ=255)
                                                                                  tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                  >  1:1840789/1‑100 (MQ=255)
                                                                                   gACCGGGGAGTATCAGCTACTCCCTGGTTTATTGATTGTCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcac                 <  1:3660833/100‑1 (MQ=255)
                                                                                     ccGGGGAGTAGCTGATACTCCACAGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTAACACCGTGACTCTATCTACCGCCAGCagac               >  2:229814/1‑100 (MQ=255)
                                                                                      cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCTCCAGCACACt              >  2:2563163/1‑100 (MQ=255)
                                                                                      cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac               >  2:1385326/1‑99 (MQ=255)
                                                                                       ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg             >  1:979611/1‑100 (MQ=255)
                                                                                       ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg             >  2:3389448/1‑100 (MQ=255)
                                                                                         ggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg           <  1:318822/100‑1 (MQ=255)
                                                                                           aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCTGACGTTACACCGTGAATCTATTTACCGACAGCACACTacg           >  1:789882/1‑98 (MQ=255)
                                                                                            gcaTCAGCTACTCCCCGGTTTATTGATTGCTTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag        <  1:938545/98‑1 (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:1019842/1‑100 (MQ=255)
                                                                                              aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca      >  1:470943/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)
                                                                                                 aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg   <  2:3016201/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGC  >  NC_000913/1666371‑1666568

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: