Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,864,410 T→A 6.0% Q488L (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,4100TA 6.0% 241.7 / 5.6 67Q488L (CAG→CTG) ‡smcChromosome partition protein Smc
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-02

TCCAGGCCGGCAAAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGGCG  >  AE015451/4864179‑4864633
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                               
tCCAGGCCGGCAAAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCa                                                                                                                                                                                                                              <  1:961702/235‑1 (MQ=255)
   aGGCCGGCAAAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCCGGt                                                                                                                                                                                                                           >  1:1080653/1‑235 (MQ=255)
     gCCGGCAAAATCCAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCg                                                                                                                                                                                                                         <  1:1015754/235‑1 (MQ=255)
            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)
                                                                 ccccGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTc                                                                                                                                                                                                                          <  1:229024/172‑1 (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                                                                                                                                         <  2:971206/223‑1 (MQ=255)
                                                                                                 ccACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTc                                                                                                                                         >  1:971209/1‑223 (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                                                                                                                                                                           <  1:1080227/172‑1 (MQ=25)
                                                                                                                  gCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCAACAGCAGGCACAGGGCGACCTGCAGCGCCTGGGTGGCCTGCTGCAGTTGCTcgc                                                                                                                                                                           >  2:1080224/1‑172 (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:1086128/236‑1 (MQ=255)
                                                                                                                                                                                                                     aGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGc        <  1:628865/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)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                               
TCCAGGCCGGCAAAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGGCG  >  AE015451/4864179‑4864633

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