Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,565,832 (C)5→6 coding (1053/1404 nt) eutA ← reactivating factor for ethanolamine ammonia‑lyase EutBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,565,8271.C99.6% 906.4 / ‑3.1 222coding (1058/1404 nt)eutAreactivating factor for ethanolamine ammonia‑lyase EutBC
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (112/109);  total (112/110)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 253 total aligned reads displayed.

CAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGC  >  NC_000913/2565727‑2565922
                                                                                                     |                                                                                               
cAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAgg                                                                                                  >  2:3322135/1‑101 (MQ=255)
   tACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccc                                                                                              >  2:3322138/1‑101 (MQ=255)
    aCCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGcccc                                                                                             <  2:3321965/101‑1 (MQ=255)
     ccGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccc                                                                                            >  2:3322142/1‑101 (MQ=255)
     ccGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccc                                                                                            <  1:3321680/101‑1 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGcccccca                                                                                          >  2:3322145/1‑99 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGcccccca                                                                                          <  2:3321753/101‑3 (MQ=255)
       gcggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGcccccca                                                                                          <  1:3321920/101‑3 (MQ=255)
        cggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccccac                                                                                         >  1:3322146/1‑98 (MQ=255)
         ggcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccccacc                                                                                        >  1:3322148/1‑97 (MQ=255)
          gcgTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGccccccacca                                                                                       <  2:3321633/101‑6 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGAt                                                                                     <  1:3321880/101‑1 (MQ=255)
             tAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGAt                                                                                     >  2:3322152/1‑101 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATcc                                                                                   <  2:3321855/101‑1 (MQ=255)
               aCGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATcc                                                                                   <  2:3321743/101‑1 (MQ=255)
                cGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCg                                                                                  >  1:3322156/1‑101 (MQ=255)
                 gCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGt                                                                                 <  1:3321744/101‑1 (MQ=255)
                 gCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGt                                                                                 >  2:3322157/1‑101 (MQ=255)
                    caGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTtc                                                                              <  1:3321899/101‑1 (MQ=255)
                     aGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGttt                                                                               <  1:3321816/99‑1 (MQ=255)
                     aGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTtca                                                                             <  2:3321845/101‑1 (MQ=255)
                      ggCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTtcat                                                                            >  2:3322161/1‑101 (MQ=255)
                        cAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCGGGCCCCCCACCAGATCCGTTtcatca                                                                          <  2:3321863/101‑1 (MQ=255)
                        cAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTtcatca                                                                          <  1:3321861/101‑1 (MQ=255)
                        cAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTtcatca                                                                          >  2:3322163/1‑101 (MQ=255)
                         aGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCaa                                                                         >  2:3322166/1‑101 (MQ=255)
                           cGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTTAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATc                                                                       <  2:3321798/101‑1 (MQ=255)
                           cGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATc                                                                       <  2:3321771/101‑1 (MQ=255)
                            gAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCg                                                                      >  2:3322168/1‑101 (MQ=255)
                            gAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCg                                                                      >  1:3322172/1‑101 (MQ=255)
                             aGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCgg                                                                     >  2:3322173/1‑101 (MQ=255)
                              ggCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCggg                                                                    <  2:3321890/101‑1 (MQ=255)
                                 gggAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATc                                                                 >  2:3322177/1‑101 (MQ=255)
                                 gggAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATc                                                                 <  1:3321701/101‑1 (MQ=255)
                                    aaGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGcc                                                              <  1:3321903/101‑1 (MQ=255)
                                    aaGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGcc                                                              <  2:3321833/101‑1 (MQ=255)
                                       cgcCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCAcc                                                           >  2:3322180/1‑101 (MQ=255)
                                        gcCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCg                                                          <  1:3321866/101‑1 (MQ=255)
                                        gcCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCg                                                          >  2:3322182/1‑101 (MQ=255)
                                        gcCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCg                                                          <  1:3321595/101‑1 (MQ=255)
                                        gcCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCg                                                          <  2:3321724/101‑1 (MQ=255)
                                         ccAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCgg                                                         >  1:3322184/1‑101 (MQ=255)
                                          cAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAgg                                                                                                  >  2:3322187/1‑59 (MQ=255)
                                           aGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa                                                       <  2:3322036/101‑1 (MQ=255)
                                           aGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCa                                                       <  1:3321784/101‑1 (MQ=255)
                                            gCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCaa                                                      <  2:3321685/101‑1 (MQ=255)
                                             cACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCaaa                                                     <  1:3322018/101‑1 (MQ=255)
                                             cACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCaaa                                                     <  1:3321925/101‑1 (MQ=255)
                                              aCGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAAt                                                    >  2:3322189/1‑101 (MQ=255)
                                                 tACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCggg                                                                    >  2:3322192/1‑82 (MQ=255)
                                                 tACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTgc                                                 <  1:3322060/101‑1 (MQ=255)
                                                 tACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTgc                                                 <  1:3321962/101‑1 (MQ=255)
                                                  aCGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTgcg                                                <  1:3321839/101‑1 (MQ=255)
                                                   cgcgTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTgcgc                                               >  1:3322194/1‑101 (MQ=255)
                                                   cgcgTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTgcgc                                               <  1:3321788/101‑1 (MQ=255)
                                                    gcgTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCa                                              >  2:3322196/1‑101 (MQ=255)
                                                     cgTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAg                                             >  2:3322198/1‑101 (MQ=255)
                                                      gTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGc                                            >  2:3322201/1‑101 (MQ=255)
                                                      gTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGc                                            <  1:3321774/101‑1 (MQ=255)
                                                       tCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCg                                           >  1:3322204/1‑101 (MQ=255)
                                                        cAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCgg                                          <  2:3321914/101‑1 (MQ=255)
                                                         aGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGc                                         >  2:3322206/1‑101 (MQ=255)
                                                          gTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCa                                        <  2:3321881/101‑1 (MQ=255)
                                                          gTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCa                                        >  2:3322210/1‑101 (MQ=255)
                                                            tttGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGt                                      >  2:3322214/1‑101 (MQ=255)
                                                            tttGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGt                                      >  1:3322212/1‑101 (MQ=255)
                                                               gggATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGt                                   <  2:3321773/101‑1 (MQ=255)
                                                               gggATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGt                                   >  2:3322218/1‑101 (MQ=255)
                                                               gggATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGt                                   >  1:3322219/1‑101 (MQ=255)
                                                                  aTCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACg                                <  2:3321820/101‑1 (MQ=255)
                                                                  aTCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACg                                >  2:3322222/1‑101 (MQ=255)
                                                                      aGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCt                            >  1:3322225/1‑101 (MQ=255)
                                                                       gATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTc                           >  2:3322228/1‑101 (MQ=255)
                                                                        aTCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCa                          >  1:3322229/1‑101 (MQ=255)
                                                                         tCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCaa                         <  2:3321806/101‑1 (MQ=255)
                                                                           cAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGc                       >  1:3322233/1‑101 (MQ=255)
                                                                           cAGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGc                       <  1:3321836/101‑1 (MQ=255)
                                                                            aGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGcc                      <  1:3321821/101‑1 (MQ=255)
                                                                            aGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGcc                      <  1:3321811/101‑1 (MQ=255)
                                                                            aGCTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGcc                      >  2:3322234/1‑101 (MQ=255)
                                                                              cTGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAg                    >  2:3322238/1‑101 (MQ=255)
                                                                               tGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGa                   >  2:3322242/1‑101 (MQ=255)
                                                                               tGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGa                   >  1:3322239/1‑101 (MQ=255)
                                                                               tGAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGa                   <  1:3321960/101‑1 (MQ=255)
                                                                                gAATCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGAt                  <  2:3321871/101‑1 (MQ=255)
                                                                                 aaTCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGAtt                 <  1:3321945/101‑1 (MQ=255)
                                                                                  aTCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtg                >  1:3322246/1‑101 (MQ=255)
                                                                                   tCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtgt               >  1:3322248/1‑101 (MQ=255)
                                                                                   tCAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtgt               <  1:3321941/101‑1 (MQ=255)
                                                                                    cAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtgtg              <  2:3322012/101‑1 (MQ=255)
                                                                                    cAGCGCCTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtgtg              >  2:3322249/1‑101 (MQ=255)
                                                                                         ccTGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGcc         >  2:3322252/1‑101 (MQ=255)
                                                                                           tGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCag       <  2:3321987/101‑1 (MQ=255)
                                                                                           tGTTGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCag       >  1:3322256/1‑101 (MQ=255)
                                                                                              tGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCagaga    >  2:3322258/1‑101 (MQ=255)
                                                                                              tGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCagaga    <  2:3321760/101‑1 (MQ=255)
                                                                                              tGCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCagaga    <  1:3321635/101‑1 (MQ=255)
                                                                                               gCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCagagag   >  1:3322259/1‑101 (MQ=255)
                                                                                               gCCAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCagagag   <  1:3321751/101‑1 (MQ=255)
                                                                                                ccAGGCCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGc  >  1:3322264/1‑101 (MQ=255)
                                                                                                     |                                                                                               
CAGTACCGCGGCGTAACGCACAGGCAGCGAGGCGGGAAGCGCCAGCACGTACGCGTCAGTTTTGGGATCAAGATCCAGCTGAATCAGCGCCTGTTGCCAGGCCCCCACCAGATCCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGC  >  NC_000913/2565727‑2565922

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