Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 2,935,792 G→A intergenic (+87/‑175) PP_2568 → / → PP_2569 Ring‑cleaving dioxygenase/Metabolite MFS transporter, MHS family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,935,7920GA98.2% 180.2 / NA 56intergenic (+87/‑175)PP_2568/PP_2569Ring‑cleaving dioxygenase/Metabolite MFS transporter, MHS family
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/25);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGAAGGGCGCAACCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGTGGTACGACTTCTACC  >  AE015451/2935501‑2936087
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                       
ggAAGGGCGCAACCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCggg                                                                                                                                                                                                                                                                                                >  1:116722/1‑301 (MQ=255)
   aGGGGGCAACCGTCACCCTTTCTGCCGGGCAATCGAGCCGGTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTgg                                                                                                                                                                                                                                                                                             <  1:340067/301‑1 (MQ=255)
         cAACCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGTGGAGGGGCTGTCGCTGTACTGCTACGACCCGGATGGTAATCAGTTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCCGCGTGTCGTCCAGGCCGCCCTCTCGGTCAGGACCCGGCGCGTGGatctgc                                                                                                                                                                                                                                                                                       >  1:169586/1‑297 (MQ=255)
          aaCCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCACAATCGTCAGCCAACCCGTTATCGCCCGCCACCCCGGTGCTGCGTGTCGTCCCGGCCGCCCTCTCGTGCAGGACCCGCAGGGTGGATATCCa                                                                                                                                                                                                                                                                                      >  2:505678/1‑301 (MQ=255)
               gCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCaa                                                                                                                                                                                                                                                                                 <  1:687563/301‑1 (MQ=255)
               gCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCaa                                                                                                                                                                                                                                                                                 <  1:540515/301‑1 (MQ=255)
                cACCATTTCGGCCTGGGCATCGAGCCGTTCGACGGGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAAt                                                                                                                                                                                                                                                                                <  2:321618/301‑1 (MQ=255)
                    ttttCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGACTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCg                                                                                                                                                                                                                                                                            <  1:56698/300‑1 (MQ=255)
                            cccgcggagcgagacgtgcgcggcggtggcgTGGACGGCCTATCTGGAGGCCGCCGGGGTAAAGTTGGAGCCTGCCGTGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTGGCTGTACTGCTACGACCCGGATGGGAATCAGGTCGAGTTGAAGGGGCCGGGAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGc                                                                                                                                                                                                                                                                     <  1:809055/273‑1 (MQ=255)
                                gcATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGCTTAttcttcct                                                                                                                                                                                                                                                                >  2:336521/1‑299 (MQ=255)
                                   aTCGAGGCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGa                                                                                                                                                                                                                                                              <  1:552036/301‑1 (MQ=255)
                                                 cgaGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGTACCCGGCGGGTGTATATCGAGCCAATccccaatctgcttccataacctttcaacgatc                                                                                                                                                                                                                                                >  2:171386/1‑272 (MQ=255)
                                                   gggCCGCGTTGGCGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGGGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCg                                                                                                                                                                                                                                              <  1:150716/300‑1 (MQ=255)
                                                                    tATCTGGAAGCCGGCGGGGTAAAGGTGGGGCATGCCGAAAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                                                                                                             <  1:686872/301‑1 (MQ=255)
                                                                           aaGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTc                                                                                                                                                                                                                      >  1:393180/1‑301 (MQ=255)
                                                                                 ccgGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCCATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCaa                                                                                                                                                                                                                >  1:162623/1‑301 (MQ=255)
                                                                                       gaaggggtgggCCTGCCGGGAAGGGTTATGGTGCGGAAGGGCAGGGGCTATCGCTGTACTGCTACGACCCGGATGGTAACAAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAACCCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCGGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaag                                                                                                                                                                                                          <  2:468262/291‑1 (MQ=255)
                                                                                                  ccTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCCCGTGCAGGACCCGGCGGGCGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCCCGCAACCAAAGCAAGCGCAAgg                                                                                                                                                                                               >  1:57269/1‑301 (MQ=255)
                                                                                                        gagaAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCACTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTGTTGCTTCGCGCAACCAAAACAAGCGCCAGGGCCTgg                                                                                                                                                                                         >  2:799030/1‑301 (MQ=255)
                                                                                                          gaAGGGTTAGGGTGCGGAAGGGGGGGGGCTTTCGCTGTACTGCTCCGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGTTAGTGCCGGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGt                                                                                                                                                                                       <  2:591044/301‑1 (MQ=255)
                                                                                                                tgaTGGTGGCGAAGGGGAGGGGCTTTGGCTGTACTACGACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACg                                                                                                                                                                                 <  1:563369/299‑1 (MQ=255)
                                                                                                                                 aggggCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCCGCTGGCGGATATCCAGTCAATCCCGATTATGCTTCATGAACATTGCTGCGATTCCAGCCGGTCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGCCGGTGAACGGGTTGTGAgg                                                                                                                                                                >  2:95737/1‑301 (MQ=255)
                                                                                                                                        ttCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTccc                                                                                                                                                         >  1:362267/1‑300 (MQ=255)
                                                                                                                                         tCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCgg                                                                                                                                                        >  1:434860/1‑301 (MQ=255)
                                                                                                                                           gCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGcc                                                                                                                                                      >  1:355283/1‑301 (MQ=255)
                                                                                                                                             tGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCTGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCTGCCACCCTCGTCCTGCTTCGCGCAACCCACGCAAGCGCCAGAGCCTGGGTGCGACTGCGACGTGGTTGTCAGGGACGCCGGCCTg                                                                                                                                                    >  2:344932/1‑301 (MQ=255)
                                                                                                                                             tGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCCCGCAACCAAAACAAGCGCCAGGCCCTGGGTGCGACGGTGAAGTGGGTGTGATGTGCTCCGGACTg                                                                                                                                                    >  1:693659/1‑301 (MQ=255)
                                                                                                                                                  tGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCACAACCAAAGCAAGCGCCAGGGCCTGGCTGCGACGGTGCAGTGGTTGTGAGGTGCTCCGCCCttaatta                                                                                                                                               >  2:633675/1‑295 (MQ=255)
                                                                                                                                                     tACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGACATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCCATCCGAGGCGGCCCGCCTCCTCGTGCTTCGTTCAACCAAAGCAAGAGCCCGGGCCTGGGTGCTACGGTGAAGTTGGTGTGACGCGCTCCGGCCTGTCTGacaa                                                                                                                                            >  2:822710/1‑301 (MQ=255)
                                                                                                                                                        gACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGCCGGAGAAGTGGGTGTGAGGTGCCCCGGCCTGTTTGGCAACAt                                                                                                                                         >  1:719957/1‑301 (MQ=255)
                                                                                                                                                                  gTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCATGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCa                                                                                                                               >  1:590992/1‑301 (MQ=255)
                                                                                                                                                                           gcggtttGAAGGTGCGGGTAGTCCCTGGCGGGCGAAACCCCCCCCGCCAAAATCGTCATCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGc                                                                                                                      <  2:393055/296‑1 (MQ=255)
                                                                                                                                                                                ttGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATTGCgggtt                                                                                                                 >  2:655815/1‑298 (MQ=255)
                                                                                                                                                                                 tGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTtct                                                                                                                >  1:331193/1‑301 (MQ=255)
                                                                                                                                                                                  gcaGGGGCCGGTGGGGCTCGGCGGGTGAAACCGTCCCCGCCAAAATCGTCAGCCAACCCGATATCGCCCGGGACCCCGGTGCTGCGTGTCGTCCAGGCCGCACTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTtctc                                                                                                               <  2:546570/299‑1 (MQ=255)
                                                                                                                                                                                  gaagtgggcgggaGTGCCTGGCGGGTGATCCAGCCCCCGCACAAGTCGGCAGCCAGCCGGCTGTCGCACGGCGCCCCGGTGCTGCGTGTCGTCCAGGCTGCCCTCTCGTGCAGGGCCCGGCGGGGGGATAGCCAGCCAAGCCAGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCTGCCCGCATCGTCGTGCTTCGCGCAACCAAGACAAGCGCCAGGGCCTGGGTGCGACAGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATAACGACATGACGGAGGCAACCCATGGCGGTtctc                                                                                                               <  2:169520/289‑1 (MQ=255)
                                                                                                                                                                                       ggCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCCTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCGGTATGACAAAAAGACGGACGCAACACAGGGAGGTTCTCTACac                                                                                                          >  1:464645/1‑300 (MQ=255)
                                                                                                                                                                                          gggTAGTGCCTGGCGGGTGAACCCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGTCACCCCGGTGCTGCGTGTCGGCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGcac                                                                                                       <  2:116683/300‑1 (MQ=255)
                                                                                                                                                                                          cGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGcac                                                                                                       >  1:89088/1‑301 (MQ=255)
                                                                                                                                                                                           ggtgttgcgtGGCGGGTGAAACGGCCCCGGCCAAAATCGTCAGCCAACACTCTATCGCCCGGCACCCCGGGGTTGCGGGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGATGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGcaca                                                                                                      <  1:4539/292‑1 (MQ=255)
                                                                                                                                                                                                 gcctggtgtggaaacccGCGCCCGCCAAACTCGTCAGCCAACCCGTTATCGGCGGGCACCCCGGTGCTGCGTGTCGTCCGGGCCGCCCTCTCGTGAAGGACCCGGGGGGTGGATAGCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACg                                                                                                <  1:302191/286‑1 (MQ=255)
                                                                                                                                                                                                     ggCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCa                                                                                            <  2:89069/301‑1 (MQ=255)
                                                                                                                                                                                                                              aTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTTCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCACCGCGGTATCAccca                                                                   >  1:548681/1‑299 (MQ=255)
                                                                                                                                                                                                                                    ggCCAACCCGCTCTCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCGCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGGACCTTGCTGTGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAggagga                                                             <  2:57257/300‑1 (MQ=255)
                                                                                                                                                                                                                                       cAACCCGCCATCGCCCGGCGCCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACACGGCGGGTGGATATACAGCCAATCCCGATTATGCTGCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAgcg                                                          <  1:670423/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gcGTGTCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCAc                            <  1:505845/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gtgtCGTCCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACgg                          >  2:203271/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                          tgttcCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATACAGCCAATCCCGATTATGCTTCATGAACCTGGCTGCGATCCGAGACGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTAt                       <  1:335569/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gtcCAGGCCGCCCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCACCGCGGTATCACCAAGGAGGAGCGCACGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTc                     >  1:517671/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                            gtcCAGCCCGCCCTGTCGTGCAGGACCCGGCGGGTGGATATCCACCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTc                     <  2:331082/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     gccCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTTGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGCCAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGCAGCGCAAGGTCATCTTCGCCTCTCCCCTGGTCAAGGTATTctgttggtac            >  2:503801/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                                     gccCTCTCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCACAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGTGGTAc            <  1:633898/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ctctCGTGCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGTGGTACGAc         >  2:582236/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           ccGTGCAGGACCCGGCGGGTGGAAATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGTGGTACGACTTc      <  1:355675/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               gCAGGCCCCGGCGGGTGGATATCCAGCCAAGCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAACAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGTGGTACGACTTCTAcc  <  1:799307/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               gCAGGACCCGGCGGGTGGATATCCAGCCAATCCCGAATATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCTCGCCTCGTCGTGCTTCGCGCAACCAAAGCAACCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCTTGGGCACGGTATTCGACGGGTACGACTTCTAcc  >  1:207255/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                       
GGAAGGGCGCAACCTGCACCATTTCTGCCTGCGCATCGAGCCGTTCGACGAGGCCGCGTTGACGGCCTATCTGGAAGCCGCCGGGGTAAAGGTGGAGCCTGCCGAGAAGCGTTATGGTGCCGAAGGGGAGGGGCTTTCGCTGTACTGCTACGACCCGGATGGTAATCAGGTCGAGTTGAAGGGGCCGGTAGTGCCTGGCGGGTGAAACCGCCCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCCCCAACGCGGTATCACCAAGGAGGAGCGCAAGGTCATCTTCGCCTCGTCCCTGGGCACGGTATTCGAGTGGTACGACTTCTACC  >  AE015451/2935501‑2936087

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