Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,346,883 T→A 4.5% F99I (TTC→ATC)  melA → alpha‑galactosidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,346,8830T→A 4.5% 33.2 0.0 ‑0.0309F99I (TTC→ATC) melAalpha‑galactosidase
Reads supporting (aligned to +/- strand):  new base (7/7):  ref base (157/138):  total (164/145)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 328 total aligned reads displayed.

AATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCA  >  NC_007779/4346783‑4346983
                                                                                                    |                                                                                                    
 aTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGGGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGAttt                                                                                                     <  1:4609667/101‑1 (MQ=255)
   cACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCg                                                                                                   >  1:10004929/1‑101 (MQ=255)
     ccTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAg                                                                                                 <  1:8633765/101‑1 (MQ=255)
     ccTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAg                                                                                                 >  1:1267824/1‑101 (MQ=255)
      cTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAgg                                                                                                >  1:3740154/1‑101 (MQ=255)
        gCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTc                                                                                              <  1:7184292/101‑1 (MQ=255)
         ccACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCt                                                                                             <  1:5843250/101‑1 (MQ=255)
          cacaCCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTg                                                                                            <  1:11462477/101‑1 (MQ=255)
           acacCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGt                                                                                           <  1:43370/101‑1 (MQ=255)
           acacCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGt                                                                                           >  1:6575657/1‑101 (MQ=255)
           acacCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGt                                                                                           >  1:2225116/1‑101 (MQ=255)
            cacCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTa                                                                                          <  1:6091736/101‑1 (MQ=255)
             acCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTaa                                                                                         >  1:492524/1‑101 (MQ=255)
               ccAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGc                                                                                       <  1:852510/101‑1 (MQ=255)
               ccAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGc                                                                                       <  1:9658874/101‑1 (MQ=255)
                cAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCg                                                                                      >  1:11077576/1‑101 (MQ=255)
                   cAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCa                                                                                   >  1:10316072/1‑101 (MQ=255)
                    agaaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCAt                                                                                  <  1:4903496/101‑1 (MQ=255)
                    agaaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGACGGTTCTGAACCTTGCACTGTGACTGATTTCGAGGTCTGTAAGCGGcct                                                                                  >  1:11280380/1‑99 (MQ=255)
                     gaaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATATCGAGGTCTGTAAGCGGCATg                                                                                 >  1:5702274/1‑101 (MQ=255)
                       aagaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGt                                                                               <  1:3557195/101‑1 (MQ=255)
                       aagaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATATCGAGGTCTGTAAGCGGCATGGt                                                                               <  1:6861132/101‑1 (MQ=255)
                         gaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCt                                                                             <  1:6259667/101‑1 (MQ=255)
                         gaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCt                                                                             >  1:1684191/1‑101 (MQ=255)
                         gaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCt                                                                             >  1:8207672/1‑101 (MQ=255)
                           agCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGCTCTgg                                                                           <  1:10620281/101‑1 (MQ=255)
                                 agagGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACaaa                                                                     <  1:9349888/101‑1 (MQ=255)
                                     gATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCAt                                                                 <  1:14069498/101‑1 (MQ=255)
                                      aTGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCAtt                                                                >  1:13653238/1‑101 (MQ=255)
                                      aTGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCAtt                                                                <  1:1199626/101‑1 (MQ=255)
                                      aTGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCAtt                                                                >  1:4986/1‑101 (MQ=255)
                                        gCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGc                                                              <  1:6155994/101‑1 (MQ=255)
                                         ccGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGcc                                                             <  1:3624811/101‑1 (MQ=255)
                                          cGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCg                                                            <  1:10780777/101‑1 (MQ=255)
                                          cGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCg                                                            <  1:11869114/101‑1 (MQ=255)
                                          cGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCg                                                            <  1:389716/101‑1 (MQ=255)
                                           gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGa                                                           >  1:12130154/1‑101 (MQ=255)
                                           gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGa                                                           <  1:12998375/101‑1 (MQ=255)
                                            aTTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGAt                                                          >  1:8406162/1‑101 (MQ=255)
                                            aTTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGAt                                                          >  1:5107058/1‑101 (MQ=255)
                                             ttttGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATa                                                         >  1:13208911/1‑101 (MQ=255)
                                             ttttGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATATCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATa                                                         <  1:11764530/101‑1 (MQ=255)
                                               ttGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACg                                                       >  1:419487/1‑101 (MQ=255)
                                                tGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGt                                                      >  1:12376357/1‑101 (MQ=255)
                                                 gTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGtt                                                     >  1:879188/1‑101 (MQ=255)
                                                 gTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCACGGTCTGGAACAAACCATTGCCGATACGtt                                                     >  1:12848109/1‑101 (MQ=255)
                                                    gtggtgGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTggg                                                  <  1:1803980/101‑1 (MQ=255)
                                                       gtggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGcc                                               >  1:11414618/1‑101 (MQ=255)
                                                       gtggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGcc                                               <  1:108142/101‑1 (MQ=255)
                                                        tggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGTCCg                                              <  1:2095452/101‑1 (MQ=255)
                                                        tggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCg                                              <  1:10551061/101‑1 (MQ=255)
                                                        tggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCg                                              <  1:1945471/101‑1 (MQ=255)
                                                         ggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATATCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCgg                                             <  1:7685511/101‑1 (MQ=255)
                                                           cATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGc                                           <  1:9459131/101‑1 (MQ=255)
                                                            aTTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGTGCg                                          >  1:12044633/1‑101 (MQ=255)
                                                             tttCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCgg                                         <  1:8944953/101‑1 (MQ=255)
                                                             tttCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCgg                                         >  1:13944003/1‑101 (MQ=255)
                                                             tttCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCgg                                         >  1:4839555/1‑101 (MQ=255)
                                                               tCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGta                                       >  1:11689942/1‑101 (MQ=255)
                                                               tCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGta                                       >  1:3233796/1‑101 (MQ=255)
                                                                cAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtat                                      >  1:8702479/1‑101 (MQ=255)
                                                                cAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtat                                      >  1:3945267/1‑101 (MQ=255)
                                                                cAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtat                                      <  1:7849376/101‑1 (MQ=255)
                                                                 aGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtatt                                     <  1:14539429/101‑1 (MQ=255)
                                                                 aGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtatt                                     >  1:3824347/1‑101 (MQ=255)
                                                                   aTTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtattat                                   >  1:10504017/1‑101 (MQ=255)
                                                                   aTTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtattat                                   <  1:2474102/101‑1 (MQ=255)
                                                                   aTTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtattat                                   <  1:1039396/101‑1 (MQ=255)
                                                                   aTTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGtattat                                   <  1:1016042/101‑1 (MQ=255)
                                                                   aTTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCCGGCGGgattat                                   >  1:3312465/1‑101 (MQ=255)
                                                                    ttGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATg                                  >  1:915243/1‑101 (MQ=255)
                                                                     tGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgc                                 <  1:7237149/101‑1 (MQ=255)
                                                                      ggCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcg                                >  1:12092681/1‑101 (MQ=255)
                                                                       gCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgc                               <  1:9915536/101‑1 (MQ=255)
                                                                        cGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcg                              <  1:2769922/101‑1 (MQ=255)
                                                                         ggTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgc                             >  1:5062008/1‑101 (MQ=255)
                                                                         ggTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgc                             <  1:13028378/101‑1 (MQ=255)
                                                                           ttATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgcgc                           <  1:12495107/101‑1 (MQ=255)
                                                                           ttATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgcgc                           <  1:14217594/101‑1 (MQ=255)
                                                                             aTGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTa                         <  1:5299453/101‑1 (MQ=255)
                                                                              tGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTAc                        <  1:9738669/101‑1 (MQ=255)
                                                                                aaCCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGt                      >  1:13822605/1‑101 (MQ=255)
                                                                                 aCCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTa                     <  1:1622696/101‑1 (MQ=255)
                                                                                    ttGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCa                  >  1:5477495/1‑101 (MQ=255)
                                                                                    ttGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCa                  >  1:6445268/1‑101 (MQ=255)
                                                                                      gCACGGTGACTGANTTNNNGNNCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCAtt                <  1:261662/101‑1 (MQ=255)
                                                                                       cACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTc               >  1:7137421/1‑101 (MQ=255)
                                                                                        ccGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTcc              <  1:14419942/100‑1 (MQ=255)
                                                                                        aCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTcc              <  1:4542762/101‑1 (MQ=255)
                                                                                            tGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCAt          >  1:11318791/1‑101 (MQ=255)
                                                                                             gACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATc         <  1:7569915/101‑1 (MQ=255)
                                                                                              aCTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCt        <  1:5219458/101‑1 (MQ=255)
                                                                                              aCTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCt        >  1:564860/1‑101 (MQ=255)
                                                                                                tGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCtgt      >  1:703266/1‑101 (MQ=255)
                                                                                                 gATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCtgtg     >  1:3073607/1‑101 (MQ=255)
                                                                                                  aTTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTgg    >  1:13371047/1‑101 (MQ=255)
                                                                                                  aTTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTgg    >  1:14273868/1‑101 (MQ=255)
                                                                                                   tttCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGc   <  1:9221680/101‑1 (MQ=255)
                                                                                                   tttCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGc   >  1:10653062/1‑101 (MQ=255)
                                                                                                    ttCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCa  <  1:3034634/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
AATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCA  >  NC_007779/4346783‑4346983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.