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% 882.5 / NA 242A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (133/109);  total (133/109)
*NC_0009133,958,8780CG100.0% 886.6 / NA 243A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (133/110);  total (133/110)

Only 100 of 246 total aligned reads displayed.

TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATG  >  NC_000913/3958790‑3958964
                                                                                       ||                                                                                      
   tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTACTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                 <  1:1141823/93‑1 (MQ=255)
     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               <  2:3922185/93‑1 (MQ=255)
      gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                              <  1:4199455/93‑1 (MQ=255)
         cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                           <  2:540391/93‑1 (MQ=255)
            tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                        >  1:4118442/1‑93 (MQ=255)
             atgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                                       >  1:3193086/1‑93 (MQ=255)
              tgatgCACAGCTGAAAGCCTGTAACGTGCTCGACTTCGATGAGCTGATGTTATGGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                      <  2:4039139/93‑1 (MQ=255)
              tgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                      >  1:5277253/1‑93 (MQ=255)
               tatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                     <  2:3334848/92‑1 (MQ=255)
               gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCTGCAAg                                                                     >  1:6142755/1‑93 (MQ=255)
                atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGGTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                    >  2:5932178/1‑93 (MQ=255)
                 tgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                                   >  1:1262602/1‑93 (MQ=255)
                 tgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                                   <  2:5037723/93‑1 (MQ=255)
                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                  >  1:4090841/1‑93 (MQ=255)
                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                  >  2:6083033/1‑93 (MQ=255)
                    acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                >  2:3232656/1‑93 (MQ=255)
                    acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                >  2:2937492/1‑93 (MQ=255)
                    acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                >  2:2027038/1‑93 (MQ=255)
                    acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                <  1:5607163/93‑1 (MQ=255)
                     cacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTgg                                                               <  2:928910/93‑1 (MQ=255)
                       ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                                             <  1:4381376/93‑1 (MQ=255)
                       ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                                             >  2:3424629/1‑93 (MQ=255)
                       ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                                             <  1:3973164/93‑1 (MQ=255)
                       ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                                             <  2:821340/93‑1 (MQ=255)
                       ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                                             >  2:5555570/1‑93 (MQ=255)
                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                           <  1:497654/93‑1 (MQ=255)
                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                           >  1:3094112/1‑93 (MQ=255)
                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                           >  2:2347697/1‑93 (MQ=255)
                          gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGaa                                                          >  1:1358530/1‑93 (MQ=255)
                          gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGaa                                                          <  2:5592176/93‑1 (MQ=255)
                            aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                        >  2:2052971/1‑93 (MQ=255)
                             aGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaa                                                       >  2:575668/1‑93 (MQ=255)
                              gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaa                                                      <  1:3786916/93‑1 (MQ=255)
                               ccTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaaa                                                     <  2:4012296/93‑1 (MQ=255)
                                cTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAt                                                    >  2:1335242/1‑93 (MQ=255)
                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  <  1:3711616/93‑1 (MQ=255)
                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  >  1:5864304/1‑93 (MQ=255)
                                    tgcGTTCTCGTCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                <  2:5485725/91‑1 (MQ=255)
                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                <  1:2648000/93‑1 (MQ=255)
                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                >  1:1923024/1‑93 (MQ=255)
                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                <  1:6110265/93‑1 (MQ=255)
                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCAGCAAGCGCTGGCAGAACAAAATTCGc                                                >  2:1141461/1‑93 (MQ=255)
                                     aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                                               <  2:4696122/93‑1 (MQ=255)
                                     aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                                               >  1:868917/1‑93 (MQ=255)
                                           tcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                         <  1:4268352/93‑1 (MQ=255)
                                           tcGACTTCGATGATCTGATTCTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                         >  2:1753819/1‑93 (MQ=255)
                                            cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                        <  2:940273/93‑1 (MQ=255)
                                            cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                        <  1:1722849/93‑1 (MQ=255)
                                             gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                                       <  1:5377814/93‑1 (MQ=255)
                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    >  1:1040369/1‑93 (MQ=255)
                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    >  1:5217194/1‑93 (MQ=255)
                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                <  1:602050/93‑1 (MQ=255)
                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                >  2:6140449/1‑93 (MQ=255)
                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                >  2:5172391/1‑93 (MQ=255)
                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                <  1:4338102/93‑1 (MQ=255)
                                                      gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                              >  1:2981461/1‑93 (MQ=255)
                                                      gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                              >  2:3511890/1‑93 (MQ=255)
                                                       atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                             >  1:3613419/1‑93 (MQ=255)
                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                           <  1:5982553/93‑1 (MQ=255)
                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                           <  1:1634385/93‑1 (MQ=255)
                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                           <  2:2305832/93‑1 (MQ=255)
                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                          >  1:1764367/1‑93 (MQ=255)
                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                          <  2:4606036/93‑1 (MQ=255)
                                                           taTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                         <  2:1707499/92‑1 (MQ=255)
                                                           gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATTAGTATCa                         >  2:1610262/1‑93 (MQ=255)
                                                            aTTTTATTGGCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        <  1:4569962/93‑1 (MQ=255)
                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        >  1:1546500/1‑93 (MQ=255)
                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        >  1:4829050/1‑93 (MQ=255)
                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        >  1:2799857/1‑93 (MQ=255)
                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        >  1:4226461/1‑93 (MQ=255)
                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       >  2:5864483/1‑93 (MQ=255)
                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       >  2:1504546/1‑93 (MQ=255)
                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       >  1:3304787/1‑93 (MQ=255)
                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      >  1:761973/1‑93 (MQ=255)
                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                     <  2:797576/93‑1 (MQ=255)
                                                                 aTTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacac                   <  1:5806252/93‑1 (MQ=255)
                                                                 aTTGCCGACGTTGCTGCCGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacac                   >  2:4329616/1‑93 (MQ=255)
                                                                  ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacacc                  >  1:4702941/1‑93 (MQ=255)
                                                                   tGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacacca                 <  2:2181815/93‑1 (MQ=255)
                                                                    gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaa                <  2:4399750/93‑1 (MQ=255)
                                                                    gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaa                >  1:2205516/1‑93 (MQ=255)
                                                                      cGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacc              >  1:4630512/1‑92 (MQ=255)
                                                                      cGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaaca              <  2:1206904/93‑1 (MQ=255)
                                                                       gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac             >  1:5623098/1‑93 (MQ=255)
                                                                       gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac             >  2:2782292/1‑93 (MQ=255)
                                                                       gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac             >  2:5709680/1‑93 (MQ=255)
                                                                        aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacc            >  1:2417310/1‑93 (MQ=255)
                                                                         cGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacca           <  2:3366333/93‑1 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  2:5311941/1‑93 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  2:1053323/1‑93 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  2:4851372/1‑93 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  1:4885856/1‑93 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  1:5260873/1‑93 (MQ=255)
                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  1:3435967/1‑93 (MQ=255)
                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc        <  2:420532/93‑1 (MQ=255)
                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc        >  1:93875/1‑93 (MQ=255)
                                                                             actgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcca       <  1:1160461/92‑1 (MQ=255)
                                                                              ctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagccag      <  1:2516138/93‑1 (MQ=255)
                                                                                gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTa    <  1:5197222/93‑1 (MQ=255)
                                                                                 ctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTAt   <  2:3685539/93‑1 (MQ=255)
                                                                                       ||                                                                                      
TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATG  >  NC_000913/3958790‑3958964

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