Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,864,422 T→A 6.0% Q484L (CAG→CTG)  smc ← Chromosome partition protein Smc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,864,4220TA 6.0% 252.1 / 5.4 67Q484L (CAG→CTG) smcChromosome partition protein Smc
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base A (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

AAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGGCG  >  AE015451/4864191‑4864633
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                   
aaaTCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGc                                                                                                                                                                                                                 >  1:599414/1‑236 (MQ=255)
    cAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctgc                                                                                                                                                                                                             <  2:599411/236‑1 (MQ=255)
                                aGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctgctg                                                                                                                                                                                                           <  1:22047/210‑1 (MQ=255)
                                 gCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctg                                                                                                                                                                                                              <  1:679738/206‑1 (MQ=255)
                                     gCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTg                                                                                                                                                                                                        >  2:420320/1‑208 (MQ=255)
                                         aTCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGt                                                                                                                                                                                               <  2:911386/213‑1 (MQ=255)
                                          tCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTcgcgcgc                                                                                                                                                                       >  1:296654/1‑236 (MQ=255)
                                                      aCCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCCCCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGCCGCCCACCCACGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGc                                                                                                                                                                                            <  1:664579/203‑1 (MQ=255)
                                                                 cGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCt                                                                                                                                                <  1:579367/236‑1 (MQ=255)
                                                                 cGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCt                                                                                                                                                <  1:469913/236‑1 (MQ=255)
                                                                  gCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATccc                                                                                                                                                  >  2:469910/1‑231 (MQ=255)
                                                                        tCCCACCCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTgggcgcctgg                                                                                                                                                                                     <  1:1273358/192‑9 (MQ=255)
                                                                         cccAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCg                                                                                                                                        >  1:146038/1‑236 (MQ=255)
                                                                         cccAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCACGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCg                                                                                                                                        >  1:662117/1‑236 (MQ=255)
                                                                               ccGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCg                                                                                                                                        >  2:673883/1‑230 (MQ=255)
                                                                                     ccACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCctgctg                                                                                                                                                                                      >  2:1166055/1‑178 (MQ=255)
                                                                                     ccACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTc                                                                                                                                         >  1:971209/1‑223 (MQ=255)
                                                                                     ccACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTc                                                                                                                                         <  2:971206/223‑1 (MQ=255)
                                                                                     ccACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGt                                                                                                                            >  1:206726/1‑236 (MQ=255)
                                                                                      cACCCGCAGGCCTGCGGCCAGCCGCGGCTGGTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctgc                                                                                                                                                                                                             <  2:301124/154‑1 (MQ=255)
                                                                                      cACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGcctgtctc                                                                                                                           >  1:66144/1‑229 (MQ=255)
                                                                                             aaGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCcag                                                                                                                    <  1:1224930/236‑1 (MQ=255)
                                                                                                   tCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCAt                                                                                                              >  1:825566/1‑236 (MQ=255)
                                                                                                     ggCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCAACAGCAGGCACAGGGCGACCTGCAGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCAtt                                                                                                             <  1:762023/235‑1 (MQ=255)
                                                                                                      gCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCAACAGCAGGCACAGGGCGACCTGCAGCGCCTGGGTGGCCTGCTGCAGTTGCTcgc                                                                                                                                                                           >  2:1080224/1‑172 (MQ=25)
                                                                                                      gCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCAACAGCAGGCACAGGGCGACCTGCAGCGCCTGGGTGGCCTGCTGCAGTTGCTcgc                                                                                                                                                                           <  1:1080227/172‑1 (MQ=25)
                                                                                                              cgGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTgg                                                                                                   <  1:627174/236‑1 (MQ=255)
                                                                                                                 cTGTTGCGCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctgc                                                                                                                                                                                                             <  2:736704/127‑1 (MQ=255)
                                                                                                                        tCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTc                                                                                         <  1:342590/236‑1 (MQ=255)
                                                                                                                                   gCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGgc                                                                                                                                                            <  1:179474/158‑1 (MQ=255)
                                                                                                                                    ccGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGCGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTccg                                                                                                                                                               >  2:299958/1‑153 (MQ=255)
                                                                                                                                    ccGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGa                                                                              <  1:968881/235‑1 (MQ=255)
                                                                                                                                          aGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGTGCTGCa                                                                                                                              >  2:932884/1‑181 (MQ=255)
                                                                                                                                          aGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTccc                                                                                                                     >  1:1000480/1‑189 (MQ=255)
                                                                                                                                          aGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCg                                                                               >  2:885242/1‑228 (MQ=255)
                                                                                                                                          aGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATgg                                                                        >  1:770933/1‑235 (MQ=255)
                                                                                                                                          aGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATgg                                                                        >  1:969145/1‑235 (MQ=255)
                                                                                                                                          aGCCATTGCGCAGCACCTGCGCCCGACTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCAACAGCAGGCACAGGGCGACCTGCAGCGCCTGGGTGGCCTGCTGCAGtcg                                                                                                                                                                                >  2:63933/1‑129 (MQ=18)
                                                                                                                                            ccATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGctct                                                                                                                                           <  1:60364/166‑1 (MQ=255)
                                                                                                                                              aTTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctg                                                                                                                                                                                                              <  2:76074/97‑1 (MQ=255)
                                                                                                                                                  cgcAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGCGGCCAGGGGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGCCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCt                                                               <  1:762087/236‑1 (MQ=255)
                                                                                                                                                  cgcAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCa                                                                           <  2:825563/224‑1 (MQ=255)
                                                                                                                                                  cgcAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATcc                                                                <  1:1078926/235‑1 (MQ=255)
                                                                                                                                                       caCCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCgg                                                           >  1:1089605/1‑235 (MQ=255)
                                                                                                                                                          cTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTc                                                                                         <  2:296653/202‑1 (MQ=255)
                                                                                                                                                           tGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCg                                                       >  1:669682/1‑235 (MQ=255)
                                                                                                                                                           tGCGCCCGGCGCCAGGGCGGCCTGCTGCACGGCTTGCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCt                                                                                                                                                                                                   <  2:272925/95‑1 (MQ=255)
                                                                                                                                                             cgcCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCctg                                                                                                                                                                                                              <  2:279769/82‑1 (MQ=255)
                                                                                                                                                                ccGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCt                                                                                                                                                <  2:66144/141‑1 (MQ=255)
                                                                                                                                                                                ctgctgCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGctcct                                 <  1:480018/236‑1 (MQ=255)
                                                                                                                                                                                 tgctgcAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGTCAGTTGCTCTGCCAGTTCGAg                                                                             >  1:241097/1‑191 (MQ=255)
                                                                                                                                                                                  gctgcAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAg                                                                             >  1:936030/1‑190 (MQ=255)
                                                                                                                                                                                             accAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTgc                     <  1:694797/234‑1 (MQ=255)
                                                                                                                                                                                                 gCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTgg                >  1:1143282/1‑236 (MQ=255)
                                                                                                                                                                                                 gCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTg                 <  1:588011/235‑1 (MQ=255)
                                                                                                                                                                                                       ccAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCg          <  1:1186620/236‑1 (MQ=255)
                                                                                                                                                                                                        cAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTcgcg                                      >  2:935266/1‑207 (MQ=255)
                                                                                                                                                                                                         aGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGc        <  1:628865/236‑1 (MQ=255)
                                                                                                                                                                                                         aGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGc        <  1:1086128/236‑1 (MQ=255)
                                                                                                                                                                                                               ccacccaGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGgc   <  1:1270004/235‑1 (MQ=255)
                                                                                                                                                                                                                cacccaGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGgcg  >  1:301459/1‑235 (MQ=255)
                                                                                                                                                                                                                         gcTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTg                                                              >  2:521998/1‑166 (MQ=255)
                                                                                                                                                                                                                         gcTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCCCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAg                                                                                               >  2:1233902/1‑133 (MQ=255)
                                                                                                                                                                                                                                gTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCGCCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGGGCATGGCGGCATCCTGCggga                                                         >  2:113653/1‑163 (MQ=255)
                                                                                                                                                                                                                                    cccTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGCTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGcc                                                  >  2:720089/1‑167 (MQ=255)
                                                                                                                                                                                                                                     ccTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTcgcg                                      >  1:979551/1‑178 (MQ=255)
                                                                                                                                                                                                                                     ccTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTcgcg                                      <  2:979548/178‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                   
AAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGGCG  >  AE015451/4864191‑4864633

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