Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,468 A→T 18.7% 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.7% 368.2 / 74.7 161Q385L (CAG→CTG) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base A (59/71);  new base T (15/15);  total (75/86)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 162 total aligned reads displayed.

TTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAG  >  NC_000913/1666370‑1666567
                                                                                                  |                                                                                                   
tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                    <  2:3479975/100‑1 (MQ=255)
tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                    <  1:2314197/100‑1 (MQ=255)
 ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                   <  2:1277763/100‑1 (MQ=255)
 ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                   <  2:2490659/100‑1 (MQ=255)
  tACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCt                                                                                                  <  2:2135175/100‑1 (MQ=255)
    cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGCGTATCAGCTAc                                                                                                >  2:1441827/1‑100 (MQ=255)
    cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGATTATCAGCTAc                                                                                                >  1:2881880/1‑100 (MQ=255)
      cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACt                                                                                               <  1:1174564/99‑1 (MQ=255)
       ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                             >  1:3721765/1‑100 (MQ=255)
         aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                           <  1:3904636/100‑1 (MQ=255)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                         <  2:373770/100‑1 (MQ=255)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:3114987/100‑1 (MQ=39)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:3656435/100‑1 (MQ=39)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:830051/100‑1 (MQ=39)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:4381550/100‑1 (MQ=255)
           tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:1658788/100‑1 (MQ=39)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                         >  1:2290385/1‑99 (MQ=255)
            gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                         <  2:3999520/99‑1 (MQ=39)
             aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGACATGACCGGGGAGTATCAGCTACTCCCCg                                                                                          >  1:4077490/1‑97 (MQ=255)
             aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                       <  2:516661/100‑1 (MQ=255)
              ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGGAGATGATGATATGTGAGATGACCGGGGGGAATCAGCTACTCCCCGGgtt                                                                                      >  1:2969928/1‑100 (MQ=255)
              ccGGGATGGCGACACTGTTGGCGGCGACCACACACGCGCCGATTATGTCGACGTTTATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtc                                                                                       <  1:1347488/99‑2 (MQ=38)
                gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTat                                                                                    >  2:1904952/1‑100 (MQ=255)
                 ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                   <  2:3414531/100‑1 (MQ=255)
                 ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                   <  2:1296003/100‑1 (MQ=255)
                  gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                  <  2:4339634/100‑1 (MQ=255)
                  gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                  <  2:1600393/100‑1 (MQ=255)
                     ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                               <  1:3026609/100‑1 (MQ=255)
                        gccacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGCATCAGCTACTCCGCGGTTTATTGATTGc                                                                             >  1:2543834/3‑99 (MQ=255)
                            cTGTTGGCGGCGACCACGCCCGCGCCGTTTATGTCGACTTTGATGATATGTGAAATGACCGGGGAGTATCAGCTATTCCCCGGTTTATTGATTGCCTGCg                                                                        <  1:2863512/100‑1 (MQ=255)
                            cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                        >  2:3907289/1‑100 (MQ=255)
                             tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGt                                                                       <  1:189298/100‑1 (MQ=255)
                              ggtGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                      >  2:2962568/3‑100 (MQ=255)
                              gTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                      <  2:752367/100‑1 (MQ=255)
                                tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                    >  1:3270249/1‑100 (MQ=255)
                                 ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTaa                                                                     >  1:3555127/1‑98 (MQ=255)
                                 ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                   >  2:2236396/1‑100 (MQ=255)
                                  gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                  >  2:2990953/1‑100 (MQ=255)
                                  gcggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                  >  2:4131384/1‑100 (MQ=255)
                                    ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCg                                                                >  2:2361494/1‑100 (MQ=255)
                                    ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCg                                                                <  2:1254177/100‑1 (MQ=255)
                                     gcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCGACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                  <  1:1444067/97‑1 (MQ=255)
                                       gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCg                                                             <  2:2649189/100‑1 (MQ=255)
                                       gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCg                                                             >  2:1558957/1‑100 (MQ=255)
                                        aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                            <  1:1904946/100‑1 (MQ=255)
                                         ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                           <  2:368654/100‑1 (MQ=255)
                                         ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                           <  2:1244411/100‑1 (MQ=255)
                                         ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                           >  2:2873426/1‑100 (MQ=255)
                                          cacgcacgCGCCGATTATGTCGACGTTGATGATCTGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTGTATTGATTGCCTGCGTAATTGAGTCGGTa                                                          >  1:2649881/1‑100 (MQ=255)
                                                 gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                   >  2:1374692/1‑100 (MQ=255)
                                                  cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCACGGTTTATTGATTGCCTGCGTAATTGCGGCGGTAATTTcgcg                                                  >  2:1081267/1‑100 (MQ=255)
                                                   gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  2:1160497/100‑1 (MQ=255)
                                                   gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                 <  2:4366550/100‑1 (MQ=255)
                                                    ccGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                                >  1:3060089/1‑100 (MQ=255)
                                                     cGATTATGTTGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                               <  1:3829940/100‑1 (MQ=255)
                                                     cGATTATGTCGACGTTGATGATATGTGGAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAc                                               <  1:956660/100‑1 (MQ=255)
                                                      gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                              >  1:3993827/1‑100 (MQ=255)
                                                        ttATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGtt                                            <  1:816890/100‑1 (MQ=255)
                                                          aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                          <  1:2509070/100‑1 (MQ=255)
                                                           tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                         <  1:53565/100‑1 (MQ=255)
                                                              cGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                      <  1:2388221/100‑1 (MQ=255)
                                                               gACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                     >  1:3611474/1‑100 (MQ=255)
                                                                aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                    <  2:3770984/100‑1 (MQ=255)
                                                                   ttGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                  <  1:3042187/99‑1 (MQ=255)
                                                                    tgatgaTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActc                                <  2:4085946/100‑1 (MQ=255)
                                                                        gatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                            >  2:1862308/1‑100 (MQ=255)
                                                                         atatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                           <  2:4162573/100‑1 (MQ=255)
                                                                         atatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                           <  2:2301517/100‑1 (MQ=255)
                                                                          tatGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatcta                          >  2:4751/1‑100 (MQ=255)
                                                                            tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  2:4222251/100‑1 (MQ=255)
                                                                            tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                        <  2:2350863/100‑1 (MQ=255)
                                                                             gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                       >  2:4312878/1‑100 (MQ=255)
                                                                              tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                      >  2:3695012/1‑100 (MQ=255)
                                                                               gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     <  1:3581193/100‑1 (MQ=255)
                                                                               gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     <  1:2666057/100‑1 (MQ=255)
                                                                               gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  1:4072741/1‑100 (MQ=255)
                                                                               gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCCTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  1:2810222/1‑100 (MQ=255)
                                                                               caataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACGCTATcgcccgcc                     >  1:1730816/7‑100 (MQ=25)
                                                                                aataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                    >  1:2328859/6‑100 (MQ=38)
                                                                                aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                     >  1:3124811/1‑99 (MQ=255)
                                                                                aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                    <  1:4243332/100‑1 (MQ=255)
                                                                                 ataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                   >  2:2219308/5‑100 (MQ=38)
                                                                                     aCCGGGGAGTAGCTGATACGCCCCGGTGTATCGATTGCCTGCGTAATTGCGTCGGTCATCTCGCCGACGTCACACCGTGACCGTATCGAACGCCAGcac                >  1:3461322/1‑99 (MQ=255)
                                                                                      ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:1032035/1‑100 (MQ=255)
                                                                                      ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:2447190/1‑100 (MQ=255)
                                                                                      ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:3794172/1‑100 (MQ=39)
                                                                                      ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  2:1418857/1‑100 (MQ=39)
                                                                                      ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  2:1398040/1‑100 (MQ=39)
                                                                                       cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac              >  1:2206997/1‑99 (MQ=39)
                                                                                        ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTCCACCGTGACTCTATCTACCGCCAGCACACt             >  1:3932596/1‑99 (MQ=255)
                                                                                        ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  2:2683951/1‑100 (MQ=255)
                                                                                        ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg            >  2:1558033/1‑100 (MQ=255)
                                                                                             gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       >  2:2863548/1‑100 (MQ=255)
                                                                                             gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       >  2:2648225/1‑100 (MQ=255)
                                                                                             gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag       <  2:860983/100‑1 (MQ=255)
                                                                                               aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGCAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca     >  1:1123555/1‑100 (MQ=255)
                                                                                               aTCAGCTACGCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCCCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca     >  1:1739011/1‑100 (MQ=255)
                                                                                                tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt    >  2:9340/1‑100 (MQ=255)
                                                                                                 cAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATa   <  2:650968/100‑1 (MQ=255)
                                                                                                  aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAg  <  1:2369262/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: