Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 474,643 G→T 48.7% V173V (GTC→GTA tesB ← acyl‑CoA thioesterase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913474,6430GT48.7% 28.8 / 248.8 165V173V (GTC→GTAtesBacyl‑CoA thioesterase 2
Reads supporting (aligned to +/- strand):  ref base G (36/48);  new base T (38/42);  total (75/90)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 167 total aligned reads displayed.

GTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCC  >  NC_000913/474544‑474742
                                                                                                   |                                                                                                   
   aaGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACg                                                                                                  >  2:1587574/1‑100 (MQ=255)
   aaGATCAGAAGCGTAACCGAGCAGATACTGATCAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACg                                                                                                  >  1:1010432/1‑100 (MQ=255)
    aGATCAGAAGCGTAACCGAGCAGATACTGATGAACGGGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGt                                                                                                 >  2:3011628/1‑100 (MQ=255)
        cAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGAcc                                                                                             <  1:2812376/100‑1 (MQ=255)
        cAGAAGCGTAACCGAGCAGATACTGATATACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTg                                                                                                <  2:3763718/97‑1 (MQ=255)
          gAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCtt                                                                                           <  1:4221360/100‑1 (MQ=255)
          gAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCtt                                                                                           <  2:456102/100‑1 (MQ=255)
          gAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCtt                                                                                           <  1:1108508/100‑1 (MQ=255)
            aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACTCCTGACGATGTGGTTCTGCTACGTGACCTTTc                                                                                         <  2:4059023/100‑1 (MQ=255)
            aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTc                                                                                         >  2:1688371/1‑100 (MQ=255)
            aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTc                                                                                         <  1:4045282/100‑1 (MQ=255)
            aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTc                                                                                         <  2:3251209/100‑1 (MQ=255)
            aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCAGTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTc                                                                                         <  1:558103/100‑1 (MQ=255)
              cGTAACCGAGCAGATACTGATGAACGCGCAGGTTATTCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAg                                                                                       <  1:3926702/100‑1 (MQ=255)
              cGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAg                                                                                       <  2:1527518/100‑1 (MQ=255)
              cGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAg                                                                                       >  2:1223021/1‑100 (MQ=255)
               gTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGt                                                                                      <  2:1027748/100‑1 (MQ=255)
               gTAACCGAGCAGATACTCATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGt                                                                                      <  1:3527223/100‑1 (MQ=255)
                   ccGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACCCCTGACGATGTGGTTCTGCGACATGACCTTTCAGTGGGt                                                                                  >  1:1730575/1‑100 (MQ=255)
                   ccGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGt                                                                                  >  1:266089/1‑100 (MQ=255)
                   ccGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGt                                                                                  <  1:3019668/100‑1 (MQ=255)
                     gAGCAGATACTGATGAACGCGCAGGTCATCCGGCAGGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTggg                                                                                   <  1:824069/97‑1 (MQ=255)
                        cAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGa                                                                             <  2:2566751/100‑1 (MQ=255)
                        cAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGa                                                                             <  2:1813604/100‑1 (MQ=255)
                         aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACTATTTGCGCGGATCCACAGCTGACGATGTGGTTCTTCTACGTTACCTTTCAGTGGGTTATGa                                                                             <  1:2649142/99‑1 (MQ=255)
                         aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGaa                                                                            <  2:3093851/100‑1 (MQ=255)
                          gATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGaaa                                                                           >  1:3290856/1‑100 (MQ=255)
                          gATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGaa                                                                            >  2:1295162/1‑99 (MQ=255)
                           aTACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAAc                                                                          >  1:3506185/1‑100 (MQ=255)
                             aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTc                                                                        >  1:1335856/1‑100 (MQ=255)
                              cTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTcc                                                                       >  2:648614/1‑100 (MQ=255)
                              cTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTcc                                                                       <  1:3900534/100‑1 (MQ=255)
                                gatgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCAc                                                                     <  2:3166847/100‑1 (MQ=255)
                                gatgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCAc                                                                     >  2:186836/1‑100 (MQ=255)
                                  tgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcg                                                                   >  1:2788843/1‑100 (MQ=255)
                                  tgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcg                                                                   <  2:3375499/100‑1 (MQ=255)
                                   gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcgg                                                                  >  2:841265/1‑100 (MQ=255)
                                    aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCCCcgga                                                                 >  1:2077403/1‑100 (MQ=255)
                                      cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcggacg                                                               <  1:3872527/100‑1 (MQ=255)
                                       gcgcAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcggacgg                                                              <  2:1742987/100‑1 (MQ=255)
                                            ggTCATCCTGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTc                                                         >  1:426107/1‑100 (MQ=255)
                                            ggTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTc                                                         >  1:1479131/1‑100 (MQ=255)
                                             gTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTcc                                                        <  1:1300443/100‑1 (MQ=255)
                                              tCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCa                                                       >  1:3719396/1‑100 (MQ=255)
                                              tCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCa                                                       <  2:343171/100‑1 (MQ=255)
                                                aTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGc                                                     <  1:3220367/100‑1 (MQ=255)
                                                 tCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCg                                                    >  2:3677370/1‑100 (MQ=255)
                                                 tCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCg                                                    <  1:3041236/100‑1 (MQ=255)
                                                  ccGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCgg                                                   >  1:3359032/1‑100 (MQ=255)
                                                  ccGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCgg                                                   >  2:3000084/1‑100 (MQ=255)
                                                    ggCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGAc                                                 <  1:695830/100‑1 (MQ=255)
                                                    ggCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGAc                                                 <  2:231181/100‑1 (MQ=255)
                                                     gCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACg                                                <  2:2007928/100‑1 (MQ=255)
                                                     gCACGCTACCAGTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGAc                                                 <  1:973796/99‑1 (MQ=255)
                                                       aCGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGAt                                              <  2:2074428/100‑1 (MQ=255)
                                                       aCGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGAt                                              <  2:3850129/100‑1 (MQ=255)
                                                       aCGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGAt                                              <  2:3588883/100‑1 (MQ=255)
                                                         gCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCg                                            <  2:426109/100‑1 (MQ=255)
                                                          cTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGc                                           >  2:826754/1‑100 (MQ=255)
                                                           tagcATCTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCa                                          <  1:1680040/97‑1 (MQ=255)
                                                           tACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTGTCAGTGGGTTATGAAGCTCCACCGCACGGACTTCCAGCGGACGATCGCa                                          <  1:2251934/100‑1 (MQ=255)
                                                              cATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGAt                                       >  1:3653688/1‑100 (MQ=255)
                                                               aTTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATg                                      <  2:723546/100‑1 (MQ=255)
                                                                 ttGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGaa                                    <  1:3715390/100‑1 (MQ=255)
                                                                 ttGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGaa                                    <  2:3468052/100‑1 (MQ=255)
                                                                 ttGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGaa                                    >  2:973240/1‑100 (MQ=255)
                                                                  tGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAt                                   <  2:3666843/100‑1 (MQ=255)
                                                                  tGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAt                                   >  1:2456023/1‑100 (MQ=255)
                                                                  tGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTCCAGTGGGTTCTGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGa                                     >  1:406743/1‑98 (MQ=255)
                                                                   gcgcGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAtt                                  <  2:2805883/100‑1 (MQ=255)
                                                                   gcgcGGATCCACACCGGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAtt                                  <  2:206324/100‑1 (MQ=255)
                                                                    cgcgGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGCCGATCGCCGATGAAttt                                 >  1:1473258/1‑100 (MQ=255)
                                                                    cgcgGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAttt                                 <  2:2196283/100‑1 (MQ=255)
                                                                     ggtgagccACACCTGGCGATGTGGTTCTGCGAGGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTa                                <  1:3949829/94‑1 (MQ=255)
                                                                       ggATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATc                              <  1:1367804/100‑1 (MQ=255)
                                                                        gATCCACACCTGACGATGTGGTTCTGCGACGTGACTCTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCt                             <  2:1062984/100‑1 (MQ=255)
                                                                         aTCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCtt                            <  2:1556775/100‑1 (MQ=255)
                                                                         aTCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCtt                            <  2:681484/100‑1 (MQ=255)
                                                                           ccACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTc                          <  1:2023681/100‑1 (MQ=255)
                                                                            cacaCCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCa                         <  2:2552245/100‑1 (MQ=255)
                                                                            cacaCCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCa                         >  2:595857/1‑100 (MQ=255)
                                                                            cacaCCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCa                         >  2:875782/1‑100 (MQ=255)
                                                                             acacCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCa                         <  1:1244225/99‑1 (MQ=255)
                                                                             acacCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAg                        <  1:1789982/100‑1 (MQ=255)
                                                                              cacCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGc                       <  1:2343142/100‑1 (MQ=255)
                                                                                ccTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCAc                     >  2:1942015/1‑100 (MQ=255)
                                                                                 cTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACt                    >  2:2298441/1‑100 (MQ=255)
                                                                                  tGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTg                   >  1:2086586/1‑100 (MQ=255)
                                                                                     cGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcg                <  1:2531705/100‑1 (MQ=255)
                                                                                       aTGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcggc              <  1:2506870/100‑1 (MQ=255)
                                                                                          tgGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagc           <  2:1297812/100‑1 (MQ=255)
                                                                                          tgGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagc           <  1:61664/100‑1 (MQ=255)
                                                                                           ggTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagca          <  2:4162628/100‑1 (MQ=255)
                                                                                           ggTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTAGCTTTCAGCACTGGCGgcagca          >  1:1067716/1‑100 (MQ=255)
                                                                                           ggTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagca          >  1:3127601/1‑100 (MQ=255)
                                                                                             ttCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAgg        >  2:1760540/1‑100 (MQ=255)
                                                                                              tCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGt       >  1:2874679/1‑100 (MQ=255)
                                                                                               cTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGt       <  1:2682527/99‑1 (MQ=255)
                                                                                               cTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTg      >  2:2419405/1‑100 (MQ=255)
                                                                                                tGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTgc     >  1:1240723/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCC  >  NC_000913/474544‑474742

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: