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.3% 886.7 / ‑5.0 285Y163C (TAC→TGC) betThigh‑affinity choline transport
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (122/161);  total (123/162)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 291 total aligned reads displayed.

AGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329084‑329263
                                                                                          |                                                                                         
aGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCt                                                                                         >  1:4357248/1‑93 (MQ=255)
aGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCt                                                                                         <  2:4974829/93‑1 (MQ=255)
aGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCt                                                                                         <  2:2169405/93‑1 (MQ=255)
 ggCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        >  1:6876619/1‑93 (MQ=255)
   ccGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTa                                                                                      >  1:2409412/1‑93 (MQ=255)
   ccGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTa                                                                                      <  1:5206674/93‑1 (MQ=255)
     gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    <  2:6130914/93‑1 (MQ=255)
     gcgcgTCAGGCGATGGACTGGACGCTGTTTCACTACGGCTTCACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGCTGCTTTAGc                                                                                    >  2:5219606/1‑93 (MQ=255)
        cgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTAt                                                                                 <  2:1239570/93‑1 (MQ=255)
         gTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        >  2:5439089/1‑85 (MQ=255)
         gTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATc                                                                                <  1:1262358/93‑1 (MQ=255)
           cAGGCGATGGTCTGGGCGCTGTTTCCCTACGGCTTAACCGCCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  1:2509441/93‑1 (MQ=255)
           cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  2:3758626/93‑1 (MQ=255)
           cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  2:3343394/93‑1 (MQ=255)
           cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  1:666105/93‑1 (MQ=255)
            aGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGtt                                                                             >  2:592419/1‑93 (MQ=255)
            aGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGtt                                                                             >  2:1649003/1‑93 (MQ=255)
             ggCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTta                                                                            <  2:1411800/93‑1 (MQ=255)
             ggCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTta                                                                            >  1:1174672/1‑93 (MQ=255)
             ggCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTta                                                                            >  1:4311159/1‑93 (MQ=255)
             ggCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTta                                                                            <  1:3265860/93‑1 (MQ=255)
              gCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTtat                                                                           <  1:3975471/93‑1 (MQ=255)
               cGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTtata                                                                          <  1:5939458/93‑1 (MQ=255)
               cGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTtata                                                                          <  1:548520/93‑1 (MQ=255)
                gATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATaa                                                                         >  2:576528/1‑93 (MQ=255)
                gATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATaa                                                                         >  2:2747865/1‑93 (MQ=255)
                 aTGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAt                                                                        <  2:4626536/93‑1 (MQ=255)
                  tGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAtt                                                                       >  1:2719450/1‑93 (MQ=255)
                   ggTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAttt                                                                      <  1:1315991/93‑1 (MQ=255)
                     tCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGc                                                                    <  1:1836679/93‑1 (MQ=255)
                     tCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGc                                                                    <  1:6814158/93‑1 (MQ=255)
                     tCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGc                                                                    <  1:6725761/93‑1 (MQ=255)
                      cTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGcc                                                                   <  2:4361702/93‑1 (MQ=255)
                      cTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGcc                                                                   <  1:7758155/93‑1 (MQ=255)
                      cTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGcc                                                                   <  1:2553519/93‑1 (MQ=255)
                       tGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCg                                                                  >  2:2504913/1‑93 (MQ=255)
                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 <  2:1544009/93‑1 (MQ=255)
                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 <  2:6641432/93‑1 (MQ=255)
                         gACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCt                                                                >  2:6433312/1‑93 (MQ=255)
                         gACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGACGCt                                                                >  2:2646421/1‑93 (MQ=255)
                           cGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCa                                                              <  2:253923/93‑1 (MQ=255)
                               gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCAt                                                          >  1:668122/1‑93 (MQ=255)
                               gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCAt                                                          <  1:1739844/93‑1 (MQ=255)
                                tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  2:3107561/93‑1 (MQ=255)
                                tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  2:5671148/93‑1 (MQ=255)
                                tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  1:6169841/93‑1 (MQ=255)
                                    aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     <  1:7050230/93‑1 (MQ=255)
                                    aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     <  1:6445151/93‑1 (MQ=255)
                                    aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     >  1:2302605/1‑93 (MQ=255)
                                      tACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCg                                                   <  1:6299694/93‑1 (MQ=255)
                                          gCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcgc                                               >  2:6138496/1‑93 (MQ=255)
                                           cTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCt                                              <  2:7402885/93‑1 (MQ=255)
                                            ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             >  2:1343494/1‑93 (MQ=255)
                                            ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             <  1:2238533/93‑1 (MQ=255)
                                            ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             <  1:555562/93‑1 (MQ=255)
                                             tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGt                                            >  1:7373715/1‑93 (MQ=255)
                                             tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGt                                            >  2:1945094/1‑93 (MQ=255)
                                              aaCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTa                                           <  1:107758/93‑1 (MQ=255)
                                                 cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAccc                                        <  2:6350375/93‑1 (MQ=255)
                                                 cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAccc                                        >  2:44988/1‑93 (MQ=255)
                                                   gCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGa                                      <  2:6042711/93‑1 (MQ=255)
                                                    gtGGTCGATGTATGCGCTGATGGGCATGGCGCAAGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGAt                                     <  2:6801162/92‑1 (MQ=255)
                                                      ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   >  1:3664532/1‑93 (MQ=255)
                                                       gTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCtt                                  <  2:6907547/93‑1 (MQ=255)
                                                       gTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCtt                                  >  1:7137604/1‑93 (MQ=255)
                                                        tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 <  1:3548302/93‑1 (MQ=255)
                                                        tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 >  2:1751418/1‑93 (MQ=255)
                                                        tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 <  2:3517508/93‑1 (MQ=255)
                                                          gATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCgg                               >  2:1087896/1‑93 (MQ=255)
                                                          gATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGCACCCGATCTTCgg                               >  2:622219/1‑93 (MQ=255)
                                                           aTGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGATGCTCACCATCCGCTCGGCGCTGTACCCGATATTCGGt                              >  2:4114424/1‑93 (MQ=255)
                                                            tGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCCCCATCCGCTCGGCGCTGTACCCGATCTTCGGTa                             <  2:4760008/93‑1 (MQ=255)
                                                            tGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTa                             >  1:21071/1‑93 (MQ=255)
                                                             gTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTaa                            >  2:7467915/1‑93 (MQ=255)
                                                             gTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTaa                            <  1:4774437/93‑1 (MQ=255)
                                                             gTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTaa                            >  1:3736177/1‑93 (MQ=255)
                                                                tGCGCTGATGGGCATGGCGCTCGGATGCTGTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACg                         >  2:4839366/1‑93 (MQ=255)
                                                                  cgcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGa                       <  1:146540/93‑1 (MQ=255)
                                                                   gcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAt                      >  1:7696442/1‑93 (MQ=255)
                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   <  1:7016309/93‑1 (MQ=255)
                                                                       aTGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAAc                  >  2:3428718/1‑93 (MQ=255)
                                                                       aTGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAAc                  <  1:502826/93‑1 (MQ=255)
                                                                        tGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACg                 >  2:6101264/1‑93 (MQ=255)
                                                                        tGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACg                 >  2:116029/1‑93 (MQ=255)
                                                                           gCAGGGCGCTCGGATGCTTTAGCTATCGTTATAATGTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGc              <  1:3091491/93‑1 (MQ=255)
                                                                            cATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGcc             <  2:2833579/93‑1 (MQ=255)
                                                                             aTGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCg            >  2:7758852/1‑93 (MQ=255)
                                                                              tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           >  1:1500146/1‑93 (MQ=255)
                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  2:4072892/93‑1 (MQ=255)
                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  1:867076/93‑1 (MQ=255)
                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  1:4403122/93‑1 (MQ=255)
                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  1:1653052/93‑1 (MQ=255)
                                                                                 cgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAg        <  1:5352701/93‑1 (MQ=255)
                                                                                  gcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAgg       >  2:2361493/1‑93 (MQ=255)
                                                                                  gcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAgg       >  1:2962677/1‑93 (MQ=255)
                                                                                  gcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAgg       >  2:6988271/1‑93 (MQ=255)
                                                                                    tCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     >  1:4702654/1‑93 (MQ=255)
                                                                                      ggATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCAc   <  1:1810298/93‑1 (MQ=255)
                                                                                       gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  2:289655/1‑93 (MQ=255)
                                                                                       gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  1:320883/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329084‑329263

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