Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,348,9550GA32.7% 68.4 / 35.4 52intergenic (‑320/+154)PP_3820/galUGroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/23);  new base A (9/8);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCACCGCAGTCGAAACGCTTGCCCTTGAACTTATAGGCAATCACGCAGCCATCCTTCGCCTG  >  AE015451/4348656‑4349240
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                              
gtgcggtgcaCCGCCACTTAGTACTGCTTCAACTGACGTCTGCACGCTCATCCCCTTGCCTTTCGACGCGCGGTAGCACCGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGACTCCGGCTAGc                                                                                                                                                                                                                                                                                               <  2:1257/292‑1 (MQ=11)
       tAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTccc                                                                                                                                                                                                                                                                                        >  1:438445/1‑301 (MQ=14)
          ctgCCACTTACGACGGCTTCGGCTGACTTCTGCACGCACATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGTGCTCCCTGGCTAATTCGCTCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCGAGTATGGGGCTTTGGCGATTTGTGCCACCGCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCtt                                                                                                                                                                                                                                                                                     <  2:271453/299‑1 (MQ=11)
          ccGCCACTTACTACGTCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCtt                                                                                                                                                                                                                                                                                     <  1:142556/301‑1 (MQ=14)
            gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGCCTCACACTTACCCtttc                                                                                                                                                                                                                                                                                   >  2:506222/1‑299 (MQ=9)
              cACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCg                                                                                                                                                                                                                                                                                 >  1:718836/1‑301 (MQ=14)
               aCTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCgg                                                                                                                                                                                                                                                                                >  1:641025/1‑301 (MQ=33)
                cTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCggg                                                                                                                                                                                                                                                                               <  1:682963/301‑1 (MQ=14)
                 ttACTACGGCGTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGgt                                                                                                                                                                                                                                                                              <  2:548237/301‑1 (MQ=14)
                  tactacGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGTCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCCGgtg                                                                                                                                                                                                                                                                             >  1:579932/1‑301 (MQ=18)
                    ctacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCTCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCCCAGCGTTCCGTGCAGGTATTGTGCTTTGGCGATTTGTTCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCCGGCCAGCATTTTGCCTCGTGCTTCCTTCAGATTCGCGGTCACCCGCGCCACCCTGCCCTCCTGCTAACACTTCCCCTTCCCGGgtgtt                                                                                                                                                                                                                                                                           >  1:34617/1‑300 (MQ=11)
                               gCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGCGGACttt                                                                                                                                                                                                                                                                >  2:68117/1‑301 (MQ=14)
                                cTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTc                                                                                                                                                                                                                                                               >  1:641263/1‑301 (MQ=14)
                                 tGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTcc                                                                                                                                                                                                                                                              >  1:547143/1‑300 (MQ=14)
                                  gACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCAc                                                                                                                                                                                                                                                             >  1:292794/1‑301 (MQ=14)
                                         gCACGCTCGTCCAGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACCGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAg                                                                                                                                                                                                                                                      <  2:710237/301‑1 (MQ=9)
                                           aCGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTc                                                                                                                                                                                                                                                    >  1:665188/1‑301 (MQ=33)
                                                 gtCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCAc                                                                                                                                                                                                                                              >  1:45873/2‑301 (MQ=14)
                                                    ccGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTg                                                                                                                                                                                                                                           <  1:131025/301‑1 (MQ=14)
                                                     cGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCCCTTgt                                                                                                                                                                                                                                          >  1:303110/1‑300 (MQ=18)
                                                             ggcgaggcgcgGTAGCGCAGGGGCAAGAGGGCAGATCTCCCAAGGGAATTCCCGCGACCTTCCTGCTTATGCCTGTCGAATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCAGCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACcaca                                                                                                                                                                                                                                  <  2:207132/295‑1 (MQ=11)
                                                               ggACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAAc                                                                                                                                                                                                                                <  2:640991/300‑1 (MQ=14)
                                                               cGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAAc                                                                                                                                                                                                                                <  1:433528/301‑1 (MQ=14)
                                                                gACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCATGGTAAGTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAAcc                                                                                                                                                                                                                               <  1:398663/301‑1 (MQ=18)
                                                                 aCGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACATCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCa                                                                                                                                                                                                                              <  1:657541/301‑1 (MQ=18)
                                                                   gcgcgGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATCCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTGGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAg                                                                                                                                                                                                                            <  1:301168/301‑1 (MQ=18)
                                                                            cacaGTTGCAAGCGTCCATATCTGCCAGTGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATGTGCGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCCCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTgcg                                                                                                                                                                                                                   <  2:159758/301‑1 (MQ=11)
                                                                               aGTGGCAAGCGTACAGAACTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCtt                                                                                                                                                                                                                <  1:511119/301‑1 (MQ=18)
                                                                                 tGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCTACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTAc                                                                                                                                                                                                              <  1:322158/301‑1 (MQ=17)
                                                                                        cgtgcagctctcCCAGGGTAATTCGCGCGATCTTCCTGCTTAGGCCTGTCGGTTCTACGTCACAGCGTTCCGTGCAAGTATTGTGCTTTGGCGATTTGTGCCACCTCACCCAGCTGTGCCGCCTCTGTCCGCGTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCACGGTCACACGCGACACCCTTACCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAc                                                                                                                                                                                                       <  2:716445/293‑1 (MQ=11)
                                                                                        cGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAc                                                                                                                                                                                                       <  2:640753/301‑1 (MQ=31)
                                                                                         gTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCCCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACTCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAATCGGATTGCGCttccccgaaacc                                                                                                                                                                                                      >  2:112845/1‑291 (MQ=9)
                                                                                                    ccAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCCGGTTGCGCTTAAGCGCAACCCGCCATGCCtc                                                                                                                                                                                           >  1:3959/1‑300 (MQ=14)
                                                                                                                gcgcgACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACaaaa                                                                                                                                                                              <  1:137144/301‑1 (MQ=31)
                                                                                                                          ccTGCTTATGCCTGTCGGATCTACGCAACAGCCTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCGCACCCCGCTGTGCCGCCTCTACCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCAGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCAcccc                                                                                                                                                                    <  2:168185/301‑1 (MQ=255)
                                                                                                                                      tGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACCTTCACTTCCAAGTCCCCCACTTGGCCACCACAACCAAGCTGGTTGCGGTTCGGCGCAACGCGCCAGTCCTGGCGACAACAAAAAAGCCCCCCCGAAGGATGGCtt                                                                                                                                                        >  1:319163/1‑301 (MQ=11)
                                                                                                                                           gATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATGTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAgg                                                                                                                                                   <  2:702988/301‑1 (MQ=39)
                                                                                                                                               tACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCcag                                                                                                                                               >  2:502255/1‑301 (MQ=40)
                                                                                                                                                           ttCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCa                                                                                                                                   <  1:502463/301‑1 (MQ=255)
                                                                                                                                                            tgcGTGGAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGGTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAg                                                                                                                                  <  2:573001/299‑1 (MQ=255)
                                                                                                                                                                            gCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACCACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCAACCAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGATTATCAGTAAGCCTTTCCAGTCt                                                                                                                  >  2:155509/1‑301 (MQ=255)
                                                                                                                                                                               ttGGCGCTTCGTGCCACCTCACCCCGCTGTGCCTCCTCTATCCCCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGTGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgt                                                                                                               <  1:40191/301‑1 (MQ=255)
                                                                                                                                                                                 ggCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtag                                                                                                             <  1:259027/301‑1 (MQ=255)
                                                                                                                                                                                                              ccgccTCTATCCGCGTCCTGTGCGTCAGGCCAGCATTTTGCCTCTGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTc                                                                                <  1:96223/301‑1 (MQ=255)
                                                                                                                                                                                                                    ctttCCGCTTCCTGTTCGGCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACACGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTa                                                                          <  1:205106/298‑1 (MQ=255)
                                                                                                                                                                                                                     ggtCCGCTTACTGTGCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAg                                                                         <  2:664905/299‑1 (MQ=255)
                                                                                                                                                                                                                                 tttCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCa                                                             <  1:262370/300‑1 (MQ=255)
                                                                                                                                                                                                                                 tttCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTGGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCa                                                             <  2:579682/300‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTGATCAGTAAGCCTTGCCAGTCTTGGAGGTGGTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCATTCGGCGccc                                                             >  2:465982/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          gggCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCACCGCAGTCGAAAAGCTTGCCCTTGa                                    >  1:295916/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           cccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCACCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCACCGCAGTCGAAACGCTTGCCCTTGAACTTATAGGCAATCACGCAGCCATCCTTCGCct   <  1:112893/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            ccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCACCGCAGTCGAAACGCTTGCCCTTGAACTTATAGGCAATCACGCAGCCATCCTTCGCctg  >  1:210907/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                              
GCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCACCGCAGTCGAAACGCTTGCCCTTGAACTTATAGGCAATCACGCAGCCATCCTTCGCCTG  >  AE015451/4348656‑4349240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: