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,9530AG93.4% 678.6 / 20.8 243T423A (ACC→GCC) betThigh‑affinity choline transport
Reads supporting (aligned to +/- strand):  ref base A (7/9);  new base G (112/115);  total (119/124)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTT  >  NC_000913/329863‑330044
                                                                                          |                                                                                           
gAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGcc                                                                                           >  2:5833414/1‑93 (MQ=255)
 aGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCt                                                                                          <  2:4033922/93‑1 (MQ=255)
 aGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCt                                                                                          <  1:6777201/93‑1 (MQ=255)
  gcgcGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTc                                                                                         <  1:2022383/93‑1 (MQ=255)
  gcgcGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTc                                                                                         >  2:5558500/1‑93 (MQ=255)
  gcgcGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTc                                                                                         <  1:7202208/93‑1 (MQ=255)
   cgcgGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTcg                                                                                        <  2:5722373/93‑1 (MQ=255)
    gcgGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcgg                                                                                       <  2:9565673/93‑1 (MQ=255)
    gcgGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcgg                                                                                       <  2:566481/93‑1 (MQ=255)
      ggCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcggcg                                                                                     >  2:8154358/1‑93 (MQ=255)
       gCTTCTACAGCCTTGTGGCGCCGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTcggcgg                                                                                    <  1:1330986/93‑1 (MQ=255)
         ttCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGAc                                                                                  <  1:7815870/93‑1 (MQ=255)
          tCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACt                                                                                 >  1:3781443/1‑93 (MQ=255)
          tCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACt                                                                                 <  2:2052957/93‑1 (MQ=255)
          tCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACt                                                                                 >  1:4055024/1‑93 (MQ=255)
           cgaCAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTc                                                                                <  2:1945454/91‑1 (MQ=255)
           cTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTc                                                                                >  1:6682629/1‑93 (MQ=255)
            tACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTcc                                                                               >  1:2271454/1‑93 (MQ=255)
            tACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTcc                                                                               >  2:6122981/1‑93 (MQ=255)
              cAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCgg                                                                             <  1:5247516/93‑1 (MQ=255)
               aGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTTGCCTGCTGTTTTATGTTGCCTCGGCGGACTCCggg                                                                            >  2:2035595/1‑93 (MQ=255)
               aGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCggg                                                                            >  1:5969942/1‑93 (MQ=255)
                 gctgctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgc                                                                          <  2:8345320/92‑1 (MQ=255)
                  ctgctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgcg                                                                         >  1:365116/1‑93 (MQ=255)
                   tgctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGgcgc                                                                        >  1:4400474/1‑93 (MQ=255)
                    gctgGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCt                                                                       <  1:6511828/93‑1 (MQ=255)
                         cgcAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGgcgc                                                                  >  1:6916347/1‑93 (MQ=255)
                         cgcAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGc                                                                  <  2:2564467/93‑1 (MQ=255)
                          gcAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCt                                                                 >  1:4460272/1‑93 (MQ=255)
                           cAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTg                                                                <  2:905078/93‑1 (MQ=255)
                             gTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTggg                                                              >  2:6742081/1‑93 (MQ=255)
                               aTCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGa                                                            >  1:792489/1‑93 (MQ=255)
                                 ccGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAt                                                          <  2:2213412/93‑1 (MQ=255)
                                   ggCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAttt                                                        >  1:466027/1‑93 (MQ=255)
                                     cGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCa                                                      >  2:279514/1‑93 (MQ=255)
                                      gTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAc                                                     <  1:5849975/93‑1 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAcc                                                    >  2:9296678/1‑93 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAcc                                                    >  1:9819859/1‑93 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAcc                                                    >  1:9259892/1‑93 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAcc                                                    >  1:1833992/1‑93 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAcc                                                    >  1:5606892/1‑93 (MQ=255)
                                       tttACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCAcc                                                    >  1:8764047/1‑93 (MQ=255)
                                         tACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTc                                                  <  1:6825392/93‑1 (MQ=255)
                                         tACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTc                                                  >  1:7086681/1‑93 (MQ=255)
                                           ccTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGc                                                >  1:9336899/1‑93 (MQ=255)
                                            cTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCa                                               <  1:8449698/93‑1 (MQ=255)
                                            cTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCa                                               >  2:7912940/1‑93 (MQ=255)
                                            cTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCa                                               <  2:8084756/93‑1 (MQ=255)
                                              ttAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGc                                             >  2:652191/1‑93 (MQ=255)
                                                aGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCtt                                           <  2:8944207/93‑1 (MQ=255)
                                                   gcCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTaaa                                        >  2:9113685/1‑93 (MQ=255)
                                                   gcCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTaaa                                        >  2:2186933/1‑93 (MQ=255)
                                                     cTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGa                                      >  1:7552234/1‑93 (MQ=255)
                                                      tCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGat                                     >  2:5351929/1‑93 (MQ=255)
                                                       ccGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGata                                    <  2:8880861/93‑1 (MQ=255)
                                                       ccGTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGata                                    <  2:4880191/93‑1 (MQ=255)
                                                         gTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATc                                  >  2:8512807/1‑93 (MQ=255)
                                                         gTCGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATc                                  <  2:7117984/93‑1 (MQ=255)
                                                           cGCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCaa                                >  1:2754498/1‑93 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                               >  2:8678714/1‑93 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                               >  1:3209785/1‑93 (MQ=255)
                                                            gCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAAc                               <  1:7475691/93‑1 (MQ=255)
                                                             ccaccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACa                              >  2:4533000/1‑93 (MQ=255)
                                                             ccaccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACa                              >  1:5062246/1‑93 (MQ=255)
                                                             ccaccaTTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACa                              >  2:4406495/1‑93 (MQ=255)
                                                              caccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAg                             >  1:3570267/1‑93 (MQ=255)
                                                               accaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGc                            <  1:7979634/93‑1 (MQ=255)
                                                                ccaTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGGGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCg                           >  1:1337274/1‑93 (MQ=255)
                                                                 caTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGa                          >  2:503654/1‑93 (MQ=255)
                                                                  aTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGAc                         <  2:10046444/93‑1 (MQ=255)
                                                                  aTTACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGAc                         >  1:4763169/1‑93 (MQ=255)
                                                                    tACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGc                       >  2:698487/1‑93 (MQ=255)
                                                                    tACTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGc                       <  2:182013/93‑1 (MQ=255)
                                                                    tACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGc                       <  2:3817236/93‑1 (MQ=255)
                                                                     aCTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcc                      <  1:9626002/93‑1 (MQ=255)
                                                                     aCTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcc                      <  1:2504519/93‑1 (MQ=255)
                                                                      cTGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGccc                     >  1:730624/1‑93 (MQ=255)
                                                                       tGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcccc                    >  2:3587121/1‑93 (MQ=255)
                                                                       tGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcccc                    <  1:9425323/93‑1 (MQ=255)
                                                                       tGGCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcccc                    >  2:1502005/1‑93 (MQ=255)
                                                                        ggCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGccccc                   <  1:6187543/93‑1 (MQ=255)
                                                                         gCCTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCg                  <  2:4502035/93‑1 (MQ=255)
                                                                          ccTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCgg                 >  2:2518256/1‑93 (MQ=255)
                                                                           ctgctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggc                <  2:638078/93‑1 (MQ=255)
                                                                           ctgctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggc                >  2:1181760/1‑93 (MQ=255)
                                                                            tgctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggct               >  2:3286098/1‑93 (MQ=255)
                                                                             gctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCGACAGCGACGCCCCCggctg              >  2:3875319/1‑93 (MQ=255)
                                                                             gctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctg              <  2:4150941/93‑1 (MQ=255)
                                                                             gctgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctg              >  1:8951315/1‑93 (MQ=255)
                                                                               tgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggc            <  1:9981376/93‑1 (MQ=255)
                                                                               tgTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggc            <  1:8039213/93‑1 (MQ=255)
                                                                                gTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCggctggcg           >  2:7586467/1‑92 (MQ=255)
                                                                                 ttttATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGgctggctg          >  2:1557094/1‑93 (MQ=255)
                                                                                 ttttATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGgctggctg          >  1:4050992/1‑93 (MQ=255)
                                                                                  tttATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTgc         <  1:1167279/93‑1 (MQ=255)
                                                                                   ttATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTgcg        >  1:3287691/1‑93 (MQ=255)
                                                                                     aTGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGcgcg      <  1:6641050/93‑1 (MQ=255)
                                                                                      tgtgGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAATATATCAACAGCGACGCCCCCGGCTGGCTGCGCGt     >  1:4369080/1‑93 (MQ=255)
                                                                                       gtgGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGtc    >  2:7378729/1‑93 (MQ=255)
                                                                                         gACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGtctt  <  2:361382/93‑1 (MQ=255)
                                                                                          |                                                                                           
GAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTT  >  NC_000913/329863‑330044

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