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% 2135.7 / NA 538A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (254/284);  total (254/284)
*NC_0009133,958,8780CG99.8% 2122.3 / NA 539A196R (GCC→CGC) reprep helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (253/285);  minor base T (0/1);  total (253/286)
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 547 total aligned reads displayed.

CATTGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAGC  >  NC_000913/3958787‑3958967
                                                                                          ||                                                                                         
cATTGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAAt                                                                                       <  1:6994879/96‑1 (MQ=255)
cATTGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAAt                                                                                       <  2:13770261/96‑1 (MQ=255)
cATTGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGATCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAAt                                                                                       >  1:1921275/1‑96 (MQ=255)
 aTTGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATg                                                                                      <  1:11977675/96‑1 (MQ=255)
  ttGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATga                                                                                     >  1:7389987/1‑96 (MQ=255)
  ttGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATga                                                                                     >  1:3978922/1‑96 (MQ=255)
   tGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                    >  2:1666905/1‑96 (MQ=255)
     ttATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                                  <  1:9494543/96‑1 (MQ=255)
     ttATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                                  >  1:12972924/1‑96 (MQ=255)
      tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                                 >  2:9423994/1‑96 (MQ=255)
       aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                                                <  1:3763520/96‑1 (MQ=255)
        tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                                               >  2:3131273/1‑96 (MQ=255)
         gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                              >  2:10012261/1‑96 (MQ=255)
         gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                              <  1:762183/96‑1 (MQ=255)
         gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                              <  2:7370643/96‑1 (MQ=255)
           gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                            <  2:7888019/96‑1 (MQ=255)
           gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                            >  1:13637439/1‑96 (MQ=255)
           gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                            >  2:10658864/1‑96 (MQ=255)
            cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                           >  2:13246723/1‑96 (MQ=255)
             tGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                                          >  2:6609853/1‑96 (MQ=255)
              gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                         >  2:6810100/1‑96 (MQ=255)
              gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                         >  1:14026098/1‑96 (MQ=255)
              gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                         >  1:4174063/1‑96 (MQ=255)
              gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                         >  2:6193593/1‑96 (MQ=255)
               tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                        >  1:2677683/1‑96 (MQ=255)
                atgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                       <  1:9286779/96‑1 (MQ=255)
                  gatgCACACCTGAAAGCCTGGAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                     >  2:93312/1‑96 (MQ=255)
                   atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCt                                                                    <  1:1329183/96‑1 (MQ=255)
                    tgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGCCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                   <  2:7719048/96‑1 (MQ=255)
                      cacaCCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGc                                                                 >  2:3361070/1‑96 (MQ=255)
                        cacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                               <  2:4897408/96‑1 (MQ=255)
                         acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAGGAAGTCCGCAAGCGCTGGCAGa                                                              >  1:5549276/1‑96 (MQ=255)
                         acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                              >  2:7587425/1‑96 (MQ=255)
                         acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                              <  2:13533324/96‑1 (MQ=255)
                         acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                              >  2:5895217/1‑96 (MQ=255)
                         acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                              >  1:12310602/1‑96 (MQ=255)
                           cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAAc                                                            <  1:11018849/96‑1 (MQ=255)
                           cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAAc                                                            >  2:12437175/1‑96 (MQ=255)
                           cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAAc                                                            <  1:12097907/96‑1 (MQ=255)
                            tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                           >  2:2662698/1‑96 (MQ=255)
                             gAAATCCTGTAACGTTCTCGACTTCGATGATCTGATTTCTTTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaa                                                          >  1:12599136/1‑96 (MQ=255)
                             gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaa                                                          <  2:7204038/96‑1 (MQ=255)
                               aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaaa                                                        <  1:6414552/96‑1 (MQ=255)
                               aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaaa                                                        >  1:13919929/1‑96 (MQ=255)
                                 gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                      >  2:2928686/1‑96 (MQ=255)
                                 gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                      >  1:1436084/1‑96 (MQ=255)
                                  ccTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                     <  1:6705473/96‑1 (MQ=255)
                                   cTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCg                                                    <  2:1300266/96‑1 (MQ=255)
                                    tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                   <  2:10545205/96‑1 (MQ=255)
                                     gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                                                  <  2:11315070/96‑1 (MQ=255)
                                      tAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTa                                                 >  2:10887990/1‑96 (MQ=255)
                                      tAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCAGCAAGCGCTGGCAGAACAAAATTCGCTa                                                 >  2:14203878/1‑96 (MQ=255)
                                       aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                                <  1:5317953/96‑1 (MQ=255)
                                        aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                               <  2:3829092/96‑1 (MQ=255)
                                           ttCTAGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGGCCGCAAGCGCTGGCAGAAGAAAATTTGCTATctgc                                            >  1:4500899/1‑96 (MQ=255)
                                               cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                                        <  2:5459658/96‑1 (MQ=255)
                                                gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                       >  2:6410099/1‑96 (MQ=255)
                                                  cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                     <  1:11596891/96‑1 (MQ=255)
                                                     cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGa                                  <  2:9003946/96‑1 (MQ=255)
                                                     cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGa                                  >  2:8543927/1‑96 (MQ=255)
                                                     cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGa                                  >  2:8335480/1‑96 (MQ=255)
                                                      gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                                 <  1:8539259/96‑1 (MQ=255)
                                                       atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                                <  2:4699624/96‑1 (MQ=255)
                                                       atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt                                <  2:225180/96‑1 (MQ=255)
                                                        tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa                               <  2:10274633/96‑1 (MQ=255)
                                                        tgatCTGATTTTATTGCCGACGTAGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGATGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa                               <  2:12637620/96‑1 (MQ=255)
                                                        tgatCTGAGTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCGGAACAAAATGCGCTATCTGCTGGTGGATGAGTa                               >  1:9930689/1‑96 (MQ=255)
                                                         gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCTATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGtct                              <  1:6227444/96‑3 (MQ=255)
                                                         gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                              <  1:2928686/96‑1 (MQ=255)
                                                          atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTaac                             <  1:10606373/96‑3 (MQ=255)
                                                          atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                             <  2:12805552/96‑1 (MQ=255)
                                                           tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                            <  1:5909874/96‑1 (MQ=255)
                                                           tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                            <  1:7776720/96‑1 (MQ=255)
                                                            cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCac                           <  1:8821692/96‑2 (MQ=255)
                                                            cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                           <  1:8216736/96‑1 (MQ=255)
                                                             tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                          <  1:8714834/96‑1 (MQ=255)
                                                               aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                        <  2:4384939/96‑1 (MQ=255)
                                                                ttttNTTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATATGCTGGTGGATGAGAATCAGGaca                       >  2:4195837/1‑96 (MQ=255)
                                                                       gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGTACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaactc                <  1:14221258/96‑3 (MQ=255)
                                                                       gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac                <  1:3233336/96‑1 (MQ=255)
                                                                         cGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacca              <  1:6064050/96‑1 (MQ=255)
                                                                         cGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAAAAAAATTCGCTATCTGCTGGTGGATGAGTATAAGGacaccaacacca              >  1:499289/1‑96 (MQ=255)
                                                                          gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccag             <  1:10338039/96‑1 (MQ=255)
                                                                          gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccag             >  2:355873/1‑96 (MQ=255)
                                                                          gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccag             <  2:11920117/96‑1 (MQ=255)
                                                                          gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccag             <  2:1165472/96‑1 (MQ=255)
                                                                             gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGTACACCAACAccagcca          >  1:230279/1‑96 (MQ=255)
                                                                             gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcca          <  2:2533894/96‑1 (MQ=255)
                                                                             gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcca          >  2:8056554/1‑96 (MQ=255)
                                                                               tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt        >  2:5155423/1‑96 (MQ=255)
                                                                                gctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTa       >  1:1832737/1‑96 (MQ=255)
                                                                                  tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTacg     <  1:11285401/96‑3 (MQ=255)
                                                                                  tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATg     >  2:5316200/1‑96 (MQ=255)
                                                                                  tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATg     <  2:12141925/96‑1 (MQ=255)
                                                                                  tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATg     >  1:13481863/1‑96 (MQ=255)
                                                                                    ctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACCCCAACCCCAGCCAGTATGAg   <  1:1006890/96‑1 (MQ=255)
                                                                                    ctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAg   <  2:981525/96‑1 (MQ=255)
                                                                                    ctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAg   <  2:1408156/96‑1 (MQ=255)
                                                                                     tgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAGc  >  1:9670482/1‑96 (MQ=255)
                                                                                     tgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAGc  >  2:9583078/1‑96 (MQ=255)
                                                                                          ||                                                                                         
CATTGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATGAGC  >  NC_000913/3958787‑3958967

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