Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,958,877 2 bp→CG coding (586‑587/2022 nt) rep → rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,958,8770GC100.0% 561.8 / NA 171A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (77/94);  total (77/94)
*NC_0009133,958,8780CG99.4% 560.0 / NA 170A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (76/93);  minor base T (0/1);  total (76/94)
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.00e+00

Only 100 of 174 total aligned reads displayed.

GTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTAT  >  NC_000913/3958791‑3958963
                                                                                      ||                                                                                     
 ttATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATga                                                                                 >  1:3555555/1‑93 (MQ=255)
  tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                >  2:3556544/1‑93 (MQ=255)
  tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                <  1:2942134/93‑1 (MQ=255)
   aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTGGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                               <  2:649736/93‑1 (MQ=255)
   aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                               <  2:4130175/93‑1 (MQ=255)
   aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                               <  2:1825527/93‑1 (MQ=255)
    tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                              <  2:548711/93‑1 (MQ=255)
     gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                             <  2:1624592/93‑1 (MQ=255)
     gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                             >  2:4527376/1‑93 (MQ=255)
      ggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                                            >  2:259836/1‑93 (MQ=255)
      ggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                                            <  2:1467713/93‑1 (MQ=255)
       gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                                           <  2:2910340/93‑1 (MQ=255)
       gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                                           <  2:3681229/93‑1 (MQ=255)
        cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                          <  1:3251208/93‑1 (MQ=255)
        cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                          <  1:2839064/93‑1 (MQ=255)
        cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                          <  2:3182800/93‑1 (MQ=255)
        cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                          <  1:96935/93‑1 (MQ=255)
        cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                          <  1:4653308/93‑1 (MQ=255)
         tGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                                         <  2:116474/93‑1 (MQ=255)
          gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                        >  1:590428/1‑93 (MQ=255)
          gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                        >  2:2795756/1‑93 (MQ=255)
          gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                        >  2:3603713/1‑93 (MQ=255)
           tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                       <  2:1978623/93‑1 (MQ=255)
            atgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                                      >  1:4518565/1‑93 (MQ=255)
             tgatgCACTCGTGAAAGCCTGTATCGTTCTCGACTTCGATAATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                     <  2:4375599/93‑1 (MQ=255)
              gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAGCTCAATGAAGAAGTCCGCAAg                                                                    <  2:809166/93‑1 (MQ=255)
              gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                    <  2:2683874/93‑1 (MQ=255)
               atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                   <  2:1724127/93‑1 (MQ=255)
               atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                   >  2:2946443/1‑93 (MQ=255)
               atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                   <  2:2474121/93‑1 (MQ=255)
               atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                   >  1:356225/1‑93 (MQ=255)
                tgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                                  <  2:4440996/93‑1 (MQ=255)
                  cacaCCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCt                                                                <  1:4153172/93‑1 (MQ=255)
                  cacaCCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCt                                                                <  2:472528/93‑1 (MQ=255)
                   acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                               >  1:4247074/1‑93 (MQ=255)
                    cacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTgg                                                              <  1:1056167/93‑1 (MQ=255)
                       cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                           >  1:107967/1‑93 (MQ=255)
                       cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                           >  1:74916/1‑93 (MQ=255)
                         gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGaa                                                         <  1:807288/93‑1 (MQ=255)
                          aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGCAc                                                        >  1:1895076/1‑93 (MQ=255)
                             gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaa                                                     >  1:2996905/1‑93 (MQ=255)
                               cTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAt                                                   >  1:627009/1‑93 (MQ=255)
                               cTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGACCGCAAGCGCTGGCAGAACAAAAt                                                   >  1:627008/1‑93 (MQ=255)
                                tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                  >  1:3303511/1‑93 (MQ=255)
                                tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                  <  1:309629/93‑1 (MQ=255)
                                tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                  <  2:1993344/93‑1 (MQ=255)
                                 gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                 >  2:799425/1‑93 (MQ=255)
                                 gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                 <  1:4409377/93‑1 (MQ=255)
                                 gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                 <  1:2552279/93‑1 (MQ=255)
                                   aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                               <  1:4735170/93‑1 (MQ=255)
                                   aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                               <  1:4270377/93‑1 (MQ=255)
                                   aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCGGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                               <  2:570761/93‑1 (MQ=255)
                                    aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                                              >  1:3798844/1‑93 (MQ=255)
                                     cGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTa                                             <  1:517635/93‑1 (MQ=255)
                                     cGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTa                                             >  1:332880/1‑93 (MQ=255)
                                      gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                            >  1:227247/1‑93 (MQ=255)
                                      gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                            >  2:757447/1‑93 (MQ=255)
                                      gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                            >  1:3494457/1‑93 (MQ=255)
                                       ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                           >  1:1527969/1‑93 (MQ=255)
                                        tctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATct                                          <  2:43321/93‑1 (MQ=255)
                                         ctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctg                                         >  2:634875/1‑93 (MQ=255)
                                          tcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                        >  2:2451269/1‑93 (MQ=255)
                                           cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                       >  1:2550164/1‑93 (MQ=255)
                                           cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                       <  1:1742428/93‑1 (MQ=255)
                                           cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                       >  1:406392/1‑93 (MQ=255)
                                            gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                                      <  2:2817244/93‑1 (MQ=255)
                                            gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                                      >  2:1257868/1‑93 (MQ=255)
                                            gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                                      >  1:3966345/1‑93 (MQ=255)
                                             aCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                                     <  2:90398/93‑1 (MQ=255)
                                              cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                                    <  2:4733897/93‑1 (MQ=255)
                                               ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                   <  2:672749/93‑1 (MQ=255)
                                                tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                  <  1:2324539/93‑1 (MQ=255)
                                                tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                  <  2:3207593/93‑1 (MQ=255)
                                                 cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                 <  2:879771/93‑1 (MQ=255)
                                                  gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt                                >  2:763206/1‑93 (MQ=255)
                                                    tgcgcTGAGTTTATTGCCGACGATGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGTAGAACAAAATTCGCTATCTGCTGGTGGATGa                              <  2:780102/89‑1 (MQ=255)
                                                    tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGa                              <  2:774947/93‑1 (MQ=255)
                                                     gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                             <  1:12206/93‑1 (MQ=255)
                                                      atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                            <  2:2078876/93‑1 (MQ=255)
                                                       tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATgggga                           >  2:2217501/1‑89 (MQ=255)
                                                        cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                          <  2:2893372/93‑1 (MQ=255)
                                                        cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                          >  1:1675765/1‑93 (MQ=255)
                                                         tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                         >  1:2678186/1‑93 (MQ=255)
                                                          gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                        <  2:1035163/93‑1 (MQ=255)
                                                           aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                       <  2:2367215/93‑1 (MQ=255)
                                                           aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                       >  1:1386265/1‑93 (MQ=255)
                                                           aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                       >  2:1384576/1‑93 (MQ=255)
                                                             tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                     <  2:1229553/93‑1 (MQ=255)
                                                              ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                    <  2:323488/93‑1 (MQ=255)
                                                                aTTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacac                  >  1:140728/1‑93 (MQ=255)
                                                                   gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaa               <  1:1773675/93‑1 (MQ=255)
                                                                      gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac            >  2:161502/1‑93 (MQ=255)
                                                                       aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacc           <  1:1788256/93‑1 (MQ=255)
                                                                        cGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacca          >  2:267953/1‑93 (MQ=255)
                                                                           tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc       >  1:21656/1‑93 (MQ=255)
                                                                           tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc       >  1:2968543/1‑93 (MQ=255)
                                                                            gctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGCCACCAACAccagcca      <  2:4048264/93‑1 (MQ=255)
                                                                             ctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagccag     >  1:3490663/1‑93 (MQ=255)
                                                                              tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt    <  2:4709045/93‑1 (MQ=255)
                                                                                ctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTAt  <  1:1946073/93‑1 (MQ=255)
                                                                                      ||                                                                                     
GTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTAT  >  NC_000913/3958791‑3958963

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