Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,256,260 C→G 9.9% intergenic (+320/‑384) 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,2600CG 9.9% 226.1 / 17.8 71intergenic (+320/‑384)asd/PP_1990Aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base C (39/25);  new base G (3/4);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-02

GCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAAC  >  AE015451/2256032‑2256487
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
gCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGc                                                                                                                                                                                                                              >  1:943486/1‑236 (MQ=255)
       aaTTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAgggg                                                                                                                                                                                                                        <  1:768242/235‑1 (MQ=255)
                 gCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGa                                                                                                                                                                                                                >  2:1144796/1‑233 (MQ=255)
                 gCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGct                                                                                                                                                                                                             >  2:756239/1‑235 (MQ=255)
                     tttAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                         <  1:239959/236‑1 (MQ=255)
                     tttAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                         >  1:528958/1‑236 (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)
                                   gTCTGCTCTGACAGAGCCGTTGCCGTGGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGAcc                                                                                                                                                                                                                     <  1:187630/210‑1 (MQ=255)
                                     cTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                         <  1:27013/220‑1 (MQ=255)
                                              cAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAg                                                                                                                                                                                                               >  2:542117/1‑205 (MQ=255)
                                                        gccgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                                                                                                                >  2:1203290/1‑178 (MQ=255)
                                                          cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGcgcag                                                                                                                                                                               <  2:1495163/225‑5 (MQ=255)
                                                          cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGGTTTTGGtt                                                                                                                                                                    >  1:429861/1‑236 (MQ=255)
                                                                cgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCCCCCGGGAGCAAGc                                                                                                                                                                                                                                     >  2:1198477/1‑165 (MQ=25)
                                                                  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)
                                                                                        gTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAgg                                                                                                                                                                                                                          <  1:747266/152‑1 (MQ=33)
                                                                                             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)
                                                                                                                              gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGc                                                                                                                                                                                                              >  2:1497341/1‑126 (MQ=33)
                                                                                                                                       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)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
GCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAAC  >  AE015451/2256032‑2256487

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: