Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,604 T→C G51G (GGT→GGC)  deoD → purine‑nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,6040TC99.7% 1034.4 / ‑5.1 300G51G (GGT→GGC) deoDpurine‑nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (143/156);  total (143/157)
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.91e-01

Only 100 of 314 total aligned reads displayed.

GACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTT  >  NC_000913/4618515‑4618693
                                                                                         |                                                                                         
gACCCGCTGCGTGCGAAGTCTATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTTGGCTTCACCGGCACt                                                                                        >  1:5565510/1‑93 (MQ=255)
gACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACt                                                                                        <  1:282780/93‑1 (MQ=255)
  cccGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTa                                                                                      >  1:4928068/1‑93 (MQ=255)
  cccGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTa                                                                                      >  1:4941204/1‑93 (MQ=255)
   ccGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTAc                                                                                     >  1:6798447/1‑93 (MQ=255)
   ccGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTAc                                                                                     >  2:270748/1‑93 (MQ=255)
    cGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACa                                                                                    <  1:8597531/93‑1 (MQ=255)
     gCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACaa                                                                                   <  2:290938/93‑1 (MQ=255)
     gCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACaa                                                                                   >  1:3506281/1‑93 (MQ=255)
      cTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACaaa                                                                                  <  2:5296533/93‑1 (MQ=255)
       tgcgtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAg                                                                                 <  1:424857/93‑1 (MQ=255)
        gcgtgcgAAGTATATTGCCGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAgg                                                                                >  2:8624525/1‑93 (MQ=255)
         cgtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGc                                                                               <  1:1167880/93‑1 (MQ=255)
          gtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGcc                                                                              <  1:3676673/93‑1 (MQ=255)
           tgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGACCAACGTTCGCGGCATGCTGGGCTTCACCGGCACTTACAAAGGACg                                                                             >  1:3876058/1‑93 (MQ=255)
            gcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACACAGGCCGc                                                                            >  1:8381140/1‑93 (MQ=255)
              gAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaa                                                                          >  2:2361986/1‑93 (MQ=255)
              gAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaa                                                                          >  1:7453746/1‑93 (MQ=255)
               aaGTTTATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaaa                                                                         >  1:2433406/1‑93 (MQ=255)
                aGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaaaa                                                                        >  2:6947354/1‑93 (MQ=255)
                 gTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAAt                                                                       <  2:790674/93‑1 (MQ=255)
                   atatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAAttt                                                                     >  2:6968848/1‑93 (MQ=255)
                   atatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAAttt                                                                     <  2:7118789/93‑1 (MQ=255)
                    tatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTc                                                                    >  2:1212022/1‑93 (MQ=255)
                    tatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTc                                                                    <  1:8140607/93‑1 (MQ=255)
                    tatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTc                                                                    >  2:8807273/1‑93 (MQ=255)
                    tatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTc                                                                    <  1:1260071/93‑1 (MQ=255)
                          tGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                              >  1:1039617/1‑93 (MQ=255)
                          tGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                              <  2:7279475/93‑1 (MQ=255)
                           gAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATg                                                             <  2:3060872/93‑1 (MQ=255)
                           gAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATg                                                             <  1:2904482/93‑1 (MQ=255)
                            aaaCTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATgg                                                            <  1:5568474/93‑1 (MQ=255)
                             aaCTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATggg                                                           >  2:6269873/1‑93 (MQ=255)
                              aCTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGt                                                          <  1:5761600/93‑1 (MQ=255)
                               cTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTc                                                         <  2:3157998/93‑1 (MQ=255)
                                tttCCTTGAAGATGCACGTGAAGTGAACAACGTTCGCGGTATGCTAGGCTTCACCGGCACTTACAAAGGCCGCCAAATTTCCGTAATAGGTcc                                                        <  1:7844143/93‑2 (MQ=255)
                                   ccTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACgg                                                     >  2:1633931/1‑93 (MQ=255)
                                    cTTGAAGATGCCCGTGAAGTAAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGGCGCACAATTTCCGTAATGGGTCACGGt                                                    >  2:7879837/1‑93 (MQ=255)
                                     ttGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTa                                                   >  1:6234711/1‑93 (MQ=255)
                                       gAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATg                                                 <  1:6543561/93‑1 (MQ=255)
                                         aGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATggg                                               <  2:7659199/93‑1 (MQ=255)
                                           aTGCCCGTGAAGTTAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTa                                             <  2:2638750/93‑1 (MQ=255)
                                           aTGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTa                                             <  2:1962573/93‑1 (MQ=255)
                                            tGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTAt                                            >  2:4783824/1‑93 (MQ=255)
                                              cccGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATcc                                          <  1:7327733/93‑1 (MQ=255)
                                              cccGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATcc                                          >  1:8554689/1‑93 (MQ=255)
                                              cccGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATcc                                          >  1:2015126/1‑93 (MQ=255)
                                                cGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCg                                        >  2:6159552/1‑93 (MQ=255)
                                                 gtgaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGt                                       <  1:7169583/93‑1 (MQ=255)
                                                  tgaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTc                                      <  1:3459161/93‑1 (MQ=255)
                                                   gaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTcc                                     >  2:3832586/1‑93 (MQ=255)
                                                   gaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTcc                                     <  1:6089118/93‑1 (MQ=255)
                                                    aagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCt                                    <  1:3970266/93‑1 (MQ=255)
                                                    aagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCt                                    <  2:7212848/93‑1 (MQ=255)
                                                    aagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCt                                    >  2:2155312/1‑93 (MQ=255)
                                                       tgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCt                                 <  2:7454444/93‑1 (MQ=255)
                                                       tgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCt                                 <  1:59071/93‑1 (MQ=255)
                                                        gaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTc                                <  2:4376948/93‑1 (MQ=255)
                                                        gaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTc                                <  2:6760378/93‑1 (MQ=255)
                                                          acaacGTTCGCGGTATGCTGGGCTTCACCGGCTCTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa                              <  1:2983384/93‑1 (MQ=255)
                                                          acaacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa                              <  2:4956872/93‑1 (MQ=255)
                                                          acaacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa                              <  1:2666190/93‑1 (MQ=255)
                                                          acaacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa                              <  2:4286732/93‑1 (MQ=255)
                                                          acaacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa                              <  1:1904331/93‑1 (MQ=255)
                                                           caacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCAt                             <  1:2350471/93‑1 (MQ=255)
                                                            aacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                            <  1:1439976/93‑1 (MQ=255)
                                                            aacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                            >  1:7756757/1‑93 (MQ=255)
                                                            aacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                            <  1:7999293/93‑1 (MQ=255)
                                                            aacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                            >  2:3121635/1‑93 (MQ=255)
                                                            aacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                            <  2:4058972/93‑1 (MQ=255)
                                                             acGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                           >  1:5280223/1‑93 (MQ=255)
                                                             acGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                           >  1:6732074/1‑93 (MQ=255)
                                                               gTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTac                         >  2:3619986/1‑93 (MQ=255)
                                                                ttCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTaca                        >  2:1660882/1‑93 (MQ=255)
                                                                     ggTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaaa                   <  1:5382816/93‑1 (MQ=255)
                                                                     ggTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaaa                   >  2:5498532/1‑93 (MQ=255)
                                                                     ggTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaaa                   <  2:1891845/93‑1 (MQ=255)
                                                                     ggTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaaa                   >  2:8152178/1‑93 (MQ=255)
                                                                     ggTATGCTGAGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTATACCaaa                   <  2:5792860/93‑1 (MQ=255)
                                                                       tATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGa                 <  1:869695/93‑1 (MQ=255)
                                                                       tATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGa                 <  2:4024764/93‑1 (MQ=255)
                                                                       tATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGa                 >  1:5175005/1‑93 (MQ=255)
                                                                       tATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGa                 >  1:6272194/1‑93 (MQ=255)
                                                                        aTGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGaa                <  2:274388/93‑1 (MQ=255)
                                                                         tGCTGGGCTTCACCGGCACTTACAAAGGCCTCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAAc               <  1:1776486/93‑1 (MQ=255)
                                                                            tGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGa            <  1:4405504/93‑1 (MQ=255)
                                                                              ggCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATc          >  2:2409033/1‑93 (MQ=255)
                                                                               gCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCa         >  2:8776222/1‑93 (MQ=255)
                                                                               gCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCa         >  1:4179306/1‑93 (MQ=255)
                                                                                gttCACCGGCACTTACAACGGCCGCAAAATTTCCGTCATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAc        <  1:7568846/92‑1 (MQ=255)
                                                                                cTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCATGTTATGGATACGGCGTCCTCCGACATCTACACCAAAGAACTGATCAc        <  2:6451174/93‑1 (MQ=255)
                                                                                cTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAc        <  1:8782200/93‑1 (MQ=255)
                                                                                cTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAc        >  1:6998691/1‑93 (MQ=255)
                                                                                 ttCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAcc       <  1:4475420/93‑1 (MQ=255)
                                                                                 ttCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAcc       <  1:2198672/93‑1 (MQ=255)
                                                                                 ttCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAcc       >  1:4664644/1‑93 (MQ=255)
                                                                                    cccGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt    <  1:1603747/92‑1 (MQ=255)
                                                                                    aCCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt    >  1:3609751/1‑93 (MQ=255)
                                                                                     ccGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAtt   >  1:312659/1‑93 (MQ=255)
                                                                                      cGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAttt  >  1:8425284/1‑93 (MQ=255)
                                                                                         |                                                                                         
GACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTT  >  NC_000913/4618515‑4618693

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