Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,466 T→G 18.0% Y384* (TAT→TAG clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4660TG18.0% 337.2 / 78.2 150Y384* (TAT→TAGclcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base T (57/66);  new base G (13/14);  total (70/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

Only 100 of 152 total aligned reads displayed.

CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAT  >  NC_000913/1666368‑1666565
                                                                                                  |                                                                                                   
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                    >  1:3075837/1‑100 (MQ=255)
 ttttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCa                                                                                                   >  1:1057513/1‑100 (MQ=255)
   ttACTCGGATTGACCGGGATGGCGACGCTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAAAGACCGGGGAGTATCAga                                                                                                 <  2:122030/100‑2 (MQ=255)
   ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                 <  2:2605883/100‑1 (MQ=255)
    tACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCt                                                                                                <  2:937744/100‑1 (MQ=255)
     aCTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                               >  1:2238625/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              <  2:2801240/100‑1 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              >  1:2583673/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              >  1:2756691/1‑100 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              >  2:1935057/1‑100 (MQ=255)
        cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                            >  2:2325457/1‑100 (MQ=255)
        cGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTc                                                                                            >  2:3611221/1‑100 (MQ=255)
         ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                           <  2:250415/100‑1 (MQ=255)
          gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTccc                                                                                          >  2:2480225/1‑100 (MQ=255)
           cgatctcGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                         <  1:431835/94‑1 (MQ=255)
           aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                         <  1:283790/100‑1 (MQ=255)
            ttGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                        <  1:2483192/100‑1 (MQ=255)
             tGACCTGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                       <  1:3481406/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                       <  2:2654467/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:99022/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  2:1104674/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  2:1204336/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:2823723/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  2:2091706/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:796684/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:1849460/100‑1 (MQ=39)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:2445411/100‑1 (MQ=39)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                      <  1:2682082/100‑1 (MQ=255)
              gACCGGGATGGCGACACTGGTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                      <  2:1047948/100‑1 (MQ=255)
               aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                     <  1:1403635/100‑1 (MQ=255)
               aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGtt                                                                                     >  2:1729373/1‑100 (MQ=255)
                ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                    >  1:3198078/1‑100 (MQ=255)
                   ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                 <  2:2845776/100‑1 (MQ=255)
                   ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                 <  1:3577092/100‑1 (MQ=255)
                   ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                 <  1:2034975/100‑1 (MQ=255)
                   ggATCGCTACACTGTTGGCGGCGACCACGCACGTGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcattt                                                                                 <  2:3110683/100‑6 (MQ=25)
                    gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                <  1:1542482/100‑1 (MQ=255)
                       ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                             >  1:2271147/1‑100 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            <  2:1889671/100‑1 (MQ=255)
                         cGACACTGTTAGCGGCGACCACGAACGCGCCGATTATCTCGACGTTGATGATATGGGAAATGACCGGGGAGTGTCAGCGACTCCCCGGTTGATTGATTGc                                                                           >  2:3539441/1‑100 (MQ=255)
                           acacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCt                                                                         <  2:3461903/100‑1 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGAGTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCTCCGGTTTATTGATTGCCTGCg                                                                      <  2:1679164/100‑1 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCAGATTATGTAGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                      >  2:2505343/1‑100 (MQ=255)
                                gTTGGCGGCGACCACGCACGCGCCGCTTTTGTCGACGTTGATGAGATGTGCAATGACCGGGGCGTCTCATCTACTCGCCGGGGTATCGGTTGCCTGc                                                                       >  2:1129896/1‑97 (MQ=255)
                                 ttGGCGGCGACAACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTa                                                                    <  1:386152/99‑1 (MQ=255)
                                  tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAt                                                                  <  1:1418384/100‑1 (MQ=255)
                                   ggcggcGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                 >  1:1075787/1‑100 (MQ=255)
                                   ggcggcGACCACGCACGCGCCGAATATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAAtt                                                                 >  1:2004592/1‑100 (MQ=255)
                                       gcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGt                                                             >  1:2263116/1‑100 (MQ=255)
                                        cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                            <  1:2595991/100‑1 (MQ=255)
                                        cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                            >  1:904998/1‑100 (MQ=255)
                                        cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                            <  2:3026560/100‑1 (MQ=255)
                                            cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCAGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                        <  2:142403/100‑1 (MQ=255)
                                                 acgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTc                                                   <  1:23252/100‑1 (MQ=255)
                                                  cgcgCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                  >  1:3303684/1‑100 (MQ=255)
                                                   ctgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                 <  1:1093088/98‑1 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                <  1:628219/100‑1 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                <  1:1866110/100‑1 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                >  2:723477/1‑100 (MQ=255)
                                                     gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                               <  1:1777926/100‑1 (MQ=255)
                                                      ccGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                              <  2:3173346/100‑1 (MQ=255)
                                                          ttATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGtt                                          <  1:3135374/100‑1 (MQ=255)
                                                            aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                        <  1:616170/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)
                                                                  aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                  <  2:2524878/100‑1 (MQ=255)
                                                                  aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATGGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                  <  2:2530695/100‑1 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                   <  2:2446251/97‑1 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACAGTGAc                                >  2:2423294/1‑100 (MQ=255)
                                                                     ttGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                 >  2:1631716/1‑98 (MQ=255)
                                                                     ttGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAct                               <  1:1699968/100‑1 (MQ=255)
                                                                        atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                            <  2:462166/100‑1 (MQ=255)
                                                                            naggtgAAATGACCGGGGAGTATCAGCGACTCCCCGGGTTATGGAGTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGGTACACCGTGACtctacct                         >  2:729543/4‑99 (MQ=255)
                                                                              tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                      <  2:2072899/100‑1 (MQ=255)
                                                                                tgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                    >  2:2009412/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  2:3480029/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  1:1510903/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTAATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTAGCTAACGcc                   >  2:1813023/1‑100 (MQ=255)
                                                                                 caataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  2:1033629/7‑100 (MQ=37)
                                                                                  aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACATTACACCGTGACTCTATCTACCGCCa                  <  1:54165/100‑1 (MQ=255)
                                                                                   aaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTCATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  2:3518072/1‑98 (MQ=255)
                                                                                    ataaCCGGGGAGTATCATCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                <  2:3001909/97‑1 (MQ=255)
                                                                                    aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                <  1:3031925/100‑1 (MQ=255)
                                                                                     aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca               >  2:14476/3‑100 (MQ=38)
                                                                                      gACCGGGGAGTATCAGCTACTCCCCGCTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGGCGTTACACCGTGACTCTATCTACCGCCAGcac              >  2:2304932/1‑100 (MQ=255)
                                                                                        ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:989926/1‑100 (MQ=255)
                                                                                        ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            <  1:1659659/100‑1 (MQ=255)
                                                                                        ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            <  1:3528976/100‑1 (MQ=255)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:1057699/1‑100 (MQ=39)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:473032/1‑100 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          >  1:3337216/1‑100 (MQ=255)
                                                                                             gAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcgc       <  1:1791730/100‑1 (MQ=255)
                                                                                               gtctCAGCTACTCCCCGGTTTATTGATTGCCCGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTGTATCTACCGCCAGCACACTGCgcan     <  1:729546/97‑2 (MQ=255)
                                                                                               gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag     <  1:1176425/100‑1 (MQ=255)
                                                                                               gTATCAGCTACTCCCCGGTTTATTGATTGCCGGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag     <  2:468619/100‑1 (MQ=255)
                                                                                                tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagc    <  2:391870/100‑1 (MQ=255)
                                                                                                  tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt  >  1:3123901/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAT  >  NC_000913/1666368‑1666565

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: