Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,538,491 C→T G92D (GGC→GAC)  nanM ← N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,4910CT100.0% 280.5 / NA 78G92D (GGC→GAC) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (45/33);  total (45/33)

TTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538356‑4538625
                                                                                                                                       |                                                                                                                                      
ttGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:359039/139‑1 (MQ=255)
ttGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:641974/139‑1 (MQ=255)
  gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  1:389731/139‑1 (MQ=255)
  gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:128157/1‑139 (MQ=255)
    gtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  1:1280067/139‑1 (MQ=255)
     tACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAac                                                                                                                                >  1:401463/1‑139 (MQ=255)
        aaaaGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAAcaca                                                                                                                             <  1:366751/139‑1 (MQ=255)
                   ggCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  <  1:723760/139‑1 (MQ=255)
                   ggCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  <  2:1161231/139‑1 (MQ=255)
                       cGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCa                                                                                                                 >  1:617693/1‑136 (MQ=255)
                       cGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCa                                                                                                                 <  2:617693/136‑1 (MQ=255)
                        gCCATGCCCATCGGCGCGCGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTTCCATCaa                                                                                                             >  1:308902/1‑139 (MQ=255)
                         ccATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            <  2:1249651/139‑1 (MQ=255)
                             gCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAAt                                                                                                        >  2:272225/1‑139 (MQ=255)
                              cccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATg                                                                                                       <  2:599643/139‑1 (MQ=255)
                                cATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCa                                                                                                     >  1:1031602/1‑139 (MQ=255)
                                   cGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  1:1033648/108‑1 (MQ=255)
                                   cGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 >  2:1033648/1‑108 (MQ=255)
                                        cgTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGc                                                                                             >  1:1313921/1‑139 (MQ=255)
                                           gCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                                                          <  1:505195/139‑1 (MQ=255)
                                                             aaCTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAgg                                                                        >  2:1154531/1‑139 (MQ=255)
                                                               cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  2:1230673/78‑1 (MQ=255)
                                                               cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:1230673/1‑78 (MQ=255)
                                                                   tGGTTTTGGGGTTGTGTTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:879965/1‑74 (MQ=255)
                                                                   tGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATAGCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  2:879965/74‑1 (MQ=255)
                                                                   tGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGc                                                                  >  2:1162540/1‑139 (MQ=255)
                                                                    ggTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCa                                                                 >  2:1236096/1‑139 (MQ=255)
                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:1038682/1‑68 (MQ=255)
                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  2:1038682/68‑1 (MQ=255)
                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATc                                                                                                               <  1:714377/88‑1 (MQ=255)
                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATc                                                                                                               >  2:714377/1‑88 (MQ=255)
                                                                          ggggTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAAc                                                           >  2:759550/1‑139 (MQ=255)
                                                                          ggggTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAAc                                                           >  2:263852/1‑139 (MQ=255)
                                                                           gggTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACg                                                          <  2:156900/139‑1 (MQ=255)
                                                                                gTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCTATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCt                                                        >  1:552694/1‑136 (MQ=255)
                                                                                gTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCTATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCt                                                        <  2:552694/136‑1 (MQ=255)
                                                                                gTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTc                                                     <  2:713444/139‑1 (MQ=255)
                                                                                      gtgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                                                          >  1:608910/1‑96 (MQ=255)
                                                                                      gtgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                                                          <  2:608910/96‑1 (MQ=255)
                                                                                      gtgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGAGAATGCAGCTAACGCTGTCCAtttt                                               >  2:208022/1‑139 (MQ=255)
                                                                                       tgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAttttt                                              <  1:272225/139‑1 (MQ=255)
                                                                                        gtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAtttttt                                             >  1:1019507/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAttttt                                              <  2:568060/135‑1 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAttttt                                              >  1:568060/1‑135 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  2:546921/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  2:986790/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:1273057/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:9450/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:1048541/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:724952/1‑139 (MQ=255)
                                                                                                ttAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     >  1:385574/1‑139 (MQ=255)
                                                                                                ttAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     <  1:1162540/139‑1 (MQ=255)
                                                                                                 tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGc                                    >  2:1024295/1‑139 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     >  1:297423/1‑35 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:297423/35‑1 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:91928/1‑37 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  2:91928/37‑1 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            <  2:151356/60‑1 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            >  1:151356/1‑60 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAtttt                                               <  1:1158422/121‑1 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAtttt                                               >  2:1158422/1‑121 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTa                             >  1:1192892/1‑139 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaa                                                                                                                                  <  1:677150/36‑1 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaa                                                                                                                                  >  2:677150/1‑36 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATc                           >  2:683722/1‑139 (MQ=255)
                                                                                                           aGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATcc                          >  1:863932/1‑139 (MQ=255)
                                                                                                             tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATc                                                                                                               >  1:417177/1‑52 (MQ=255)
                                                                                                             tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATc                                                                                                               <  2:417177/52‑1 (MQ=255)
                                                                                                               aaGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCt                      >  1:253837/1‑139 (MQ=255)
                                                                                                                aGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCtt                     >  2:583480/1‑139 (MQ=255)
                                                                                                                aGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCtt                     <  2:701536/139‑1 (MQ=255)
                                                                                                                  cccTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGt                   >  1:944796/1‑139 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGa                                                                      <  2:470500/84‑1 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGa                                                                      >  1:470500/1‑84 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    <  2:1012199/86‑1 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    >  1:1012199/1‑86 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     <  2:728335/117‑1 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     >  1:728335/1‑117 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  1:134964/1‑139 (MQ=255)
                                                                                                                        cTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATg             <  2:253837/139‑1 (MQ=255)
                                                                                                                            ttttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGt         <  2:1048541/139‑1 (MQ=255)
                                                                                                                                   aaTGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTgcg  <  2:385574/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538356‑4538625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑