Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 329,953 A→G T423A (ACC→GCC)  betT → high‑affinity choline transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913329,9530AG99.6% 733.5 / ‑5.7 230T423A (ACC→GCC) betThigh‑affinity choline transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (99/130);  total (99/131)
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 232 total aligned reads displayed.

GAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTC  >  NC_000913/329863‑330042
                                                                                          |                                                                                         
 aGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCt                                                                                        <  2:5821716/93‑1 (MQ=255)
  gcgcGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTc                                                                                       <  2:2423929/93‑1 (MQ=255)
  gcgcGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTc                                                                                       <  1:4423175/93‑1 (MQ=255)
   cgcgGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcg                                                                                      >  1:1934426/1‑93 (MQ=255)
    gcgGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcgg                                                                                     >  1:4955891/1‑93 (MQ=255)
       gCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcggcgg                                                                                  >  1:6306357/1‑93 (MQ=255)
        cTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGa                                                                                 <  2:3319931/93‑1 (MQ=255)
         ttCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGAc                                                                                >  1:3107423/1‑93 (MQ=255)
         ttCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGAc                                                                                >  1:168636/1‑93 (MQ=255)
          tCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACt                                                                               <  1:3797290/93‑1 (MQ=255)
          tCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACt                                                                               >  2:779289/1‑93 (MQ=255)
            tACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTcc                                                                             >  1:6049270/1‑93 (MQ=255)
             aCAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCg                                                                            <  1:6021937/93‑1 (MQ=255)
             aCAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCg                                                                            >  2:65631/1‑93 (MQ=255)
             aCAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCg                                                                            <  1:5112146/93‑1 (MQ=255)
             aCAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCAACATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCg                                                                            >  2:5846424/1‑93 (MQ=255)
              cAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCgg                                                                           <  2:2913044/93‑1 (MQ=255)
              cAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCgg                                                                           <  1:2061216/93‑1 (MQ=255)
               aGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCggg                                                                          <  2:5488078/93‑1 (MQ=255)
                gCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCgggg                                                                         >  2:2564176/1‑93 (MQ=255)
                 ccTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgc                                                                        <  2:1371824/93‑1 (MQ=255)
                 ccTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgc                                                                        >  1:965757/1‑93 (MQ=255)
                   tgctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgcgc                                                                      >  1:5219718/1‑93 (MQ=255)
                   tgctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgcgc                                                                      <  2:307495/93‑1 (MQ=255)
                   tgctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgcgc                                                                      <  1:3852167/93‑1 (MQ=255)
                    gctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCt                                                                     <  2:2171850/93‑1 (MQ=255)
                      tgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTgg                                                                   <  1:1065387/93‑1 (MQ=255)
                       ggCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGt                                                                  <  2:2049535/93‑1 (MQ=255)
                        gcgcAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTg                                                                 <  2:408633/93‑1 (MQ=255)
                            aGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGcgctggcgctgg                                                             >  1:2952847/1‑93 (MQ=255)
                            aGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTgg                                                             >  2:3603406/1‑93 (MQ=255)
                            aGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTgg                                                             <  1:1460736/93‑1 (MQ=255)
                            aGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTgg                                                             <  1:5168902/93‑1 (MQ=255)
                             gTATCCGGCTTTTCCCTTTAGCGCCTCCGTCGCCTCCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTgtg                                                            <  2:4220959/93‑3 (MQ=255)
                             gTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTggg                                                            >  1:2015302/1‑93 (MQ=255)
                                tCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGaa                                                         <  1:1620822/93‑1 (MQ=255)
                                tCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGaa                                                         <  1:2827041/93‑1 (MQ=255)
                                 ccGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAt                                                        <  2:1245776/93‑1 (MQ=255)
                                 ccGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAt                                                        <  2:3439192/93‑1 (MQ=255)
                                  cgtcGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAtt                                                       <  1:4878100/90‑1 (MQ=255)
                                  cGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAtt                                                       <  2:4367605/93‑1 (MQ=255)
                                  cGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAtt                                                       <  2:3386817/93‑1 (MQ=255)
                                   ggCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAttt                                                      >  2:1312986/1‑93 (MQ=255)
                                     cGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCa                                                    <  1:661139/93‑1 (MQ=255)
                                     cGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCa                                                    >  2:1708145/1‑93 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGGGCTGGGGTATTTCAct                                                  >  2:3909174/1‑92 (MQ=255)
                                        ttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCt                                                 >  2:1911724/1‑93 (MQ=255)
                                         tccatttAGCGCTTCCGCCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTc                                                <  2:159048/89‑1 (MQ=255)
                                         tACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTc                                                <  2:6101238/93‑1 (MQ=255)
                                           ccTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGc                                              <  1:1180883/93‑1 (MQ=255)
                                            cTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCa                                             >  1:5823664/1‑93 (MQ=255)
                                            cTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCa                                             >  2:4381485/1‑93 (MQ=255)
                                            cTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCa                                             <  2:5312350/93‑1 (MQ=255)
                                             tttAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAg                                            >  1:4121250/1‑93 (MQ=255)
                                                   gcCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTaaa                                      <  1:5305381/93‑1 (MQ=255)
                                                   gcCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTaaa                                      >  2:885518/1‑93 (MQ=255)
                                                     cTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGa                                    <  2:10445/93‑1 (MQ=255)
                                                      tCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGat                                   <  1:4272732/93‑1 (MQ=255)
                                                       ccGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGata                                  <  2:2783396/93‑1 (MQ=255)
                                                       ccGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGata                                  <  2:32810/93‑1 (MQ=255)
                                                       ccGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGata                                  <  2:5069712/93‑1 (MQ=255)
                                                       ccGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGata                                  >  1:896485/1‑93 (MQ=255)
                                                         gTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATc                                <  1:4569845/93‑1 (MQ=255)
                                                         gTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATc                                >  1:5374759/1‑93 (MQ=255)
                                                           cGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCaa                              >  1:4781476/1‑93 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGGCTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                             >  2:1281021/1‑93 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                             >  2:4991372/1‑93 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                             <  1:1730858/93‑1 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                             <  1:1302455/93‑1 (MQ=255)
                                                             ccaccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACa                            <  1:2396616/93‑1 (MQ=255)
                                                             ccaccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACa                            <  2:2277829/93‑1 (MQ=255)
                                                              caccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAg                           <  1:3657115/93‑1 (MQ=255)
                                                              caccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAg                           >  2:3719458/1‑93 (MQ=255)
                                                              caccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAg                           >  1:5085921/1‑93 (MQ=255)
                                                               accaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATCTCAACAGc                          >  1:70010/1‑93 (MQ=255)
                                                               accaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGc                          <  2:3780252/93‑1 (MQ=255)
                                                                ccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCg                         >  2:6261899/1‑93 (MQ=255)
                                                                ccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCg                         >  1:6128265/1‑93 (MQ=255)
                                                                 caTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGa                        <  1:4475350/93‑1 (MQ=255)
                                                                 caTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGa                        >  2:2051351/1‑93 (MQ=255)
                                                                 caTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGa                        <  2:1748564/93‑1 (MQ=255)
                                                                  aTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGAc                       >  1:5900642/1‑93 (MQ=255)
                                                                  aTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGAc                       >  2:2666132/1‑93 (MQ=255)
                                                                  aTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGAc                       <  2:2374258/93‑1 (MQ=255)
                                                                  aTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGAc                       <  1:3431606/93‑1 (MQ=255)
                                                                   ttACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACg                      >  2:1530051/1‑93 (MQ=255)
                                                                   ttACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACg                      >  1:1545292/1‑93 (MQ=255)
                                                                      cTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGccc                   >  1:3615133/1‑93 (MQ=255)
                                                                      cTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGccc                   <  2:827251/93‑1 (MQ=255)
                                                                       tGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcccc                  >  2:4021142/1‑93 (MQ=255)
                                                                       tGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcccc                  <  2:115621/93‑1 (MQ=255)
                                                                           ctgctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggc              >  1:2233897/1‑93 (MQ=255)
                                                                           ctgctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggc              >  1:1666255/1‑93 (MQ=255)
                                                                             gctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctg            >  1:2352283/1‑93 (MQ=255)
                                                                               tgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggc          <  1:3365945/93‑1 (MQ=255)
                                                                               tgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggc          <  1:2786860/93‑1 (MQ=255)
                                                                               tgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggc          <  1:5891509/93‑1 (MQ=255)
                                                                                gTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggct         <  1:3129040/93‑1 (MQ=255)
                                                                                gTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggct         <  1:1149531/93‑1 (MQ=255)
                                                                                     aTGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGcgcg    >  2:3195827/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTC  >  NC_000913/329863‑330042

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: