Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,864,416 T→A 6.1% D486V (GAC→GTC) ‡ smc ← Chromosome partition protein Smc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,864,4160TA 6.1% 232.2 / 5.6 66D486V (GAC→GTC) ‡smcChromosome partition protein Smc
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base A (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

AGGCCGGCAAAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGGCG  >  AE015451/4864182‑4864633
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                         
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                                                                                                                                                                           >  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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                         
AGGCCGGCAAAATCAAGGTCGTTGAAATCGTCCACCAGCACAGCCTGCAGATCGGCGCCAAGCACCGTCTCTACCGCCAGCTCCCAACCCGGCTCCACCCGCAAGCCTTCGGCCAGCCGCGGCTGTTGCTCCAGCCCTTGGCCGTGCAGCCATTGCGCAGCACCTGCGCCCGGCTCCAGGGCGGCCTGCTGCAAGGCTTCCAGCGAGGCCAGGCGCCCACCCAGGCGCTGCAGGTCGCCCTGTGCCTGCTGTTGCGCCTGGGTGGCCTGCTGCAGTTGCTCGCGCGCGCTTTCCAGGCGCTCGATCACCTGCTCTTCGCAGAGCTGCAGTTCTTCCAGCAGCATTTCGCTGCTGGCCAGTTGCTCGGCCAGTTCGAGCATGGCGGCATCCTGCGGGTCGCTGCCCAGTTGCTCGCGCTCCTCGACCAACTTGCGCTGGCGCTCGGCCTGGCG  >  AE015451/4864182‑4864633

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