Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,085 Δ1 bp coding (2069/2664 nt) aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,0850T.97.4% 904.0 / 15.1 192coding (2069/2664 nt)aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base . (99/88);  total (102/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 192 total aligned reads displayed.

TTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAA  >  NC_000913/124990‑125179
                                                                                               |                                                                                              
   cGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGc                                                                                        >  1:2323388/1‑100 (MQ=255)
      cccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGa                                                                                     <  1:4267760/100‑1 (MQ=255)
       ccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGaa                                                                                    >  1:2018063/1‑100 (MQ=255)
        cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAAc                                                                                   >  1:7650/1‑100 (MQ=255)
             tacgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGaaaa                                                                              >  1:1738509/1‑100 (MQ=255)
              acgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAAc                                                                             <  1:2668071/100‑1 (MQ=255)
               cgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCt                                                                                       >  2:5655554/1‑89 (MQ=255)
                gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTg                                                                                      <  2:117128/89‑1 (MQ=255)
                gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTa                                                                           <  1:1374275/100‑1 (MQ=255)
                gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTa                                                                           <  1:1042464/100‑1 (MQ=255)
                  ttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTAcc                                                                         <  1:5472554/100‑1 (MQ=255)
                   tacgAAGTTGCTGTCATCATGCCTGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAAc                                                                                   <  2:1873004/89‑1 (MQ=255)
                   tacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAAc                                                                                   >  2:3454204/1‑89 (MQ=255)
                   tacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACca                                                                        <  1:3225749/100‑1 (MQ=255)
                    acgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACca                                                                        >  1:3715744/1‑100 (MQ=255)
                    acgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACg                                                                                  >  2:5573298/1‑89 (MQ=255)
                     cgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGa                                                                                 >  2:2056741/1‑89 (MQ=255)
                     cgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACcaca                                                                      >  1:212000/1‑100 (MQ=255)
                      gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGaa                                                                                <  2:1974369/89‑1 (MQ=255)
                       aaGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGaaa                                                                               <  2:3045902/89‑1 (MQ=255)
                       aaGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATg                                                                    >  1:4518623/1‑100 (MQ=255)
                        aGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGc                                                                   <  1:3132217/100‑1 (MQ=255)
                          ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCg                                                                 <  1:5125664/100‑1 (MQ=255)
                           tGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCgg                                                                <  1:1830353/100‑1 (MQ=255)
                              tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACca                                                                        <  2:2348720/89‑1 (MQ=255)
                              tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCaa                                                             >  1:3389776/1‑100 (MQ=255)
                               gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACcac                                                                       <  2:2531929/89‑1 (MQ=255)
                                tcatcaTGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACcaca                                                                      >  2:2091580/1‑89 (MQ=255)
                                tcatcaTGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACcaca                                                                      <  2:3983108/89‑1 (MQ=255)
                                tcatcaTGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACcaca                                                                      <  2:4753051/89‑1 (MQ=255)
                                 catcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACAt                                                                     >  2:999082/1‑89 (MQ=255)
                                   tcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCg                                                        >  1:1220633/1‑100 (MQ=255)
                                    catgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCgg                                                       >  1:4643174/1‑100 (MQ=255)
                                      tgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCgg                                                                >  2:3288525/1‑89 (MQ=255)
                                      tgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCgg                                                                >  2:5284356/1‑89 (MQ=255)
                                       gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGc                                                               >  2:1246190/1‑89 (MQ=255)
                                       gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                    <  1:927040/100‑1 (MQ=255)
                                       gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                    >  1:2200052/1‑100 (MQ=255)
                                       gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                    >  1:297262/1‑100 (MQ=255)
                                        catgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAgg                                                   >  1:5061554/1‑100 (MQ=255)
                                        catgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAgg                                                   <  1:1001376/100‑1 (MQ=255)
                                         atgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCaa                                                             >  2:1193606/1‑89 (MQ=255)
                                         atgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGt                                                  >  1:3474962/1‑100 (MQ=255)
                                          tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTg                                                 <  1:5593901/100‑1 (MQ=255)
                                           gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATg                                                           >  2:958489/1‑89 (MQ=255)
                                            aCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGc                                                <  1:3541308/100‑1 (MQ=255)
                                              ggTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCg                                                        >  2:630530/1‑89 (MQ=255)
                                              ggTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCg                                                        >  2:742330/1‑89 (MQ=255)
                                              ggTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGa                                             >  1:4491950/1‑100 (MQ=255)
                                               gTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAg                                            >  1:253268/1‑100 (MQ=255)
                                               gTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAg                                            <  1:5303037/100‑1 (MQ=255)
                                               gTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAg                                            >  1:4184745/1‑100 (MQ=255)
                                                tCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAATGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAg                                            >  1:4033436/1‑100 (MQ=255)
                                                  tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                    >  2:1696210/1‑89 (MQ=255)
                                                  tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGGGCTGAGGaa                                         <  1:4764482/100‑1 (MQ=255)
                                                   ggAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAgg                                                   >  2:3396032/1‑89 (MQ=255)
                                                   ggAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAgg                                                   >  2:3396271/1‑89 (MQ=255)
                                                   ggAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                        <  1:2350698/100‑1 (MQ=255)
                                                    gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAgg                                       <  1:4541420/100‑1 (MQ=255)
                                                       cGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCt                                               <  2:374613/89‑1 (MQ=255)
                                                        gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTg                                              <  2:1753304/89‑1 (MQ=255)
                                                        gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTg                                              <  2:3674361/89‑1 (MQ=255)
                                                        gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATc                                   <  1:4130867/100‑1 (MQ=255)
                                                          aTGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCg                                 <  1:4437867/100‑1 (MQ=255)
                                                           tGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAgg                                           >  2:3876040/1‑89 (MQ=255)
                                                           tGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAgg                                           >  2:420895/1‑89 (MQ=255)
                                                           tGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGt                                >  1:4688112/1‑100 (MQ=255)
                                                             tACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaa                              >  1:1672607/1‑100 (MQ=255)
                                                              aCGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                        >  2:3763137/1‑89 (MQ=255)
                                                               cGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAgg                                       <  2:1549766/89‑1 (MQ=255)
                                                               cGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAgg                                       <  2:5239922/89‑1 (MQ=255)
                                                               cGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                            >  1:3340878/1‑100 (MQ=255)
                                                                ggTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAgg                           <  1:1026431/100‑1 (MQ=255)
                                                                ggTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAgg                           <  1:5673975/100‑1 (MQ=255)
                                                                   gAAAAACAAGAGAACGTTTACTACTACACACTATGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATc                                   >  2:4355520/1‑89 (MQ=255)
                                                                   gAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTAt                        >  1:2257498/1‑100 (MQ=255)
                                                                     aaaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCg                                 >  2:4937328/1‑89 (MQ=255)
                                                                     aaaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCg                                 >  2:1825547/1‑89 (MQ=255)
                                                                      aaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGt                                <  2:2915470/89‑1 (MQ=255)
                                                                      aaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGt                                >  2:256535/1‑89 (MQ=255)
                                                                      aaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTa                     >  1:1577193/1‑100 (MQ=255)
                                                                       aaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAAGTACCACATGCCGGCAATGCCGGAAGGGGCTGAGGAAGCTGTCCGTa                               <  2:5368244/89‑1 (MQ=255)
                                                                       aaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTa                               <  2:4603729/89‑1 (MQ=255)
                                                                       aaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTa                               >  2:4142334/1‑89 (MQ=255)
                                                                         cAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                             >  2:2876626/1‑89 (MQ=255)
                                                                         cAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaa                  >  1:1094992/1‑100 (MQ=255)
                                                                          aaGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                            <  2:4813177/89‑1 (MQ=255)
                                                                          aaGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                            <  2:2687073/89‑1 (MQ=255)
                                                                          aaGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                 >  1:3904230/1‑100 (MQ=255)
                                                                           agagAACGTTTACTACTACACACTACGCTGAACGAAAACTACGACATGTCGGCAATGCCGGAAGGTGCCTAGGAAGGTATCTGGAAAGGTATCTatcacg                >  1:2845211/1‑95 (MQ=255)
                                                                           agagAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAgg                           >  2:3179273/1‑89 (MQ=255)
                                                                                aCGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaa           <  1:939071/100‑1 (MQ=255)
                                                                                 cGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaaa          <  1:3559135/100‑1 (MQ=255)
                                                                                    ttACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaa                  <  2:623480/89‑1 (MQ=255)
                                                                                       ctactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACt               >  2:2933972/1‑89 (MQ=255)
                                                                                       ctactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACt               <  2:1019498/89‑1 (MQ=255)
                                                                                       ctactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACt               >  2:2304802/1‑89 (MQ=255)
                                                                                       ctactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg    >  1:1166138/1‑100 (MQ=255)
                                                                                        tactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGa   <  1:3501659/100‑1 (MQ=255)
                                                                                         actacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGaa  >  1:2862279/1‑100 (MQ=255)
                                                                                               |                                                                                              
TTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAA  >  NC_000913/124990‑125179

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