Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,342,193 A→G T95A (ACG→GCG)  melA → alpha‑galactosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,342,1930AG100.0% 441.7 / NA 138T95A (ACG→GCG) melAalpha‑galactosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (54/84);  total (54/84)

Only 100 of 138 total aligned reads displayed.

GGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATG  >  NC_000913/4342048‑4342340
                                                                                                                                                 |                                                                                                                                                   
ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAAACTTTCGCGGTg                                                                                                                                                >  1:903701/1‑151 (MQ=255)
ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGt                                                                                                                                                 <  1:70915/150‑1 (MQ=255)
  aaGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGa                                                                                                                                               <  4:1928422/150‑1 (MQ=255)
  aaGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCgg                                                                                                                                                  >  1:1895894/1‑147 (MQ=255)
    gagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCCTGCGCgggggct                                                                                                                                             >  3:113051/1‑145 (MQ=255)
     agTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTTCGCgggg                                                                                                                                                >  1:1758381/1‑144 (MQ=255)
     agTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGa                                                                                                                                           <  2:125390/151‑1 (MQ=255)
     agTCGCATATTGTGGTGCGTAAGCTGATGGATTAAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGa                                                                                                                                           <  3:161440/151‑1 (MQ=255)
      gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGAt                                                                                                                                          <  1:982097/151‑1 (MQ=255)
        cGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGAttt                                                                                                                                        <  1:461205/151‑1 (MQ=255)
         gCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGAttt                                                                                                                                        <  4:418518/150‑1 (MQ=255)
         gCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTc                                                                                                                                       <  4:1334988/151‑1 (MQ=255)
          cATATTGTGGTGCGTAAGTTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTc                                                                                                                                       <  4:546005/150‑1 (MQ=255)
           atatTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGa                                                                                                                                     <  4:580891/151‑1 (MQ=255)
           atatTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTCgg                                                                                                                                     >  4:1387341/1‑150 (MQ=255)
            tatTGTGGTGCGTAAGCTGATGGATTCGGCAGTGGCCAGGGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTc                                                                                                                                       <  4:1308335/148‑1 (MQ=255)
                 tggtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCt                                                                                                                                <  2:462124/150‑1 (MQ=255)
                  ggtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGt                                                                                                                              <  3:321460/151‑1 (MQ=255)
                  ggtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAACCTTGCGCGGGGGCTGATTTTGAgg                                                                                                                                   >  4:660438/1‑146 (MQ=255)
                   gtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTa                                                                                                                             <  4:672854/151‑1 (MQ=255)
                   gtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTa                                                                                                                             <  1:1255277/151‑1 (MQ=255)
                      cGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGc                                                                                                                          <  1:542809/151‑1 (MQ=255)
                      cGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGACTGATTTCGCGGGCTGt                                                                                                                              >  4:399991/1‑147 (MQ=255)
                       gTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCg                                                                                                                         <  4:806445/151‑1 (MQ=255)
                        tAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCg                                                                                                                         <  1:617672/150‑1 (MQ=255)
                         aaGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCgt                                                                                                                        <  4:214093/150‑2 (MQ=255)
                         aaGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCgg                                                                                                                        <  4:4850/150‑1 (MQ=255)
                                 ggATCCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGATGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTATg                                                                                                               <  4:19221/151‑1 (MQ=255)
                                    ttCAGCAGGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGa                                                                                                             <  4:969597/150‑1 (MQ=255)
                                     tCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGaa                                                                                                            >  2:51934/1‑150 (MQ=255)
                                      cagcagGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGaa                                                                                                            <  1:659563/149‑1 (MQ=255)
                                         cagGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTTGAGGTCTGTAAGCGGCATGGTCTGGAACaaa                                                                                                        >  2:140237/1‑150 (MQ=255)
                                          agGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAACGGCCAGGTCTGGGAAAAAcc                                                                                                      >  4:947584/1‑151 (MQ=255)
                                          agGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGCTATGAAACATGCGCGGTGGCTGGTTTCGAGGGCTGTAAGCGGCATGGTGTGGGACAAAcc                                                                                                      >  3:397526/1‑151 (MQ=255)
                                           ggggCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAAcc                                                                                                      <  2:492914/150‑1 (MQ=255)
                                           ggggCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCa                                                                                                     <  1:1282165/151‑1 (MQ=255)
                                           ggggCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGGTTTTGAGGTCTTTAAGCGGCCTGGTCTGGAACAAAcc                                                                                                      >  3:422400/1‑150 (MQ=255)
                                             ggCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAAACAtt                                                                                                   >  4:1721924/1‑151 (MQ=255)
                                                  gCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCg                                                                                               <  1:933181/150‑1 (MQ=255)
                                                    ggCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCGCGGGGACTGATTTTGAGGTCTGTAAGCGGCATGGTCTGGGACAAAACATTGCCg                                                                                               >  4:768215/1‑148 (MQ=255)
                                                    ggCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGc                                                                                                 <  2:263644/146‑1 (MQ=255)
                                                     gCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAACCTTTCGCGGTGAATGATTTTGAGGTCTGTAAGCGGCATGGGCTGGAACAAACCATTTccccta                                                                                            >  4:1494035/1‑146 (MQ=255)
                                                        aaaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGGCTGATTTCGAGGGCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGtt                                                                                        >  1:1869212/1‑151 (MQ=255)
                                                          aTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTgg                                                                                      <  4:1423115/151‑1 (MQ=255)
                                                          aTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTgg                                                                                      <  1:364575/151‑1 (MQ=255)
                                                            cACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTgggg                                                                                    <  1:243142/151‑1 (MQ=255)
                                                              ccTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGcc                                                                                  <  3:134741/151‑1 (MQ=255)
                                                              ccTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACCTCTGCGGTGGCTGATTTTGGGGTCTGTAAACGGCCTGGGCTGGGACAAACCATTGCCGATACGTTgggg                                                                                    >  1:1751284/1‑149 (MQ=255)
                                                              ccTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCACTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGcc                                                                                  <  2:183604/151‑1 (MQ=255)
                                                                tGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGACTGGTTTTGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTTCCGATACGTTTGGGCCgg                                                                                >  3:46575/1‑151 (MQ=255)
                                                                 gCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCgg                                                                                <  2:380511/150‑1 (MQ=255)
                                                                  ccACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGc                                                                              <  4:90215/151‑1 (MQ=255)
                                                                    acacCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACCTGCGCGGGGACTGGTTTCGTGGTCTGTAAGCGGCATGGGCTGGAAAAAAACATTTCCGATACGTTGGGGGCgggggg                                                                            >  1:477908/1‑151 (MQ=255)
                                                                        ccAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtatt                                                                        <  4:1143133/151‑1 (MQ=255)
                                                                          aaCAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTTCGCGGGGACTGATTTTGAGGGCTGTAAGCGGCATGGGCTGGGACAAACCATTTCCGATACGTTGGGGCCGGGCGGgattat                                                                      >  4:1900704/1‑151 (MQ=255)
                                                                           aCAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATg                                                                     <  4:1464881/151‑1 (MQ=255)
                                                                             agaaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAACCTTGCGCGGGGGCTGATTTCGAGGTCTGTAAGCGGCATGGTGTGGGACAAACCATTGCCGGTACGTTGGGGCCGGGGGGTATTATgcg                                                                   >  1:395279/1‑151 (MQ=255)
                                                                               aaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGTGACTGGTTTTGAGGTCTTTAAGCGGCATGGTCTGGAACAAAACATTTCCGATAAGTTGGGGGCCGGCGGGGTTTTgcg                                                                   >  1:1052465/1‑149 (MQ=255)
                                                                                aagaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgc                                                                <  1:674388/151‑1 (MQ=255)
                                                                                  gaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTGCGCGGGGACTGATTTTGAGGGCTGTAAGCGGCATGGGCTGGAACAAAACATTGCCGATACGTTGGGGCCGGGGGGGATTATtcgcgcg                                                               >  4:800756/1‑150 (MQ=255)
                                                                                  gaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgcgc                                                              <  2:397539/151‑1 (MQ=255)
                                                                                    agCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTa                                                            <  1:907940/151‑1 (MQ=255)
                                                                                    agCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTa                                                            <  1:748329/151‑1 (MQ=255)
                                                                                     gCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGATTGGTTTCGGGGTTTGTAAGAGGGATGGTTTGGGACAAACCCTTGCCGGTACGATGGGGCCGGGCGGGGTTATtcgcgccgtcc                                                           >  4:650246/1‑146 (MQ=255)
                                                                                         tAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTa                                                        <  3:164039/150‑1 (MQ=255)
                                                                                         tAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTa                                                        <  4:1926110/150‑1 (MQ=255)
                                                                                             ggATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCAtt                                                   <  4:302262/151‑1 (MQ=255)
                                                                                             ggATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCAtt                                                   <  2:284097/151‑1 (MQ=255)
                                                                                             ggATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCAtt                                                   <  3:308400/151‑1 (MQ=255)
                                                                                                 gCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTcct                                                <  3:360453/150‑2 (MQ=255)
                                                                                                 gCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCg                                                <  1:1915337/150‑1 (MQ=255)
                                                                                                  ccGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAAACTTGCGCGGGGACTGATTTTGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTTCCGATACGTTGGGGCCGGGCGGGATTATGCGCGCGCCACGTACCATTcccc                                               >  4:908104/1‑148 (MQ=255)
                                                                                                    gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTTAACCTTGTGCGGGGGGTGATTTTGAGGGCTGTAAGCGGCTCGGTCTGGAACAAAACCTTTCCGGTACGTTGGGGCCCGGGGGTAGTATGCGCGCGCTCCGTACCCTTCCGCAt                                             >  1:1593096/1‑150 (MQ=255)
                                                                                                    gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAACCTTGCGCGGTGGCTGGTTTTGGGGGCTGTAAGCGGGATGGGCTGGAACAAAACCTTGCCGATACGGTTGGGCCGGGGGGGATTATGCGCGCGGTACGTACCCTTTCCCATc                                            >  1:60111/1‑151 (MQ=255)
                                                                                                    gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATc                                            <  1:1900923/151‑1 (MQ=255)
                                                                                                        ttGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGGATGAATTCGTGGGTTGTAAGCGGCATGGTCTGGAAGAAAACATTGCCGCTACGTTTGGGCCGGGGCGGATTTTTCGCGGGCCACCGAACAATCCGCATc                                            >  1:1681644/1‑147 (MQ=255)
                                                                                                          gTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGc                                      <  1:292305/151‑1 (MQ=255)
                                                                                                            cGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCa                                     <  1:1107120/150‑1 (MQ=255)
                                                                                                            cGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGACTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCaa                                    <  4:1016638/151‑1 (MQ=255)
                                                                                                             gtggtgGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCaaa                                   <  1:960345/151‑1 (MQ=255)
                                                                                                              tggtggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCaa                                    <  1:441881/149‑1 (MQ=255)
                                                                                                               ggtggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAAtt                                 <  4:34714/151‑1 (MQ=255)
                                                                                                               ggtggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAAt                                  <  1:1350421/150‑1 (MQ=255)
                                                                                                                gtggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAAttt                                <  4:213185/151‑1 (MQ=255)
                                                                                                                 tggCATTTCAGATTGGCGGTTATTAACCTTGCGCGGGGACCGATTTCGAGGGCTGTAAGCGGCATGGTCTGGAACAAACCAATGCCCGTACGTTTGGGCCGGGCGGTATTATGTGCGCGCCACGTACCCTTCCGCATCTTTGGCAAAatt                                >  1:1187390/1‑150 (MQ=255)
                                                                                                                 tggCATTTCAGATTGGCGGTTATGAACCCTTCGCGGGGACTGATTTCGAGGTCTTGTAGCGGCATGGTTTTGAACAAAACAATTCCGATACGTTTGGGGCGGGGGGTATTATGGGGGCGCTAAGGACCATTCCCCATCTTTGGGAAAattt                               >  1:1562332/1‑150 (MQ=255)
                                                                                                                   gCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCg                             <  4:904093/151‑1 (MQ=255)
                                                                                                                   gCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGGTTTTGAGGGCTGTAAGCGGCATGGTCTGGAACAAAACATTGCCGATACGGTGGGGCCGGGCGGGATTATTCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCg                             >  4:45369/1‑151 (MQ=255)
                                                                                                                        tCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGGAAAAAACATTGCCGGTACGTTGGGGCCGGGCGGGATTATGCGCGCGCTACGTACCATTCCGCCTCTGTGGCAAATTTGCGAgg                          >  1:567857/1‑149 (MQ=255)
                                                                                                                         cAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGAc                        <  4:1372641/150‑1 (MQ=255)
                                                                                                                            aTTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATTCGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATg                     <  4:1155818/150‑1 (MQ=255)
                                                                                                                             ttGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGAc                   >  2:501185/1‑151 (MQ=255)
                                                                                                                              tGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACg                  >  1:1472792/1‑151 (MQ=255)
                                                                                                                               ggCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACgg                 <  4:1555216/151‑1 (MQ=255)
                                                                                                                                 cGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGaa               <  1:148399/151‑1 (MQ=255)
                                                                                                                                   gTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGt             <  4:1820930/151‑1 (MQ=255)
                                                                                                                                      aTGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTg          >  4:1711410/1‑151 (MQ=255)
                                                                                                                                         aaCCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGccc       >  2:376712/1‑151 (MQ=255)
                                                                                                                                           ccTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCg     <  3:452074/151‑1 (MQ=255)
                                                                                                                                            cTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCg     >  4:1199988/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATG  >  NC_000913/4342048‑4342340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: