Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,205,033 T→C A26A (GCA→GCG purH ← phosphoribosylaminoimidazolecarboxamideformyltra nsferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,205,0330TC100.0% 705.7 / NA 221A26A (GCA→GCGpurHphosphoribosylaminoimidazolecarboxamideformyltra nsferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (119/102);  total (119/102)

Only 100 of 234 total aligned reads displayed.

ATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCT  >  NC_000913/4204943‑4205122
                                                                                          |                                                                                         
 tCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGcc                                                                                        <  1:4485569/93‑1 (MQ=255)
 tCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGcc                                                                                        >  2:4031644/1‑93 (MQ=255)
 tCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGcc                                                                                        >  1:852932/1‑93 (MQ=255)
  tgtaaaCTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCt                                                                                       <  1:2107419/90‑1 (MQ=255)
  cGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCt                                                                                       <  1:610501/93‑1 (MQ=255)
   ggAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTg                                                                                      >  1:3885450/1‑93 (MQ=255)
   ggAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTg                                                                                      <  2:2355347/93‑1 (MQ=255)
    gAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTgg                                                                                     >  2:2119161/1‑93 (MQ=255)
     aaaCTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTggg                                                                                    >  2:1460335/1‑93 (MQ=255)
      aaCTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGc                                                                                   >  2:1268070/1‑93 (MQ=255)
       aCTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCg                                                                                  >  1:2946323/1‑93 (MQ=255)
       aCTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCg                                                                                  >  1:2614576/1‑93 (MQ=255)
       aCTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCg                                                                                  >  1:4597780/1‑93 (MQ=255)
        cTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGa                                                                                 >  1:708147/1‑93 (MQ=255)
        cTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGa                                                                                 >  1:1198891/1‑93 (MQ=255)
         ttCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGaa                                                                                >  2:2725741/1‑93 (MQ=255)
         ttCGGTTACCGGCAGACCTTTTTCTGCTAACAGACCGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGac                                                                                >  1:3422895/1‑92 (MQ=255)
           cGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAAtt                                                                              >  1:4291530/1‑93 (MQ=255)
            ggTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTc                                                                             >  2:3956331/1‑93 (MQ=255)
              ttACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTcga                                                                           >  2:1997096/1‑93 (MQ=255)
              ttACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTcga                                                                           <  2:5649/93‑1 (MQ=255)
              ttACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTcga                                                                           >  1:3785986/1‑93 (MQ=255)
                 ccGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACCCCGCGTGCGGAAAGCCCCTGGGCGAATTccacga                                                                        <  2:2249533/93‑1 (MQ=255)
                 ccGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTcgacga                                                                        <  1:3778629/93‑1 (MQ=255)
                 ccGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTcgacga                                                                        <  2:174917/93‑1 (MQ=255)
                  cGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGAt                                                                       <  1:922754/93‑1 (MQ=255)
                   ggCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATa                                                                      >  1:1034007/1‑93 (MQ=255)
                    gCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATAc                                                                     >  1:2549136/1‑93 (MQ=255)
                    gCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATAc                                                                     <  1:1965188/93‑1 (MQ=255)
                    gCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATAc                                                                     <  1:1683333/93‑1 (MQ=255)
                     cAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATAcg                                                                    >  1:1975935/1‑92 (MQ=255)
                      aGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCg                                                                   >  1:3256236/1‑93 (MQ=255)
                       gACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCgg                                                                  >  1:4218585/1‑93 (MQ=255)
                       gACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCgg                                                                  >  2:959960/1‑93 (MQ=255)
                       gACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCgg                                                                  >  1:1617563/1‑93 (MQ=255)
                         ccTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCt                                                                >  2:739454/1‑93 (MQ=255)
                          cTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCtt                                                               <  2:4255643/93‑1 (MQ=255)
                          cTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCtt                                                               >  1:2380511/1‑93 (MQ=255)
                             tttCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGt                                                            <  1:4149387/93‑1 (MQ=255)
                             tttCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGt                                                            <  2:256836/93‑1 (MQ=255)
                                 tGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGa                                                        <  1:2425269/93‑1 (MQ=255)
                                    tAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAac                                                     <  2:4139772/93‑1 (MQ=255)
                                    tAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAac                                                     >  1:1342001/1‑93 (MQ=255)
                                    tAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAac                                                     <  1:4307133/93‑1 (MQ=255)
                                    tAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAac                                                     >  2:1898079/1‑93 (MQ=255)
                                     aaCAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAaca                                                    >  1:3819359/1‑93 (MQ=255)
                                      aCAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAacac                                                   >  1:2711732/1‑93 (MQ=255)
                                       cAGACGGGCAGTGCCCCCCGGAGACAGCAGCCCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACt                                                  <  1:151058/93‑1 (MQ=255)
                                        aGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTg                                                 >  2:2797458/1‑93 (MQ=255)
                                        aGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTg                                                 >  2:1703447/1‑93 (MQ=255)
                                         gACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGa                                                >  2:4388975/1‑93 (MQ=255)
                                           cGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGc                                              <  2:2900693/93‑1 (MQ=255)
                                           cGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGc                                              <  2:320234/93‑1 (MQ=255)
                                              gCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCaga                                           <  2:3239265/93‑1 (MQ=255)
                                              gCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCaga                                           >  1:1780049/1‑93 (MQ=255)
                                              gCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCaga                                           <  1:3345701/93‑1 (MQ=255)
                                                 gTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAgcg                                        <  2:3809495/93‑1 (MQ=255)
                                                 gTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAgcg                                        >  1:4085535/1‑93 (MQ=255)
                                                    cccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCgcgg                                     <  1:2321279/93‑1 (MQ=255)
                                                     ccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGACTTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCgcggc                                    >  1:2010060/1‑93 (MQ=255)
                                                     ccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCgcggc                                    >  2:2551785/1‑93 (MQ=255)
                                                     ccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCgcggc                                    >  2:3309225/1‑93 (MQ=255)
                                                     ccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCgcggc                                    >  1:1327261/1‑93 (MQ=255)
                                                     ccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCgcggc                                    <  2:1844290/93‑1 (MQ=255)
                                                       ccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGcggcgg                                  >  2:392612/1‑93 (MQ=255)
                                                       ccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGcggcgg                                  >  2:4449269/1‑93 (MQ=255)
                                                         tGTAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGAc                                >  1:4532844/1‑93 (MQ=255)
                                                           tAGACAGCAGCTCCACACCTCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTg                              <  2:4098205/93‑1 (MQ=255)
                                                           tAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTg                              >  1:774542/1‑93 (MQ=255)
                                                           tAGACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTg                              >  2:842730/1‑93 (MQ=255)
                                                             gACAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGga                            >  1:2899554/1‑93 (MQ=255)
                                                              aCAGCAGCTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGgac                           <  1:3482423/93‑1 (MQ=255)
                                                                agcagcTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGgacga                         <  1:2678398/93‑1 (MQ=255)
                                                                   agcTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGgacgactt                      <  2:4617200/93‑3 (MQ=255)
                                                                   agcTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgt                      <  2:2201471/93‑1 (MQ=255)
                                                                   agcTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgt                      <  2:3438209/93‑1 (MQ=255)
                                                                    gcTCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgtt                     >  2:3690727/1‑93 (MQ=255)
                                                                      tCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgttgt                   >  1:50051/1‑93 (MQ=255)
                                                                      tCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgttgt                   <  1:4623302/93‑1 (MQ=255)
                                                                      tCCACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgttgt                   >  1:1944264/1‑93 (MQ=255)
                                                                       ccACACCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgttgtt                  <  2:2996162/93‑1 (MQ=255)
                                                                        cacaCCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTTTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGttgttg                 >  1:385792/1‑93 (MQ=255)
                                                                        cacaCCGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGttgttg                 >  1:657192/1‑93 (MQ=255)
                                                                            ccGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATg             <  1:4368832/93‑1 (MQ=255)
                                                                            ccGCGTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATg             <  1:1942565/93‑1 (MQ=255)
                                                                             cgcgTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATgg            <  2:3295912/93‑1 (MQ=255)
                                                                             cgcgTGCGGAAAGCGCCCGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATgg            <  2:1338949/93‑1 (MQ=255)
                                                                              gcgTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGGAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGt           <  2:312507/93‑1 (MQ=255)
                                                                              gcgTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGt           >  2:3512419/1‑93 (MQ=255)
                                                                               cgTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTa          <  2:2541221/93‑1 (MQ=255)
                                                                               cgTGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTa          <  2:449650/93‑1 (MQ=255)
                                                                                 tGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTaaa        >  1:2973042/1‑93 (MQ=255)
                                                                                 tGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTaaa        >  2:4531483/1‑93 (MQ=255)
                                                                                 tGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTaaa        <  1:3679575/93‑1 (MQ=255)
                                                                                 tGCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTaaa        >  1:285430/1‑93 (MQ=255)
                                                                                  gCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAAt       <  1:3486254/93‑1 (MQ=255)
                                                                                  gCGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAAt       >  2:935220/1‑93 (MQ=255)
                                                                                   cGGAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATc      <  1:3152867/93‑1 (MQ=255)
                                                                                    cgAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATcc     <  1:3967773/92‑1 (MQ=255)
                                                                                     gAAAGCGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATccc    >  1:1159005/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCT  >  NC_000913/4204943‑4205122

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