Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,684,470 T→C Y206H (TAC→CAC)  hmpA → dihydropteridine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,4700TC100.0% 452.5 / NA 141Y206H (TAC→CAC) hmpAdihydropteridine reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (81/60);  total (81/60)

Only 100 of 146 total aligned reads displayed.

GACGGTGGCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGTACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGCACAA  >  NC_000913/2684380‑2684558
                                                                                          |                                                                                        
  cGGTGGCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCActc                                                                                      >  1:3230015/1‑93 (MQ=255)
  cGGTGGCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCActc                                                                                      <  1:3170132/93‑1 (MQ=255)
    gTGGCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCtt                                                                                    >  1:2570124/1‑93 (MQ=255)
     tGGCGCAGTGGCAGAATCCCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCttt                                                                                   <  2:70960/93‑1 (MQ=255)
     tGGCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCttt                                                                                   <  2:4095713/93‑1 (MQ=255)
      ggCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTg                                                                                  >  1:4053975/1‑93 (MQ=255)
           aGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCg                                                                             >  1:4428163/1‑93 (MQ=255)
            gTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGc                                                                            <  1:798968/93‑1 (MQ=255)
            gTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGc                                                                            >  1:1454134/1‑93 (MQ=255)
              ggCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCaa                                                                          <  1:4249541/93‑1 (MQ=255)
               gCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCaaa                                                                         >  2:3171793/1‑93 (MQ=255)
                cAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAAc                                                                        <  2:2015441/93‑1 (MQ=255)
                 cgAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAAcc                                                                       <  2:2346320/92‑1 (MQ=255)
                  gAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCg                                                                      >  2:1005971/1‑93 (MQ=255)
                  gAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGGTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTTGCACAccc                                                                      >  1:4468024/1‑92 (MQ=255)
                    aTACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGa                                                                    >  1:1095631/1‑93 (MQ=255)
                    aTACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGa                                                                    <  1:591793/93‑1 (MQ=255)
                     tACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGAt                                                                   <  1:2241662/93‑1 (MQ=255)
                     tACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGAt                                                                   <  2:3073878/93‑1 (MQ=255)
                     tACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGAt                                                                   >  2:1130922/1‑93 (MQ=255)
                     tACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGAt                                                                   >  1:2494552/1‑93 (MQ=255)
                      aCCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATg                                                                  >  1:1694253/1‑93 (MQ=255)
                      aCCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATg                                                                  <  2:3470998/93‑1 (MQ=255)
                       ccGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATgg                                                                 <  2:3823065/93‑1 (MQ=255)
                        cGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGc                                                                >  1:2939793/1‑93 (MQ=255)
                          tCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCaa                                                              >  1:1484371/1‑93 (MQ=255)
                               ggCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCt                                                         <  1:3904429/93‑1 (MQ=255)
                                gCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTa                                                        <  2:5328/93‑1 (MQ=255)
                                gCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTa                                                        <  1:2870324/93‑1 (MQ=255)
                                gCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTa                                                        <  2:1552705/93‑1 (MQ=255)
                                gCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTa                                                        >  2:3487549/1‑93 (MQ=255)
                                 cAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTAt                                                       <  2:3900785/93‑1 (MQ=255)
                                  aaTATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTTGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTAtg                                                      >  1:2450219/1‑92 (MQ=255)
                                  aaTATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATc                                                      >  2:2989500/1‑93 (MQ=255)
                                    tatCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGt                                                    >  2:4087958/1‑93 (MQ=255)
                                    tatCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGt                                                    >  1:4138128/1‑93 (MQ=255)
                                      tctcGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTAt                                                  <  2:1140917/93‑1 (MQ=255)
                                      tctcGGCGTCTGGCTGAAGCCGGAAGGTTTCACACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTAt                                                  <  1:4158379/93‑1 (MQ=255)
                                       ctcGGCGTCTGGCTGAAGCCGGAAGTTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTAtt                                                 <  1:1047451/93‑1 (MQ=255)
                                       ctcGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTAtt                                                 <  2:2857348/93‑1 (MQ=255)
                                       ctcGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTAtt                                                 <  1:3431846/93‑1 (MQ=255)
                                         cGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGc                                               >  1:3627845/1‑93 (MQ=255)
                                             gTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGGGGTg                                           >  2:4178878/1‑93 (MQ=255)
                                             gTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCggag                                           >  1:703621/1‑91 (MQ=255)
                                             gTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTg                                           >  1:2098999/1‑93 (MQ=255)
                                             gTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTg                                           >  1:2541425/1‑93 (MQ=255)
                                              tCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGa                                          >  1:3754848/1‑93 (MQ=255)
                                              tCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGa                                          <  2:3029404/93‑1 (MQ=255)
                                              tCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGa                                          >  1:761317/1‑93 (MQ=255)
                                                tGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGaaa                                        >  1:3774814/1‑93 (MQ=255)
                                                 ggCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGaaaa                                       >  2:167791/1‑92 (MQ=255)
                                                 ggCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAAc                                       >  1:2809778/1‑93 (MQ=255)
                                                   cTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAAcgc                                     <  2:300697/93‑1 (MQ=255)
                                                   cTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAAcgc                                     <  2:4548217/93‑1 (MQ=255)
                                                       aGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAag                                 <  1:1152207/93‑1 (MQ=255)
                                                         ccGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAagag                               <  1:2968763/93‑1 (MQ=255)
                                                         ccGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAagag                               <  2:814829/93‑1 (MQ=255)
                                                         ccGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAagaa                               <  1:1599801/93‑2 (MQ=255)
                                                          cGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAgg                              >  1:4613718/1‑93 (MQ=255)
                                                          cGGAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAgg                              <  1:60775/93‑1 (MQ=255)
                                                           ggAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAggg                             <  2:4311766/93‑1 (MQ=255)
                                                           ggAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAggg                             <  1:1315529/93‑1 (MQ=255)
                                                           ggAAGGTTTCCCACATCAGGAAATTCGTCAGCAATCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAggg                             >  2:4009998/1‑93 (MQ=255)
                                                            gAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGACGAgggg                            >  1:2699692/1‑92 (MQ=255)
                                                            gAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAgggg                            >  1:1055118/1‑92 (MQ=255)
                                                            gAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAgggg                            >  2:4347479/1‑92 (MQ=255)
                                                            gAAGGTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGt                            >  1:4530205/1‑93 (MQ=255)
                                                               ggTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTggg                         >  2:1004064/1‑93 (MQ=255)
                                                               ggTTTCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTggg                         >  1:4533467/1‑93 (MQ=255)
                                                                   tCCCACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAgg                     >  1:774782/1‑93 (MQ=255)
                                                                    cccACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGt                    >  2:680447/1‑93 (MQ=255)
                                                                     ccACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTa                   >  2:2899043/1‑93 (MQ=255)
                                                                     ccACATCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTa                   >  2:1746028/1‑93 (MQ=255)
                                                                       acaTCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATc                 >  2:237580/1‑93 (MQ=255)
                                                                       acaTCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATc                 <  2:2249953/93‑1 (MQ=255)
                                                                       acaTCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATc                 <  2:2069296/93‑1 (MQ=255)
                                                                         aTCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCa               >  2:85092/1‑93 (MQ=255)
                                                                         aTCAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCa               >  1:865009/1‑93 (MQ=255)
                                                                           cAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAAc             <  2:2400792/93‑1 (MQ=255)
                                                                           cAGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAAc             <  1:2943910/93‑1 (MQ=255)
                                                                            aGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACt            >  2:2891119/1‑93 (MQ=255)
                                                                            aGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACt            >  1:1958567/1‑93 (MQ=255)
                                                                            aGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACt            >  2:3839696/1‑93 (MQ=255)
                                                                            aGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACt            >  1:2060100/1‑93 (MQ=255)
                                                                            aGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACt            <  1:436585/93‑1 (MQ=255)
                                                                            aGGAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACt            <  2:432741/93‑1 (MQ=255)
                                                                             ggAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGCACCGCGAAGTGGGTGGGCCGGTATCCAACTg           >  2:1419672/1‑93 (MQ=255)
                                                                             ggAAATTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTg           >  1:4586629/1‑93 (MQ=255)
                                                                               aaaTTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGt         >  2:1435457/1‑93 (MQ=255)
                                                                               aaaTTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGt         >  1:4273351/1‑93 (MQ=255)
                                                                               aaaTTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGt         >  2:3929754/1‑93 (MQ=255)
                                                                                aaTTTGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGTCTATCGTATTGCGGTTAAACGCGAAGAGTGTGGGCAGGTATCCAACTGGtt        >  2:4684564/1‑93 (MQ=255)
                                                                                aaTTCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGtt        <  2:1689426/93‑1 (MQ=255)
                                                                                  ttCGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGc      <  2:2278282/93‑1 (MQ=255)
                                                                                    cGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGcac    >  1:305618/1‑93 (MQ=255)
                                                                                    cGTCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGcaa    >  2:1391722/1‑92 (MQ=255)
                                                                                      tCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGCACaa  <  1:3656576/93‑1 (MQ=255)
                                                                                      tCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGCACaa  >  2:3804427/1‑93 (MQ=255)
                                                                                      tCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGCACaa  >  1:1295030/1‑93 (MQ=255)
                                                                                      tCAGCACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGCACaa  >  1:1148256/1‑93 (MQ=255)
                                                                                          |                                                                                        
GACGGTGGCGCAGTGGCAGAATACCGTCCGGGGCAATATCTCGGCGTCTGGCTGAAGCCGGAAGGTTTCCCACATCAGGAAATTCGTCAGTACTCTTTGACTCGCAAACCGGATGGCAAAGGCTATCGTATTGCGGTGAAACGCGAAGAGGGTGGGCAGGTATCCAACTGGTTGCACAA  >  NC_000913/2684380‑2684558

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