Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,347,339 C→T 2.8% R251C (CGC→TGC)  melA → alpha‑galactosidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077794,347,3390C→T 2.8% 12.2 ‑0.3 ‑0.6281R251C (CGC→TGC) melAalpha‑galactosidase
Reads supporting (aligned to +/- strand):  new base (4/4):  ref base (107/166):  total (111/170)
Fisher's exact test strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.64e-01

Only 100 of 296 total aligned reads displayed.

CAGTTATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTT  >  NC_007779/4347239‑4347439
                                                                                                    |                                                                                                    
cAGTTATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTgc                                                                                                      <  1:5994739/101‑1 (MQ=255)
cAGTTATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTgc                                                                                                      <  1:5108508/101‑1 (MQ=255)
 aGTTATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTgcg                                                                                                     >  1:8685880/1‑101 (MQ=255)
   ttATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCt                                                                                                   >  1:3906362/1‑101 (MQ=255)
   ttATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCt                                                                                                   <  1:14278029/101‑1 (MQ=255)
     aTGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTAc                                                                                                 <  1:1799846/101‑1 (MQ=255)
     aTGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTAc                                                                                                 <  1:12405984/101‑1 (MQ=255)
       gtgAATCTCTACCCGGAACTGCTGGCGGCTTATTAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGa                                                                                               <  1:1879232/101‑1 (MQ=255)
         gAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGTGCTACGaaa                                                                                             >  1:3396269/1‑101 (MQ=255)
         gAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGaaa                                                                                             >  1:1058352/1‑101 (MQ=255)
         gAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGaaa                                                                                             <  1:11642585/101‑1 (MQ=255)
          aaTCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAAt                                                                                            >  1:13653324/1‑101 (MQ=255)
             ctctACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGtt                                                                                         >  1:669224/1‑101 (MQ=255)
             ctctACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGtt                                                                                         >  1:2802908/1‑101 (MQ=255)
             ctctACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGtt                                                                                         >  1:14039338/1‑101 (MQ=255)
               ctACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCa                                                                                       <  1:590491/101‑1 (MQ=255)
               ctACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCa                                                                                       <  1:8098022/101‑1 (MQ=255)
               ctACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCa                                                                                       <  1:6419051/101‑1 (MQ=255)
                tACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaa                                                                                      >  1:9690918/1‑101 (MQ=255)
                 aCCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaa                                                                                     <  1:4170722/101‑1 (MQ=255)
                 aCCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaa                                                                                     >  1:2088101/1‑101 (MQ=255)
                   ccGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaaaa                                                                                   <  1:13153401/101‑1 (MQ=255)
                    cGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAg                                                                                  >  1:7497654/1‑101 (MQ=255)
                    cGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAg                                                                                  <  1:3181062/101‑1 (MQ=255)
                    cGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAg                                                                                  <  1:1739274/101‑1 (MQ=255)
                     ggAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGc                                                                                 >  1:6926567/1‑101 (MQ=255)
                     ggAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGc                                                                                 >  1:3016897/1‑101 (MQ=255)
                      gAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGTGCTACGAAATGTTCAAAAAGCt                                                                                <  1:14628175/101‑1 (MQ=255)
                      gAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCt                                                                                <  1:5412187/101‑1 (MQ=255)
                      gAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCt                                                                                <  1:12750821/101‑1 (MQ=255)
                      gAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCt                                                                                <  1:10444651/101‑1 (MQ=255)
                         ctgctgGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTggg                                                                             >  1:242437/1‑101 (MQ=255)
                          tgctgGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGc                                                                            <  1:945031/101‑1 (MQ=255)
                           gctgGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCt                                                                           >  1:2319910/1‑101 (MQ=255)
                           gcaggcggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCt                                                                           >  1:10363794/4‑101 (MQ=255)
                            ctgGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTa                                                                          <  1:913448/101‑1 (MQ=255)
                              ggcggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTAtt                                                                        <  1:3044709/101‑1 (MQ=255)
                                cggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTc                                                                      <  1:4609327/101‑1 (MQ=255)
                                 ggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCg                                                                     >  1:13429966/1‑101 (MQ=255)
                                 ggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCg                                                                     >  1:7956052/1‑101 (MQ=255)
                                 ggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCg                                                                     <  1:5729708/101‑1 (MQ=255)
                                     tATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGTGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCAc                                                                 <  1:13875976/101‑1 (MQ=255)
                                      aTGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACg                                                                <  1:4679522/101‑1 (MQ=255)
                                      aTGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACg                                                                >  1:6864955/1‑101 (MQ=255)
                                        gAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGa                                                              >  1:8054649/1‑101 (MQ=255)
                                        gAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGa                                                              <  1:8019024/101‑1 (MQ=255)
                                        gAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGa                                                              >  1:11432622/1‑101 (MQ=255)
                                        gAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGa                                                              <  1:6546064/101‑1 (MQ=255)
                                           gcagggcaggCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATc                                                           <  1:9358050/101‑1 (MQ=255)
                                             agggcaggCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGt                                                         >  1:11951977/1‑101 (MQ=255)
                                                 caggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGa                                                     <  1:8964425/101‑1 (MQ=255)
                                                  aggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGaa                                                    <  1:10178455/101‑1 (MQ=255)
                                                  aggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGaa                                                    <  1:1471154/101‑1 (MQ=255)
                                                   ggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAAc                                                   <  1:12248010/101‑1 (MQ=255)
                                                    gcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACa                                                  <  1:7669167/101‑1 (MQ=255)
                                                      aCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAtt                                                <  1:4915508/101‑1 (MQ=255)
                                                      aCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAtt                                                <  1:4365432/101‑1 (MQ=255)
                                                       ccGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAttt                                               <  1:8360676/101‑1 (MQ=255)
                                                        cGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAtttt                                              <  1:10680623/101‑1 (MQ=255)
                                                        cGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAtttt                                              >  1:8581501/1‑101 (MQ=255)
                                                            gCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTg                                          <  1:10911466/101‑1 (MQ=255)
                                                            gCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTg                                          >  1:13285110/1‑101 (MQ=255)
                                                              cGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAg                                        <  1:5043160/101‑1 (MQ=255)
                                                               gAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGt                                       >  1:11269473/1‑101 (MQ=255)
                                                                 atatTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTac                                     <  1:4757539/101‑1 (MQ=255)
                                                                  tatTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTaca                                    >  1:8845179/1‑101 (MQ=255)
                                                                  tatTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTaca                                    >  1:4142840/1‑101 (MQ=255)
                                                                   atTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTacac                                   >  1:3580826/1‑101 (MQ=255)
                                                                   atTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTacac                                   <  1:5807005/101‑1 (MQ=255)
                                                                      cATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACAcc                                <  1:5505411/101‑1 (MQ=255)
                                                                        tGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGt                              <  1:1554376/101‑1 (MQ=255)
                                                                         ggCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTg                             <  1:12669752/101‑1 (MQ=255)
                                                                           cAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGt                           >  1:7294349/1‑101 (MQ=255)
                                                                            aaTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGtt                          <  1:4989994/101‑1 (MQ=255)
                                                                             aTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGttt                         >  1:2610145/1‑101 (MQ=255)
                                                                              tACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTtta                        <  1:14162858/101‑1 (MQ=255)
                                                                                 tCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTttatta                     <  1:6713703/101‑1 (MQ=255)
                                                                                  cGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTaa                    >  1:981859/1‑101 (MQ=255)
                                                                                  cGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTaa                    <  1:13959778/101‑1 (MQ=255)
                                                                                     tGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGcc                 <  1:12018064/101‑1 (MQ=255)
                                                                                      tccAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCa                <  1:13823507/100‑1 (MQ=255)
                                                                                      gCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCa                <  1:1189880/101‑1 (MQ=255)
                                                                                       ccAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAg               <  1:2195916/101‑1 (MQ=255)
                                                                                       ccAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCNATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAg               <  1:11820278/101‑1 (MQ=255)
                                                                                        cAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAgg              <  1:5067295/101‑1 (MQ=255)
                                                                                        cAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAgg              >  1:11743436/1‑101 (MQ=255)
                                                                                        cAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAgg              >  1:11027329/1‑101 (MQ=255)
                                                                                         aGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGt             >  1:3774626/1‑101 (MQ=255)
                                                                                          gAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTc            <  1:13566598/101‑1 (MQ=255)
                                                                                          gAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTc            <  1:12986562/101‑1 (MQ=255)
                                                                                          gAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTc            <  1:12813759/101‑1 (MQ=255)
                                                                                           aaTATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCg           >  1:3955381/1‑101 (MQ=255)
                                                                                            atatTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGt          >  1:8407763/1‑101 (MQ=255)
                                                                                             tatTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTg         >  1:13078562/1‑101 (MQ=255)
                                                                                                tgtgCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAgg      <  1:1005116/101‑1 (MQ=255)
                                                                                                tgtgCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAgg      >  1:13216431/1‑101 (MQ=255)
                                                                                                tgtgCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAgg      >  1:5613019/1‑101 (MQ=255)
                                                                                                  tgCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGAt    >  1:2027338/1‑101 (MQ=255)
                                                                                                  tgCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGAt    <  1:6503894/101‑1 (MQ=255)
                                                                                                   gcgcTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGAtt   <  1:7454894/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CAGTTATGTGAATCTCTACCCGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTT  >  NC_007779/4347239‑4347439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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.