Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 552,300 C→T A8A (GCG→GCA purE ← phosphoribosylaminoimidazole carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913552,3000CT85.7% 625.2 / 62.1 225A8A (GCG→GCApurEphosphoribosylaminoimidazole carboxylase
Reads supporting (aligned to +/- strand):  ref base C (7/25);  new base T (83/109);  total (91/134)
Fisher's exact test for biased strand distribution p-value = 3.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 234 total aligned reads displayed.

ACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAG  >  NC_000913/552208‑552388
                                                                                            |                                                                                        
  gtggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACgtgc                                                                                        >  1:5067123/1‑93 (MQ=255)
   tggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGcc                                                                                       <  1:5033891/93‑1 (MQ=255)
    ggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCg                                                                                      <  1:3408583/93‑1 (MQ=255)
      tgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGa                                                                                    <  1:3682392/93‑1 (MQ=255)
        cGGGACATTCAGGTTTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGTTTTTGGCCCCCACCAGGATGGCGCCACGTGCCGGatt                                                                                  <  1:6698753/93‑1 (MQ=255)
        cGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGatt                                                                                  <  2:5958308/93‑1 (MQ=255)
         gggACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGatta                                                                                 <  2:7287178/93‑1 (MQ=255)
         gggACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGatta                                                                                 >  2:5151264/1‑93 (MQ=255)
          ggACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGattat                                                                                <  1:338391/93‑1 (MQ=255)
           gACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGattatt                                                                               >  2:3011845/1‑93 (MQ=255)
            aCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCGGTCGCTTTTGGCCCCCACCCCGATGGCGACACGTGCCGGATTATTg                                                                              <  1:3641405/93‑1 (MQ=255)
            aCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTg                                                                              <  1:4298107/93‑1 (MQ=255)
            aCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTg                                                                              <  2:6843374/93‑1 (MQ=255)
             cATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGc                                                                             <  1:4906089/93‑1 (MQ=255)
              aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCg                                                                            <  2:4150743/93‑1 (MQ=255)
              aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCg                                                                            >  2:703113/1‑93 (MQ=255)
               ttCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCgg                                                                           >  2:2749400/1‑93 (MQ=255)
               ttCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCgg                                                                           >  1:1093045/1‑93 (MQ=255)
               ttCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCgg                                                                           <  2:3936502/93‑1 (MQ=255)
                  aGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGaa                                                                        <  1:2796963/93‑1 (MQ=255)
                   ggATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAg                                                                       >  1:831680/1‑93 (MQ=255)
                    gATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGa                                                                      >  2:1366244/1‑93 (MQ=255)
                    gATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGa                                                                      <  1:446021/93‑1 (MQ=255)
                      tttCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACa                                                                    <  1:370580/93‑1 (MQ=255)
                      tttCGAAGATGTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACa                                                                    <  2:4836721/93‑1 (MQ=255)
                       ttCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACAt                                                                   >  2:1474469/1‑93 (MQ=255)
                        tCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATg                                                                  >  1:460801/1‑93 (MQ=255)
                         cGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGc                                                                 <  1:4629000/93‑1 (MQ=255)
                          gAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCg                                                                >  1:4147770/1‑93 (MQ=255)
                           aaGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGt                                                               <  1:5634069/93‑1 (MQ=255)
                           aaGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGt                                                               <  2:1299956/93‑1 (MQ=255)
                            aGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTc                                                              >  2:381639/1‑93 (MQ=255)
                            aGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTc                                                              >  1:1338201/1‑93 (MQ=255)
                             gATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCt                                                             >  2:1507950/1‑93 (MQ=255)
                             gATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCt                                                             <  1:3956744/93‑1 (MQ=255)
                             gATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCt                                                             <  1:1007354/93‑1 (MQ=255)
                              aTTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCtt                                                            >  2:2200704/1‑93 (MQ=255)
                               tttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTa                                                           >  2:2256720/1‑93 (MQ=255)
                               tttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTa                                                           >  1:6402272/1‑93 (MQ=255)
                                ttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTaa                                                          >  2:5324227/1‑93 (MQ=255)
                                 tCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTaaa                                                         <  2:6212018/93‑1 (MQ=255)
                                 tCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTaaa                                                         <  1:594037/93‑1 (MQ=255)
                                  cggcggCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTaaaa                                                        <  1:6256915/93‑1 (MQ=255)
                                   ggcggcGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAAc                                                       <  2:2670074/93‑1 (MQ=255)
                                   ggcggcGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAAc                                                       <  1:6818925/93‑1 (MQ=255)
                                    gcggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACt                                                      <  1:589022/93‑1 (MQ=255)
                                     cggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTc                                                     <  1:7114605/93‑1 (MQ=255)
                                       gcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCt                                                   >  1:189189/1‑93 (MQ=255)
                                         gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCtgt                                                 <  1:1492552/93‑1 (MQ=255)
                                          aaCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCtgtg                                                >  2:6280360/1‑93 (MQ=255)
                                           aCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtgg                                               <  1:3836096/93‑1 (MQ=255)
                                            cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggt                                              >  1:7060854/1‑93 (MQ=255)
                                            cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggt                                              <  1:259537/93‑1 (MQ=255)
                                            cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggt                                              >  1:2378376/1‑93 (MQ=255)
                                             tGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggtg                                             <  2:1032122/93‑1 (MQ=255)
                                               cATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGca                                           <  1:5579365/93‑1 (MQ=255)
                                                aTGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGcac                                          >  1:1154910/1‑93 (MQ=255)
                                                  ggTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACaa                                        >  1:6230221/1‑93 (MQ=255)
                                                     aGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActc                                     <  1:321002/93‑1 (MQ=255)
                                                       cccAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACtctc                                   <  1:6206368/93‑1 (MQ=255)
                                                       cccAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACtctc                                   <  1:3435902/93‑1 (MQ=255)
                                                        ccAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCg                                  <  2:4974849/93‑1 (MQ=255)
                                                         cAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                 >  1:2286453/1‑93 (MQ=255)
                                                         cAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                 >  1:1403555/1‑93 (MQ=255)
                                                          aGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGc                                >  1:5805398/1‑93 (MQ=255)
                                                           gTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                               <  1:6377269/93‑1 (MQ=255)
                                                           gTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                               >  2:7452839/1‑93 (MQ=255)
                                                           gTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                               >  2:4287795/1‑93 (MQ=255)
                                                            tCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCtt                              <  2:2219394/93‑1 (MQ=255)
                                                            tCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCtt                              <  1:3085884/93‑1 (MQ=255)
                                                             cGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCttt                             <  1:3917760/93‑1 (MQ=255)
                                                              gCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTa                            <  1:889769/93‑1 (MQ=255)
                                                              gCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTa                            <  2:3361338/93‑1 (MQ=255)
                                                               cTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTag                           <  2:5432990/93‑1 (MQ=255)
                                                                ttttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTaga                          <  2:3847922/93‑1 (MQ=255)
                                                                ttttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTaga                          <  2:4474256/93‑1 (MQ=255)
                                                                 tttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTagag                         <  2:7349933/93‑1 (MQ=255)
                                                                  ttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                        <  1:2998774/93‑1 (MQ=255)
                                                                   tGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAggg                       <  1:2494970/93‑1 (MQ=255)
                                                                      aCCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCAGCTTTAGAGGGcac                    <  1:5290541/93‑1 (MQ=255)
                                                                        cccATCACGATGGCGATACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACag                  >  2:3403182/1‑93 (MQ=255)
                                                                        cccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACag                  <  2:1715141/93‑1 (MQ=255)
                                                                         ccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACaga                 <  2:5602375/93‑1 (MQ=255)
                                                                          cATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACagag                <  2:6434318/93‑1 (MQ=255)
                                                                          cATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACagag                <  1:4950085/93‑1 (MQ=255)
                                                                            tCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGaa              >  2:4189397/1‑93 (MQ=255)
                                                                            tCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGaa              <  2:473515/93‑1 (MQ=255)
                                                                              aCGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATa            >  1:6514178/1‑93 (MQ=255)
                                                                                gATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGc          >  2:5934771/1‑93 (MQ=255)
                                                                                  tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc        >  2:12380/1‑93 (MQ=255)
                                                                                  tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc        <  1:2193047/93‑1 (MQ=255)
                                                                                  tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc        <  2:1185894/93‑1 (MQ=255)
                                                                                  tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc        >  1:897997/1‑93 (MQ=255)
                                                                                   ggCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACg       <  2:5721352/93‑1 (MQ=255)
                                                                                     cGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGa     <  2:3582562/93‑1 (MQ=255)
                                                                                      gACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGaa    >  1:3156228/1‑93 (MQ=255)
                                                                                       acacGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGaaa   <  2:2085467/93‑1 (MQ=255)
                                                                                       acacGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGaaa   <  2:1648664/93‑1 (MQ=255)
                                                                                       acacGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGaaa   >  2:3540666/1‑93 (MQ=255)
                                                                                        cacGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAag  >  1:1629254/1‑93 (MQ=255)
                                                                                            |                                                                                        
ACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAG  >  NC_000913/552208‑552388

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