Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 329,174 A→G Y163C (TAC→TGC)  betT → high‑affinity choline transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913329,1740AG99.6% 886.0 / ‑5.1 278Y163C (TAC→TGC) betThigh‑affinity choline transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (147/130);  total (147/131)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

Only 100 of 288 total aligned reads displayed.

AGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329084‑329263
                                                                                          |                                                                                         
  gCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGGGCTGATGGGCATGGCGTTCGGATGCttt                                                                                       >  1:2944391/1‑93 (MQ=255)
  gCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCttt                                                                                       >  1:1801828/1‑93 (MQ=255)
   ccGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTa                                                                                      <  1:2889521/93‑1 (MQ=255)
    cgcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAg                                                                                     >  1:5310405/1‑93 (MQ=255)
    cgcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAg                                                                                     >  1:5220562/1‑93 (MQ=255)
     gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    >  2:1242202/1‑93 (MQ=255)
     gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    >  1:8479632/1‑93 (MQ=255)
     gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    >  1:1698912/1‑93 (MQ=255)
     gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    >  1:2709912/1‑93 (MQ=255)
      cgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCt                                                                                   >  1:3158658/1‑93 (MQ=255)
         gTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATc                                                                                <  1:9085774/93‑1 (MQ=255)
          tCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCg                                                                               <  2:2082098/93‑1 (MQ=255)
          tCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCg                                                                               >  1:1389449/1‑93 (MQ=255)
           cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  2:596861/93‑1 (MQ=255)
           cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  1:4143638/93‑1 (MQ=255)
           cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              >  2:1706191/1‑93 (MQ=255)
               cGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTtata                                                                          <  2:8560888/93‑1 (MQ=255)
               cGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTtata                                                                          >  2:4114134/1‑93 (MQ=255)
                   ggTCTGGACGCTGTTTCACTACGGCTTATCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAttt                                                                      <  2:2649635/93‑1 (MQ=255)
                   ggTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAttt                                                                      >  1:6331225/1‑93 (MQ=255)
                       tGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCg                                                                  <  1:4676055/93‑1 (MQ=255)
                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 >  2:9375938/1‑93 (MQ=255)
                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 <  1:2130761/93‑1 (MQ=255)
                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 >  2:6820399/1‑93 (MQ=255)
                         gACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCt                                                                <  2:3601404/93‑1 (MQ=255)
                           cGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTcc                                                              >  2:7545890/1‑92 (MQ=255)
                           cGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCa                                                              <  2:4918804/93‑1 (MQ=255)
                           cGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCa                                                              <  1:6116432/93‑1 (MQ=255)
                            gCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCAc                                                             >  2:7926212/1‑93 (MQ=255)
                             cTGTTTCTCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCAcc                                                            <  2:3327442/93‑1 (MQ=255)
                             cTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCAcc                                                            <  2:4270486/93‑1 (MQ=255)
                             cTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCAcc                                                            >  1:8061314/1‑93 (MQ=255)
                              tGTTTCGCTACGGCTTACCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCa                                                           <  2:90255/93‑1 (MQ=255)
                               gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCac                                                          >  2:4320401/1‑92 (MQ=255)
                               gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCAt                                                          >  2:6075673/1‑93 (MQ=255)
                                tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  2:8618112/93‑1 (MQ=255)
                                tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  1:1140623/93‑1 (MQ=255)
                                tttCACTACGGCGTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  1:6383525/93‑1 (MQ=255)
                                 ttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATcc                                                        <  1:5278517/93‑1 (MQ=255)
                                 ttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATcc                                                        <  2:5089432/93‑1 (MQ=255)
                                  tCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCg                                                       >  1:5529791/1‑93 (MQ=255)
                                    aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     >  1:7727540/1‑93 (MQ=255)
                                    aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     >  1:1603796/1‑93 (MQ=255)
                                     cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    >  2:2029789/1‑93 (MQ=255)
                                        cGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgc                                                 <  1:9076089/93‑1 (MQ=255)
                                        cGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgc                                                 >  2:9513242/1‑93 (MQ=255)
                                          gCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcgc                                               <  1:3540228/93‑1 (MQ=255)
                                            ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             >  1:2563398/1‑93 (MQ=255)
                                            ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             >  1:7015994/1‑93 (MQ=255)
                                              aaCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTa                                           <  2:7484303/93‑1 (MQ=255)
                                              aaCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTa                                           >  2:2346047/1‑93 (MQ=255)
                                               aCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAc                                          <  1:9229732/93‑1 (MQ=255)
                                               aCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAc                                          <  1:6547012/93‑1 (MQ=255)
                                                ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  2:9268814/93‑1 (MQ=255)
                                                ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  2:3959483/93‑1 (MQ=255)
                                                ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  2:3774167/93‑1 (MQ=255)
                                                ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  2:2829025/93‑1 (MQ=255)
                                                   gCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGa                                      <  2:6405138/93‑1 (MQ=255)
                                                   gCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGa                                      >  1:5659809/1‑93 (MQ=255)
                                                      ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   >  1:4491507/1‑93 (MQ=255)
                                                      ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   <  1:7328661/93‑1 (MQ=255)
                                                       gTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCtt                                  >  2:6509247/1‑93 (MQ=255)
                                                       gTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCtt                                  <  1:391647/93‑1 (MQ=255)
                                                        tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 >  1:8309334/1‑93 (MQ=255)
                                                         cGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCg                                <  1:7966694/93‑1 (MQ=255)
                                                           aTGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGt                              >  2:3126895/1‑93 (MQ=255)
                                                            tGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTa                             <  1:8749979/93‑1 (MQ=255)
                                                            tGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTa                             >  2:9601593/1‑93 (MQ=255)
                                                               aTGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          >  2:4430596/1‑93 (MQ=255)
                                                                tGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACg                         >  1:6592662/1‑93 (MQ=255)
                                                                   gcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAt                      >  2:7160837/1‑93 (MQ=255)
                                                                    cTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAtt                     >  1:5093871/1‑93 (MQ=255)
                                                                    cTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAtt                     <  2:8074284/93‑1 (MQ=255)
                                                                     tGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTa                    >  2:8013583/1‑93 (MQ=255)
                                                                     tGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTa                    >  2:5334118/1‑93 (MQ=255)
                                                                     tGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTa                    >  2:5831840/1‑93 (MQ=255)
                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCGCCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   <  2:1422228/93‑1 (MQ=255)
                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATTCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   >  2:7088516/1‑93 (MQ=255)
                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   <  2:1548341/93‑1 (MQ=255)
                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   >  1:8846004/1‑93 (MQ=255)
                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   >  2:885373/1‑93 (MQ=255)
                                                                        tGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACg                 <  2:4741224/93‑1 (MQ=255)
                                                                        tGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACg                 <  2:9780309/93‑1 (MQ=255)
                                                                         gggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                <  1:9745618/93‑1 (MQ=255)
                                                                         gggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                <  1:4234170/93‑1 (MQ=255)
                                                                            cATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGcc             >  1:10020176/1‑93 (MQ=255)
                                                                             aTGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCg            <  2:2324584/93‑1 (MQ=255)
                                                                              tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           >  2:8227062/1‑93 (MQ=255)
                                                                              tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           >  1:305969/1‑93 (MQ=255)
                                                                              tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           >  1:1454546/1‑93 (MQ=255)
                                                                              tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           <  1:1074138/93‑1 (MQ=255)
                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  1:7154774/93‑1 (MQ=255)
                                                                                  gcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAgg       >  1:5949087/1‑93 (MQ=255)
                                                                                   cTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTCAACGGATTAACGGGCCGATAGGt      >  2:5520347/1‑93 (MQ=255)
                                                                                     cGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCa    <  1:9618624/93‑1 (MQ=255)
                                                                                      ggATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCAc   >  2:10002950/1‑93 (MQ=255)
                                                                                      ggATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCAc   >  1:3398484/1‑93 (MQ=255)
                                                                                       gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  2:240576/1‑93 (MQ=255)
                                                                                       gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  1:4791936/1‑93 (MQ=255)
                                                                                       gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  1:3634015/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329084‑329263

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: