Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 353,067 A→G intergenic (+34/+84) hisF → / ← cbcV Imidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451353,0670AG98.0% 166.9 / NA 50intergenic (+34/+84)hisF/cbcVImidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/26);  minor base C (1/0);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATCCTGCTGACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTGCCCATTTCTGCATGTCCAG  >  AE015451/352779‑353360
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                       
aagaTCCTGCTGACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGGGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                                                                                                                                                                                             <  3:523741/300‑1 (MQ=255)
       tgctgACCAGCATGGACCAGGACGGCATGAAGAACAGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGATTTGCAGCCCTTCGCGGGCACGCCTGCCCCcaccg                                                                                                                                                                                                                                                                                      >  1:147451/1‑299 (MQ=255)
         ctgACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGGTGAGCCATTGAACGCGGCTGCTTTGCAGCCCTTCGCGGGCACGGCCGCGCCCACAggg                                                                                                                                                                                                                                                                                    >  1:570554/1‑300 (MQ=255)
          tgACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGACTACATGGCTTCGCGAGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGAGGGCACGCACGCTCACACAGGtg                                                                                                                                                                                                                                                                                   >  1:338088/1‑300 (MQ=255)
            cccAGCATGGACCAGGACGGCATGAAGATCGGCTTCGACCTGTGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCTCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                                                                                                                                                                                                 <  3:338089/300‑1 (MQ=255)
                              ggCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAg                                                                                                                                                                                                                                                               <  3:257836/301‑1 (MQ=255)
                               cagacgtgtgctcttccgatcTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                                                                                                                                              <  3:98802/281‑1 (MQ=255)
                                  tGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTTCAGCCCTTCGCCGGCACCCCCCCCCCCACCTGCTTGCTgagcactcctgtggagcggg                                                                                                                                                                                                                                                           >  2:309286/1‑284 (MQ=255)
                                  ggaagaaCGGCTTGGACCTGGGCGTGACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCggg                                                                                                                                                                                                                                                           <  1:161204/300‑1 (MQ=255)
                                          ggcttagcccTGGGCGTTCCCCGGGCTATAAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcg                                                                                                                                                                                                                                                   <  3:570556/293‑1 (MQ=255)
                                          ggCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcg                                                                                                                                                                                                                                                   <  1:363507/301‑1 (MQ=255)
                                            cTTCGACCTGGGCTTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGTTGATTGCCTCGTGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGAATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                                                                                                                                                                                 <  3:147453/301‑1 (MQ=255)
                                                 cccTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACGTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                                                                                                                                            <  4:629781/300‑1 (MQ=255)
                                                  ccTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGCTCTCTCTGGGCGCGGGCGCGCCcgctgaggcg                                                                                                                                                                                                                                           >  2:629578/1‑294 (MQ=255)
                                                     gggCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACTGGTTTTCTGCGTTCTCtgttggtgcagctcggaagaccggtgtgttt                                                                                                                                                                                                                                        >  1:98799/1‑273 (MQ=255)
                                                        ttccgatctggctcgcggTGATGCGCTGAGCATCCGGGTGGTGGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGAGCGTGCCCGCGAAGGGCCGCGaa                                                                                                                                                                                                                                     <  3:112524/284‑1 (MQ=255)
                                                         gTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGTTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGCTCTCTGTGGGAGCCGGCGTGCCCGCGAAGGGCCGCGCAGc                                                                                                                                                                                                                                   >  4:256850/1‑301 (MQ=255)
                                                           tACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCg                                                                                                                                                                                                                                  >  1:231559/1‑301 (MQ=255)
                                                            ccgctgtccctacacgactctctgccgatcccGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGTCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACGTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                                                                                                                                                                                 <  1:110532/271‑1 (MQ=255)
                                                                 ggCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGaaagatcggag                                                                                                                                                                                                                            >  1:112521/1‑292 (MQ=255)
                                                                                      cgatctggTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                                                                                                                                                                                       <  4:518022/295‑1 (MQ=255)
                                                                                            ggTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCCCCGTGGGAGCGGGCGTGCCCGTGAAGGGCCGCGCAGCGGCCCCAATTGCTTCCACCgctgttcctattgtc                                                                                                                                                                                                 >  2:517872/1‑287 (MQ=255)
                                                                                              ggATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                                                                                                                                                                               <  3:231557/300‑1 (MQ=255)
                                                                                                    ccTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcgg                                                                                                                                                                                         >  1:471878/1‑300 (MQ=255)
                                                                                                                              atcttGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCGAGTATcc                                                                                                                                                              <  4:106148/297‑1 (MQ=255)
                                                                                                                                  tGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGTCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCCGCGCGTGGTAAAGCCCGGGATcccggcgatccagca                                                                                                                                                           >  2:106124/1‑289 (MQ=255)
                                                                                                                                   ggCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTGGGTATCCCCGAGTAtatccac                                                                                                                                                          >  4:103431/1‑295 (MQ=255)
                                                                                                                                          ggCATTATGGAAGGGCAGGCCAGCGCGGTGCTGGCTGCCAGCAGCTTCCACTTTGGCGAGTACACAGGGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTg                                                                                                                                                  <  2:125481/301‑1 (MQ=255)
                                                                                                                                             cgcaccctctccgatccaaGCGCGGTGCGGGCGGCCACCATCTTCACCTTGGCGTAGTACACCGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                                                                                                                               <  2:653130/283‑1 (MQ=255)
                                                                                                                                                         cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGCTCTTGCCGAGCAGCTCGTGGTACAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                                                                                                                                   >  4:246376/1‑301 (MQ=255)
                                                                                                                                                              cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCgagatcggaagagcacac                                                                                                                              >  4:653341/1‑284 (MQ=255)
                                                                                                                                                                ggggggggcgggcggcCAGCAAGTTCAACTGTGGCGGGTACACAGTGCCGGAGGCCAAGGCCAACATGGCTTCGCGCGGGATCGTCGTTGGCTGATCCATTGACCGGGGCTGTTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCt                                                                                                                           <  2:256792/291‑1 (MQ=255)
                                                                                                                                                                 ggcgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt                                                                                                                          <  4:309359/300‑1 (MQ=255)
                                                                                                                                                                        cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCAGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCACCGCGGGGTAAAGCCCGGTATCCCCGAGTATCCCCACCACCCGGTCGTTagatcggaagagaacctgtctgaac                                                                                                                    >  4:83840/1‑277 (MQ=255)
                                                                                                                                                                             ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTg                                                                                                              <  3:471880/301‑1 (MQ=255)
                                                                                                                                                                                 aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCAAATTGCTTAAACCGTGGTTTTTGCCGAGCAGCGCGTGGTAAAGCCCGGTTTCCCAGCGTACCCCACCCACCCTGTCGTTATCTTGCACACCCAGCTTATTGCCGGTTTGATaa                                                                                                          >  1:80754/1‑301 (MQ=255)
                                                                                                                                                                                        tCCACTTTGGCGAGTACACATGGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTTCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCt                                                                                                   <  3:80756/301‑1 (MQ=255)
                                                                                                                                                                                         tttcccctccacgccgcctttcctatctggCCAAGGCCTACAGGGCTTCGCGCGGGTTCGTCGTTCGCTGATCCATTGCCCGGGGCGGCTTTGCAGCCTTCCGCGGGCACGCCCGTTCCAACAGGTTTGCTGGGATCTCTGTGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTg                                                                                                  <  1:447841/273‑1 (MQ=255)
                                                                                                                                                                                         cgtaatgtgtagagaaccgagggCCGGAGGCAAAGGCCTACATGGCGTCGCGATGGATCGTGGTTCGCTGATCCATTGACCGGAGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCt                                                                                                   <  2:103407/279‑1 (MQ=255)
                                                                                                                                                                                                                    aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCgg                                                                        >  4:531854/1‑301 (MQ=255)
                                                                                                                                                                                                                     ggCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCCTGTTGCCGGTTTGATAACGGATCTgagatcggaagagcaaacgtctgaactcc                                                                      >  3:447843/1‑273 (MQ=255)
                                                                                                                                                                                                                     atctaaGGCCTACCTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGc                                                                       <  4:161920/297‑1 (MQ=255)
                                                                                                                                                                                                                         aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTCTCCTGCAGCACCAGCTTGTTGCCGGTTTGATCACTAATCTGCAGTCCCTCCCGCATGCCGATGTTGGccgat                                                                   >  2:161882/1‑297 (MQ=255)
                                                                                                                                                                                                                          aGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGca                                                                  >  4:736154/1‑301 (MQ=255)
                                                                                                                                                                                                                                         gctcttccaatctcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTgc                                                   <  4:589189/289‑1 (MQ=255)
                                                                                                                                                                                                                                              gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCCGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACCCATGCCGATGTCGGGGTGCACCCCGGGCGGCCtcggcccc                                              >  2:335905/1‑294 (MQ=255)
                                                                                                                                                                                                                                               ggaccgtggtgcgcggttCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCca                                             <  2:735894/284‑1 (MQ=255)
                                                                                                                                                                                                                                                      cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATACCGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCTTGccccccgtctcgccctccactcggtcgcac                                      >  2:589004/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                 ccagagcacgggcatgattcgCAGCCCGTCGCGGGCACGCCCGTGCCCACAGGTTGGCTCAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGcc                           <  2:531704/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTGCCCATTTCTGCATGTCCac  >  1:13365/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                       
GAGATCCTGCTGACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTTGCCCCGGTGCCCATTTCTGCATGTCCAG  >  AE015451/352779‑353360

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