Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,343,031 A→G Q374R (CAG→CGG)  melA → alpha‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,343,0310AG89.0% 1148.9 / 112.1 355Q374R (CAG→CGG) melAalpha‑galactosidase
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (156/159);  total (177/178)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 357 total aligned reads displayed.

AGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAAT  >  NC_000913/4342931‑4343131
                                                                                                    |                                                                                                    
   gtgATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGcc                                                                                                   <  1:8004507/101‑1 (MQ=255)
   gtgATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGcc                                                                                                   >  1:8005141/1‑101 (MQ=255)
     gATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGa                                                                                                 >  1:8005146/1‑101 (MQ=255)
      aTTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGAc                                                                                                >  2:8005149/1‑101 (MQ=255)
      aTTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGAc                                                                                                <  2:8004698/101‑1 (MQ=255)
        ttACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCa                                                                                              >  1:8005151/1‑101 (MQ=255)
        ttACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCa                                                                                              <  2:8004579/101‑1 (MQ=255)
        ttACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCa                                                                                              <  1:8004797/101‑1 (MQ=255)
          aCGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCaaa                                                                                            <  1:8004436/101‑1 (MQ=255)
           cGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAg                                                                                           <  1:8004836/101‑1 (MQ=255)
             gCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTc                                                                                         >  1:8005156/1‑101 (MQ=255)
                 cGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTa                                                                                     <  1:8004702/101‑1 (MQ=255)
                  gTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTc                                                                                         <  1:8004682/96‑1 (MQ=255)
                  gTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTAc                                                                                    >  1:8005163/1‑101 (MQ=255)
                  gTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTAc                                                                                    >  2:8005165/1‑101 (MQ=255)
                    ccGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGc                                                                                  >  2:8005169/1‑101 (MQ=255)
                     cGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCt                                                                                 <  1:8004584/101‑1 (MQ=255)
                     cGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCt                                                                                 <  1:8004518/101‑1 (MQ=255)
                     cGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCt                                                                                 >  1:8005173/1‑101 (MQ=255)
                       tAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTAc                                                                               >  2:8005177/1‑101 (MQ=255)
                       tAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTAc                                                                               >  1:8005179/1‑101 (MQ=255)
                        aaCGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTAcc                                                                              >  2:8005184/1‑101 (MQ=255)
                         aCGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCt                                                                             <  1:8004781/101‑1 (MQ=255)
                         aCGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCt                                                                             >  1:8005185/1‑101 (MQ=255)
                            aTGGTTTGATTGATACCCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCg                                                                          <  2:8004544/101‑1 (MQ=255)
                            aTGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCaa                                                                                             <  2:8004604/82‑1 (MQ=255)
                             tGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTAc                                                                                    <  2:8004687/90‑1 (MQ=255)
                             tGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGc                                                                         >  2:8005191/1‑101 (MQ=255)
                              ggTTTGATTGATACCCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCa                                                                        <  1:8004539/101‑1 (MQ=255)
                              ggTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCa                                                                        >  1:8005194/1‑101 (MQ=255)
                               gTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCAt                                                                       <  1:8004828/101‑1 (MQ=255)
                                tttGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTc                                                                           <  2:8004947/96‑1 (MQ=255)
                                tttGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATc                                                                      >  1:8005198/1‑101 (MQ=255)
                                 ttgattgaTAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCt                                                                     >  2:8005201/1‑101 (MQ=255)
                                  tgattgatAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTg                                                                    <  1:8004807/101‑1 (MQ=255)
                                   gattgatAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTgg                                                                   >  1:8005207/1‑101 (MQ=255)
                                     ttgatAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGgcc                                                                 >  1:8005211/1‑101 (MQ=255)
                                     ttgatAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGgcc                                                                 >  2:8005210/1‑101 (MQ=255)
                                      tgatAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTgg                                                                   <  1:8004735/98‑1 (MQ=255)
                                      tgatAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGgccg                                                                >  1:8005215/1‑101 (MQ=255)
                                         tAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGccc                                                             >  1:8005223/1‑101 (MQ=255)
                                         tAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGccc                                                             >  2:8005222/1‑101 (MQ=255)
                                           aCCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTCCCTTCGCATCTGGCCGCCCTg                                                           <  2:8004433/101‑1 (MQ=255)
                                            ccTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTa                                                                                <  2:8004514/79‑1 (MQ=255)
                                            ccTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGgcc                                                                 <  2:8004497/94‑1 (MQ=255)
                                            ccTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGa                                                          >  2:8005226/1‑101 (MQ=255)
                                             cTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGAt                                                         <  2:8004718/101‑1 (MQ=255)
                                              tGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATg                                                        >  2:8005231/1‑101 (MQ=255)
                                                 cacaAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCaa                                                     <  1:8004297/101‑1 (MQ=255)
                                                 cacaAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCaa                                                     <  1:8004941/101‑1 (MQ=255)
                                                 cacaAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCaa                                                     <  2:8004641/101‑1 (MQ=255)
                                                 cacaAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCaa                                                     >  1:8005235/1‑101 (MQ=255)
                                                   caAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAAc                                                   >  2:8005238/1‑101 (MQ=255)
                                                    aaGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGgccgc                                                               <  2:8004795/88‑1 (MQ=255)
                                                     aGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCa                                                 >  2:8005243/1‑101 (MQ=255)
                                                       gATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAAc                                               >  1:8005244/1‑101 (MQ=255)
                                                        aTGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACa                                              <  2:8004847/101‑1 (MQ=255)
                                                         tGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACAt                                             >  2:8005248/1‑101 (MQ=255)
                                                           ttGCGTGGAATTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCa                                                 <  2:8004746/95‑1 (MQ=255)
                                                           ttGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCa                                           >  1:8005253/1‑101 (MQ=255)
                                                           ttGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCa                                           <  2:8004776/101‑1 (MQ=255)
                                                            tGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCaa                                          <  2:8004902/101‑1 (MQ=255)
                                                            tGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCaa                                          >  2:8005254/1‑101 (MQ=255)
                                                               gTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGt                                       <  2:8004673/101‑1 (MQ=255)
                                                               gTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGt                                       <  2:8004388/101‑1 (MQ=255)
                                                                   aaGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAg                                   >  1:8005261/1‑101 (MQ=255)
                                                                   aaGTAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAg                                   >  2:8005262/1‑101 (MQ=255)
                                                                     gTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGAc                                 <  2:8004765/101‑1 (MQ=255)
                                                                      tAGCCTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACg                                >  2:8005266/1‑101 (MQ=255)
                                                                         ccTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACgctg                             >  1:8005273/1‑101 (MQ=255)
                                                                         ccTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACgctg                             >  2:8005271/1‑101 (MQ=255)
                                                                          cTGTCTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACgctgc                            >  2:8005275/1‑101 (MQ=255)
                                                                             tCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGa                         >  1:8005280/1‑101 (MQ=255)
                                                                              cTGGTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGAc                        >  1:8005284/1‑101 (MQ=255)
                                                                               tGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGAcc                       <  1:8004617/101‑1 (MQ=255)
                                                                                ggTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCg                      <  2:8005040/101‑1 (MQ=255)
                                                                                 gTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGa                     >  1:8005288/1‑101 (MQ=255)
                                                                                 gTTGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGa                     >  2:8005290/1‑101 (MQ=255)
                                                                                  ttGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGaa                    <  2:8004658/101‑1 (MQ=255)
                                                                                  ttGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGaa                    <  1:8004664/101‑1 (MQ=255)
                                                                                   tGATGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTATCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGa                         >  1:8005294/1‑95 (MQ=255)
                                                                                   tGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAg                   <  1:8004747/101‑1 (MQ=255)
                                                                                     aTGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCt                 <  2:8004803/101‑1 (MQ=255)
                                                                                      tGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCt                 <  2:8004475/100‑1 (MQ=255)
                                                                                      tGCTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTa                >  2:8005297/1‑101 (MQ=255)
                                                                                        cTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTAtt              >  1:8005301/1‑101 (MQ=255)
                                                                                        cTAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTAtt              <  1:8004601/101‑1 (MQ=255)
                                                                                         tAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTc             >  2:8005308/1‑101 (MQ=255)
                                                                                         tAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTc             >  1:8005306/1‑101 (MQ=255)
                                                                                         tAATGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTc             <  1:8004866/101‑1 (MQ=255)
                                                                                         tAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTc             <  2:8004754/101‑1 (MQ=255)
                                                                                           aTGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCtt           <  2:8004650/101‑1 (MQ=255)
                                                                                            tGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTa          <  2:8004871/101‑1 (MQ=255)
                                                                                            tGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTa          <  1:8004858/101‑1 (MQ=255)
                                                                                            tGGCATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTa          <  1:8004753/101‑1 (MQ=255)
                                                                                               cATTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACgg       >  1:8005312/1‑101 (MQ=255)
                                                                                                aTTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTa          >  1:8005318/1‑97 (MQ=255)
                                                                                                aTTCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGa      <  1:8005047/101‑1 (MQ=255)
                                                                                                 ttCGGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGaa     <  2:8004726/101‑1 (MQ=255)
                                                                                                    aGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAAt  <  1:8004554/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
AGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAAT  >  NC_000913/4342931‑4343131

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