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,8770GC99.8% 1442.0 / NA 439A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (196/242);  minor base A (1/0);  total (197/242)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,958,8780CG99.8% 1436.9 / ‑6.6 440A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (197/242);  total (197/243)
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.99e-01

Only 100 of 447 total aligned reads displayed.

TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAGCTGG  >  NC_000913/3958790‑3958970
                                                                                       ||                                                                                            
  ttATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATga                                                                                        <  1:6380361/93‑1 (MQ=255)
   tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                       >  2:9229977/1‑93 (MQ=255)
    aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                                      <  1:1201403/93‑1 (MQ=255)
    aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                                      <  1:3470108/93‑1 (MQ=255)
     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                                     >  1:5558918/1‑93 (MQ=255)
     tGGGCTGTATGATGCACACCGGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                                     <  2:7165466/93‑1 (MQ=255)
      gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                                    <  1:7981934/93‑1 (MQ=255)
      gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                                    <  2:3766219/93‑1 (MQ=255)
       ggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                                                   <  2:6033028/93‑1 (MQ=255)
        gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                                                  <  2:8878573/93‑1 (MQ=255)
        gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                                                  >  2:6667462/1‑93 (MQ=255)
         cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                                 <  2:5678639/93‑1 (MQ=255)
          tGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                                                <  1:6262193/93‑1 (MQ=255)
          tGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                                                >  1:2828287/1‑93 (MQ=255)
          tGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                                                <  1:455249/93‑1 (MQ=255)
             atgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                                             >  1:1959369/1‑93 (MQ=255)
                 tgCACCCCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                                         <  2:2718445/93‑1 (MQ=255)
                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                        <  2:1793418/93‑1 (MQ=255)
                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                        <  2:865668/93‑1 (MQ=255)
                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                        <  2:2758184/93‑1 (MQ=255)
                   cacaCCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCt                                                                       >  1:7136285/1‑93 (MQ=255)
                   cacaCCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCt                                                                       <  1:9546299/93‑1 (MQ=255)
                     cacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTgg                                                                     >  1:3743439/1‑93 (MQ=255)
                      acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGc                                                                    <  2:3322764/93‑1 (MQ=255)
                       ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                                                   <  1:2280637/93‑1 (MQ=255)
                        cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                                  <  2:6276907/93‑1 (MQ=255)
                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                                 >  2:5789022/1‑93 (MQ=255)
                          gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGGCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGaa                                                                >  2:1060207/1‑93 (MQ=255)
                           aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAAc                                                               >  1:8550631/1‑93 (MQ=255)
                            aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                              <  1:8433120/93‑1 (MQ=255)
                             aGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaa                                                             >  2:2139026/1‑93 (MQ=255)
                               ccTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaaa                                                           <  2:3113224/93‑1 (MQ=255)
                                cTGTAACGTGCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAt                                                          <  2:5522505/93‑1 (MQ=255)
                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                         <  1:6148456/93‑1 (MQ=255)
                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                         <  1:6523672/93‑1 (MQ=255)
                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                         >  1:9938179/1‑93 (MQ=255)
                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                        >  2:1355987/1‑93 (MQ=255)
                                   tAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCg                                                       <  1:5281875/93‑1 (MQ=255)
                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                      >  1:87598/1‑93 (MQ=255)
                                     aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                                                     >  1:7090676/1‑93 (MQ=255)
                                      cGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTa                                                    <  1:3985816/93‑1 (MQ=255)
                                       gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                                   <  2:9710799/93‑1 (MQ=255)
                                        ttCTCGACTTCGATGCTCTGATTTTATTGCCGACGTTGCTGCTGCTACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                                  <  1:8937844/93‑1 (MQ=255)
                                        ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                                  <  1:495623/93‑1 (MQ=255)
                                        ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                                  >  1:1413436/1‑93 (MQ=255)
                                         tctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATct                                                 <  2:5073983/93‑1 (MQ=255)
                                           tcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                               >  1:7724333/1‑93 (MQ=255)
                                           tcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                               <  1:3284598/93‑1 (MQ=255)
                                            cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                              <  1:9062100/93‑1 (MQ=255)
                                             gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                                             <  2:9506503/93‑1 (MQ=255)
                                               cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                                           <  1:596714/93‑1 (MQ=255)
                                               cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                                           >  2:1242226/1‑93 (MQ=255)
                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACACAATTCGCTATCTGCtggag                                          >  1:6624832/1‑91 (MQ=255)
                                                 tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                         >  2:4387686/1‑93 (MQ=255)
                                                 tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                         >  1:5137974/1‑93 (MQ=255)
                                                 tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                         >  1:2075455/1‑93 (MQ=255)
                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                        <  2:7841050/93‑1 (MQ=255)
                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                        >  2:4717499/1‑93 (MQ=255)
                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                        <  2:2363246/93‑1 (MQ=255)
                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                        <  1:8289486/93‑1 (MQ=255)
                                                   gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt                                       >  2:5412491/1‑93 (MQ=255)
                                                   gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt                                       <  1:1392369/93‑1 (MQ=255)
                                                   gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGCGGAt                                       >  2:4498053/1‑93 (MQ=255)
                                                      gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                                    >  1:1718085/1‑93 (MQ=255)
                                                       atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                                   >  2:4160924/1‑93 (MQ=255)
                                                       atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                                   >  2:1482326/1‑93 (MQ=255)
                                                       atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                                   >  2:2967690/1‑93 (MQ=255)
                                                        tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa                                  <  1:4435340/93‑1 (MQ=255)
                                                        tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa                                  >  2:4124781/1‑93 (MQ=255)
                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                                 <  1:4328733/93‑1 (MQ=255)
                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                                 <  1:2535437/93‑1 (MQ=255)
                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                                 >  1:8064771/1‑93 (MQ=255)
                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                                >  1:4141778/1‑93 (MQ=255)
                                                           gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                               <  2:6488376/93‑1 (MQ=255)
                                                           gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGAATCa                               >  2:80402/1‑93 (MQ=255)
                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                              <  1:1084134/93‑1 (MQ=255)
                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                              >  1:7346292/1‑93 (MQ=255)
                                                            aTTTTATTGCCGACGTTGATGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                              >  1:270262/1‑93 (MQ=255)
                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                            <  2:6957287/93‑1 (MQ=255)
                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                           >  2:9149587/1‑93 (MQ=255)
                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                           <  2:1768132/93‑1 (MQ=255)
                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                           >  1:5815362/1‑93 (MQ=255)
                                                                  ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacacc                        <  2:7437760/93‑1 (MQ=255)
                                                                  ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaca                        >  2:8994940/1‑92 (MQ=255)
                                                                  ttGCCGACGTTGCTGCTGCAACCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacacc                        <  1:3010965/93‑1 (MQ=255)
                                                                   tGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacacca                       >  2:5917211/1‑93 (MQ=255)
                                                                    gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaa                      <  2:6227464/93‑1 (MQ=255)
                                                                       gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac                   <  1:7606423/93‑1 (MQ=255)
                                                                        aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTTGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacc                  <  1:9444143/93‑1 (MQ=255)
                                                                         cGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGgcaccaacacca                 >  2:4044905/1‑93 (MQ=255)
                                                                         cGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGAacaccaacacca                 >  1:1001502/1‑93 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc               >  1:9689856/1‑93 (MQ=255)
                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc              <  2:8536616/93‑1 (MQ=255)
                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc              <  1:7772128/93‑1 (MQ=255)
                                                                             gctgctgcAACGCAATGAAGAAGTCCGCAAGCCCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACAAACAccagcca             <  1:5289529/93‑1 (MQ=255)
                                                                              ctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagccag            >  1:6877172/1‑93 (MQ=255)
                                                                              ctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagccag            <  1:2146992/93‑1 (MQ=255)
                                                                               tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt           >  2:3462581/1‑93 (MQ=255)
                                                                               tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt           >  2:5244188/1‑93 (MQ=255)
                                                                                gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTa          <  2:1195158/93‑1 (MQ=255)
                                                                                       ||                                                                                            
TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAGCTGG  >  NC_000913/3958790‑3958970

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