Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,256,285 T→C 8.3% intergenic (+345/‑359) asd → / → PP_1990 Aspartate‑semialdehyde dehydrogenase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,256,2850TC 8.3% 231.6 / 10.3 72intergenic (+345/‑359)asd/PP_1990Aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base T (40/26);  new base C (3/3);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-02

TTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGACACA  >  AE015451/2256053‑2256510
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                 
tttAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                                                >  1:528958/1‑236 (MQ=255)
tttAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                                                <  1:239959/236‑1 (MQ=255)
 ttAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                                                >  1:1453553/1‑235 (MQ=255)
      ggCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGcgcagcgcag                                                                                                                                                                                                                            <  1:566705/235‑1 (MQ=17)
             tGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCgg                                                                                                                                                                                                                     >  2:1038667/1‑234 (MQ=255)
             tGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgat                                                                                                                                                                                                                   >  1:521632/1‑236 (MQ=255)
                cTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                                                <  1:27013/220‑1 (MQ=255)
                                     cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGcgcag                                                                                                                                                                                                      <  2:1495163/225‑5 (MQ=255)
                                     cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGGTTTTGGtt                                                                                                                                                                                           >  1:429861/1‑236 (MQ=255)
                                             ttgcttttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGGCCCCGGGGCGTAGCGCCGGGGCCGgag                                                                                                                                                                                                                   >  2:147873/1‑203 (MQ=255)
                                             ttgcttttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCtttt                                                                                                                                                                                     >  2:639650/1‑234 (MQ=40)
                                                         ttgcttCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                    <  2:103089/223‑1 (MQ=16)
                                                           gcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCg                                                                                                                                                                                                                      >  2:94547/1‑187 (MQ=35)
                                                                tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGTGCTGAGCCCCAGGATGGGGCCTGCAGCGCGTGCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                    >  1:1421918/1‑215 (MQ=11)
                                                                tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGAcc                                                                                                                                                                 >  1:1034875/1‑235 (MQ=255)
                                                                tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGTGCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACccg                                                                                                                                                               >  1:1078278/1‑236 (MQ=255)
                                                                tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGCCCCAGGATGGGGCCTGCAGCGCGTGCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                    <  2:1421915/215‑1 (MQ=11)
                                                                 aaGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTg                                                                                                                                                                                   <  2:1272253/216‑1 (MQ=255)
                                                                        cgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAg                                                                                                                                                        >  1:1034099/1‑236 (MQ=255)
                                                                                 aGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTa                                                                                                                                                          >  2:499729/1‑225 (MQ=255)
                                                                                   gCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCa                                                                                                                                                                                <  2:1160650/201‑1 (MQ=255)
                                                                                      aCGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGAc                                                                                                                                          >  1:222888/1‑236 (MQ=255)
                                                                                        gCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGAcc                                                                                                                                                                           >  1:474475/1‑200 (MQ=255)
                                                                                          aGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGCCGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                                               >  2:335237/1‑195 (MQ=255)
                                                                                           gATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATgc                                                                                                                                     >  1:172438/1‑236 (MQ=255)
                                                                                                      cAGGGGGCCGGGCGCCGGGCGTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCACGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGc                                                                                                                                                                                  <  2:357784/180‑1 (MQ=255)
                                                                                                                  cgcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATAc                                                                                                                <  2:890825/234‑1 (MQ=255)
                                                                                                                    cAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTAt                                                                                                                                       <  2:1417677/209‑1 (MQ=255)
                                                                                                                      gggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                                             <  2:1323854/201‑1 (MQ=255)
                                                                                                                           ggCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGGTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAg                                                                                                                                                                                                                                 <  2:1019609/111‑1 (MQ=17)
                                                                                                                               gagggaggTGACGGGCATGGGTGCCCGTCAAGCGCGGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                                                        <  1:1178678/213‑1 (MQ=255)
                                                                                                                                        gACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCCAGGGGACCCCGGAGCGTAGCGCCGGGGCCGGATGATGGGAGCCGACGGTCTTGGGTTcc                                                                                                                                                                                         >  2:173079/1‑139 (MQ=255)
                                                                                                                                          cGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCcacct                                                                                         >  2:1145669/1‑231 (MQ=255)
                                                                                                                                           gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAgg                                                                                                                                              <  2:523208/179‑1 (MQ=255)
                                                                                                                                             gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACaaaa                                                                                                                                                                        >  1:1392700/1‑151 (MQ=255)
                                                                                                                                               aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaaga                                                                                   >  2:588114/1‑234 (MQ=255)
                                                                                                                                               aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaag                                                                                    >  1:773212/1‑233 (MQ=255)
                                                                                                                                                       ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaag                                                                          <  1:18967/235‑1 (MQ=255)
                                                                                                                                                          tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGAc                                                                                                                                                                            >  2:1511020/1‑134 (MQ=255)
                                                                                                                                                             aGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                                                        <  1:1414231/183‑1 (MQ=255)
                                                                                                                                                                        cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                                                        >  1:432818/1‑173 (MQ=255)
                                                                                                                                                                                     cTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCg                                             >  1:1308213/1‑235 (MQ=255)
                                                                                                                                                                                       gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATAcc                                                                                                               >  2:535617/1‑167 (MQ=255)
                                                                                                                                                                                         aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGccc                                                                                             <  1:1270040/183‑1 (MQ=255)
                                                                                                                                                                                         aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAAt                                                >  2:661636/1‑228 (MQ=255)
                                                                                                                                                                                          gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGCGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGGGACTATGCGCCACCTGCGTGAATGAATa                                                                                                                 >  1:479569/1‑162 (MQ=255)
                                                                                                                                                                                          gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGc                                                                                                                                                     <  1:1497344/126‑1 (MQ=255)
                                                                                                                                                                                          gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGCATGAATACCTACTCGGTGTGAAt                                                                                                 >  2:156584/1‑178 (MQ=255)
                                                                                                                                                                                          gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaa                                                                              >  1:1102117/1‑197 (MQ=255)
                                                                                                                                                                                            gcgGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGCCCCGCCGTAAGGGCGGAAAGGTGACTATgcgc                                                                                                                                   <  1:1038670/142‑1 (MQ=255)
                                                                                                                                                                                               gTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGa                                   <  1:1402691/235‑1 (MQ=255)
                                                                                                                                                                                                tcctccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATc                                                                                        >  2:1330259/1‑181 (MQ=255)
                                                                                                                                                                                                tcctccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaa                                                                    >  2:664617/1‑201 (MQ=255)
                                                                                                                                                                                                tcctccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaaga                                                                            >  2:1109389/1‑193 (MQ=255)
                                                                                                                                                                                                     ccGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAAt                                                                                                 <  2:1183841/167‑1 (MQ=255)
                                                                                                                                                                                                     ccGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGc                             <  1:588116/235‑1 (MQ=255)
                                                                                                                                                                                                       gggAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGgtg                                                                                                      <  1:1425380/160‑1 (MQ=255)
                                                                                                                                                                                                        ggAGCAGGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGAc                                                                                                                                                                            <  2:660142/89‑1 (MQ=255)
                                                                                                                                                                                                        ggAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAAc                         >  1:1416082/1‑236 (MQ=255)
                                                                                                                                                                                                           gcgaGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAAc                         <  2:943484/230‑1 (MQ=255)
                                                                                                                                                                                                             aaGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaagca                                                          >  1:31581/1‑198 (MQ=255)
                                                                                                                                                                                                             aaGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAgaagcaagaagaaagacgcaagaagcaagaagcaac                                                 >  2:430092/1‑206 (MQ=255)
                                                                                                                                                                                                                gccGGGGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGAGGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGccc                                                                                             <  2:52534/159‑1 (MQ=255)
                                                                                                                                                                                                                  cGGGGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGGGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGCCccgc                                                                                                                                                              <  2:74298/93‑1 (MQ=255)
                                                                                                                                                                                                                  cGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                                                        >  2:1221202/1‑131 (MQ=255)
                                                                                                                                                                                                                       cGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACg             >  2:516299/1‑233 (MQ=255)
                                                                                                                                                                                                                         aGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGAt           <  1:1004865/233‑1 (MQ=255)
                                                                                                                                                                                                                            ggACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAAtt                                               <  2:1033718/194‑1 (MQ=255)
                                                                                                                                                                                                                             gACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGACTGTTTGGa                                   >  2:212650/1‑205 (MQ=255)
                                                                                                                                                                                                                               ccccGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGacac   >  1:1293826/1‑235 (MQ=255)
                                                                                                                                                                                                                                  cGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaagcaa                                                         >  2:691111/1‑178 (MQ=255)
                                                                                                                                                                                                                                  cGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGc        >  2:1336218/1‑227 (MQ=255)
                                                                                                                                                                                                                                  cGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGAcaca  >  2:521305/1‑233 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                 
TTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGACACA  >  AE015451/2256053‑2256510

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: