Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,470 C→A 18.0% 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,4700CA18.0% 376.9 / 66.2 150L386I (CTA→ATA) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (59/64);  new base A (13/14);  total (72/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.68e-01

Only 100 of 152 total aligned reads displayed.

TTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666371‑1666567
                                                                                                   |                                                                                                 
 tACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCt                                                                                                  <  2:937744/100‑1 (MQ=255)
  aCTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                 >  1:2238625/1‑100 (MQ=255)
  aCTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                 >  1:1133698/1‑100 (MQ=255)
   cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  1:2756691/1‑100 (MQ=255)
   cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  2:1935057/1‑100 (MQ=255)
   cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                >  2:33502/1‑100 (MQ=255)
     cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                              >  2:2325457/1‑100 (MQ=255)
       gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTccc                                                                                            >  2:2480225/1‑100 (MQ=255)
        aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                           <  1:283790/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                         <  2:2654467/100‑1 (MQ=255)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:2846792/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:2091706/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:1849460/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:1204336/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:1104674/100‑1 (MQ=39)
          tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  1:796684/100‑1 (MQ=39)
           gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                        >  1:462168/1‑100 (MQ=255)
           gACCGGGATGGCGACACTGTTAGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGCCCGGGGAGTAGCTGATACTCCCCGGt                                                                                        <  2:2358916/100‑1 (MQ=38)
             ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                      >  1:3198078/1‑100 (MQ=255)
              cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTccccagatcg                                                                                     >  2:431833/1‑94 (MQ=255)
                ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                   <  2:2845776/100‑1 (MQ=255)
                 gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                  <  1:1071102/100‑1 (MQ=255)
                 gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                  <  1:1542482/100‑1 (MQ=255)
                     gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                              <  1:2532310/100‑1 (MQ=255)
                     gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                              <  2:1889671/100‑1 (MQ=255)
                      cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                             >  2:3031902/1‑100 (MQ=255)
                      cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTACTGATTGc                                                                             >  1:1753467/1‑100 (MQ=255)
                      cGACACTGTTAGCGGCGACCACGAACGCGCCGATTATCTCGACGTTGATGATATGGGAAATGACCGGGGAGTGTCAGCGACTCCCCGGTTGATTGATTGc                                                                             >  2:3539441/1‑100 (MQ=255)
                           cTGTTGGCGGCGACCACGCACGCGCCGAGTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCTCCGGTTTATTGATTGCCTGCg                                                                        <  2:1679164/100‑1 (MQ=255)
                           cTGTTGGCGGCGACCACGCACGCGCAGATTATGTAGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                        >  2:2505343/1‑100 (MQ=255)
                               tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                    <  1:1418384/100‑1 (MQ=255)
                                ggcggcGACCACGCACGCGCCGAATATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                   >  1:2004592/1‑100 (MQ=255)
                                   ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGc                                                                 >  2:3580342/1‑99 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              <  1:2595991/100‑1 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              <  2:3026560/100‑1 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              >  1:3577275/1‑100 (MQ=255)
                                     cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                              <  2:2480548/100‑1 (MQ=255)
                                         cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCAGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                          <  2:142403/100‑1 (MQ=255)
                                          acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTaa                                                         >  2:2773346/1‑100 (MQ=255)
                                              acgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTc                                                     <  1:23252/100‑1 (MQ=255)
                                                ctgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                   <  1:1093088/98‑1 (MQ=255)
                                                 cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                  >  2:723477/1‑100 (MQ=255)
                                                  gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcag                                                 >  2:1093080/1‑98 (MQ=255)
                                                  gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  2:1170345/100‑1 (MQ=255)
                                                  gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  1:1777926/100‑1 (MQ=255)
                                                   ccGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                                <  2:3173346/100‑1 (MQ=255)
                                                     gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                              <  1:3513018/100‑1 (MQ=255)
                                                       ttATGTCGACGTTGATGATATGTAAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCCGACg                                              <  2:350917/98‑1 (MQ=255)
                                                         aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                          <  1:616170/100‑1 (MQ=255)
                                                         aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                          <  1:723480/100‑1 (MQ=255)
                                                          tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                         >  1:123174/1‑100 (MQ=255)
                                                          tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                         <  1:2265576/100‑1 (MQ=255)
                                                           ntCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTacac                                        >  2:1883868/2‑100 (MQ=255)
                                                            tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAcc                                       <  1:3611248/100‑1 (MQ=255)
                                                             cGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                      <  1:393194/100‑1 (MQ=255)
                                                               aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATGGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                    <  2:2530695/100‑1 (MQ=255)
                                                                 gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                  <  2:2135159/100‑1 (MQ=255)
                                                                  ttGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                   >  2:1631716/1‑98 (MQ=255)
                                                                     atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                              <  1:3100776/100‑1 (MQ=255)
                                                                      tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                             >  1:922013/1‑100 (MQ=255)
                                                                      tcatgtgttAAATGACCGGGGAGTATCAGCTACTCCCCGGCTTGTTGATTGCCTGCGTGATTGCGGCGGGGATTGCGCGGACGTTACACCCTGAgtc                                >  2:2428561/6‑97 (MQ=255)
                                                                         naggtgAAATGACCGGGGAGTATCAGCGACTCCCCGGGTTATGGAGTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGGTACACCGTGACtctacct                           >  2:729543/4‑99 (MQ=255)
                                                                           tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  1:1938293/100‑1 (MQ=255)
                                                                           tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  2:2072899/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                      >  2:2009412/1‑100 (MQ=255)
                                                                             tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                      <  1:3027030/100‑1 (MQ=255)
                                                                              gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  2:3480029/1‑100 (MQ=255)
                                                                              gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     <  1:998029/100‑1 (MQ=255)
                                                                              gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTAATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTAGCTAACGcc                     >  2:1813023/1‑100 (MQ=255)
                                                                              caataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  2:1033629/7‑100 (MQ=37)
                                                                               aataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                    >  1:330804/6‑100 (MQ=37)
                                                                                ataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                   >  2:2517539/5‑100 (MQ=38)
                                                                                  aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                 >  2:14476/3‑100 (MQ=38)
                                                                                    aCCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcaca               >  1:2163667/1‑100 (MQ=255)
                                                                                     ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              <  1:1659659/100‑1 (MQ=255)
                                                                                     ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:3136594/1‑100 (MQ=255)
                                                                                     ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              <  1:3528976/100‑1 (MQ=255)
                                                                                     ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:989926/1‑100 (MQ=255)
                                                                                     ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACAGTGACTCTATCTACCGCCAGCacac              >  1:2674026/1‑100 (MQ=255)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTCATTGCGTCGGTAATTTCGCGGACGATCCACCGTGACTCTATCTACCGCCAGCacac              >  2:2544345/1‑100 (MQ=38)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:1057699/1‑100 (MQ=39)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:473032/1‑100 (MQ=255)
                                                                                     ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:2821277/1‑100 (MQ=39)
                                                                                      cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt             >  1:2283682/1‑100 (MQ=255)
                                                                                       ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  1:139804/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)
                                                                                            gtctCAGCTACTCCCCGGTTTATTGATTGCCCGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTGTATCTACCGCCAGCACACTGCgcan       <  1:729546/97‑2 (MQ=255)
                                                                                            gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       >  1:2263831/1‑100 (MQ=255)
                                                                                            gTATCAGCTACTCCCCGGTTTATTGATTGCCGGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       <  2:468619/100‑1 (MQ=255)
                                                                                             tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagc      <  2:391870/100‑1 (MQ=255)
                                                                                              aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca     <  1:1485913/100‑1 (MQ=255)
                                                                                               tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt    >  1:3123901/1‑100 (MQ=255)
                                                                                                 aacTACTCCCCGGATTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg  >  1:2454405/3‑100 (MQ=255)
                                                                                                 aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTACTTTCGCGGACGTTACACCGTGACGCTATCTACCGCCGGCACACTGCGCAGCAtcg  >  2:2225374/1‑98 (MQ=255)
                                                                                                 aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCCGCACACTCCGCAGCATAg  >  2:1777807/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666371‑1666567

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