Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 4,305,124 G→A G42G (GGG→GGA LY180_20645 → rhsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065844,305,1240GA97.1% 177.6 / NA 70G42G (GGG→GGALY180_20645rhsA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (31/36);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACA  >  CP006584/4304829‑4305312
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                            
gcagcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAGAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGt                                                                                                                                                                                            <  1:849352/298‑1 (MQ=255)
   gcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCGGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGAc                                                                                                                                                                                         <  1:691444/298‑1 (MQ=255)
      attattccggtatgtCCGAAAAAAGCGCGTGACGGAGCGTTTTTTGTGCCGGTAATTTGTCTTTTATAAATTCTGAAAGATTCCGGGGTGTGAAATCATTATCATTTCCGGTTGGTGACTGGCGCTGGGGTTCTGGTTCTTATGTAGTTCAAAAGGATTTTACGATGAGCGGAAAAGCGGCGGCGCGTCAGGGCGGCATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATGGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc                                                                                                                                                                                      <  2:550117/284‑1 (MQ=21)
      aaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc                                                                                                                                                                                      <  1:179368/298‑1 (MQ=255)
      aaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGCAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc                                                                                                                                                                                      >  1:566002/1‑298 (MQ=255)
       atgagagCGGAGTTCCGACAAAAGTGCCTGCCGGTGCTTTTTTTGTGCCTGTAATTTGCCGTTTATAAATTCTGAAAAGTTCCGGGGTGTTATCTCAGTATCATTCACTGTTGTTGACTGGCGCTAGGCTTCTTGTTCTTATTTATTTCAAAGTGATTTTAATATGAGCGGAAAACCGTCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGGCGTcc                                                                                                                                                                                     <  2:562925/295‑1 (MQ=21)
          aGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGCCGCATTATGGCGGAAGCATTGTCCAGGGTGCAGCCGGCGTCCGTATTGGTGCCCCCACCGGGGTGGCCCGTTCGTCGTGCACCCGCGGAGTGCCGCCCgcc                                                                                                                                                                                  >  2:442056/1‑296 (MQ=255)
                 aTTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCgg                                                                                                                                                                           >  1:346458/1‑298 (MQ=255)
                 aTTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCgg                                                                                                                                                                           <  1:897992/298‑1 (MQ=255)
                        aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGGTGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATcc                                                                                                                                                                    <  2:100276/298‑1 (MQ=255)
                         aaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATccc                                                                                                                                                                   >  1:221078/1‑297 (MQ=255)
                           caaGCGCCGCAGGGTTCTTTTTTTGTGCCCTTATTTTGTCTGTTATAACTTCTGAAAAATTCCGTGGTGGTTACTCATTATCATTTCCTGGTGTTTGCTGGCGTTAGGTTTCTGGTTGTTATTTAGTTCAAAAGGATTTTACGAGGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTGCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCTGAGTGACGTCCGGCCATCCGGTCAATCCcct                                                                                                                                                                 <  2:690729/297‑1 (MQ=18)
                            aaGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCcctg                                                                                                                                                                >  1:121581/1‑298 (MQ=255)
                             aGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCcctgc                                                                                                                                                               <  1:695141/298‑1 (MQ=255)
                              gcgcCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCcctgct                                                                                                                                                              >  1:375760/1‑298 (MQ=255)
                                gccgcACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCg                                                                                                                                                            >  1:583779/1‑298 (MQ=255)
                                  cgcACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGt                                                                                                                                                          >  1:43026/1‑298 (MQ=255)
                                              ttttgGTGCCTTTAATTTGCTGTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGGCTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCtt                                                                                                                                              <  2:346458/298‑1 (MQ=255)
                                              ttgtgtgtccTTTAATTTGCTTTGTATAACTTCTGAAAAGTTCCGGGGTGTTACCTCCTTTTCATTTCTGGTTGTTGACTGGAGGTAGGCTTCTGGTTCTTAGTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCtt                                                                                                                                              <  2:376898/290‑1 (MQ=18)
                                                   gtgCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCgg                                                                                                                                         <  1:347839/298‑1 (MQ=255)
                                                    ggccgggaaTTTGCCTGTGATAAATTCTGAAAAAGTCCGGGGTGTTAACTCATTATCTTTTGCTGTTGTTGACTGGCGCTAGGCTTTTGGTTCTTATTGATTTCAAAAGGAGTTTACGAGGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGt                                                                                                                                        <  2:600212/291‑1 (MQ=18)
                                                      ccTGTAATTTGCCTTTTTTAAATTCTGAAAAATTCCTGGGGTTTAACTCATTTTCATTTCCTGTTGTTGGCTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGa                                                                                                                                      <  2:487416/298‑1 (MQ=255)
                                                      ccTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATCACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTCCGTTTTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGGGGAGGGCCGTCCGGCCACCCGGTCAATCCCCCCCCCGGGCCAAACGCCCTTCCCCGtca                                                                                                                                      >  2:981540/1‑296 (MQ=25)
                                                         gTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAc                                                                                                                                   >  1:433644/1‑298 (MQ=255)
                                                         gTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTCCGGTGTGCCCCGGCGGAGTGACGTCCGCCCATCCGGTCAACCCCCCGCTCGGTGCAAAGGCCCTTCCCGGTGacac                                                                                                                                   >  2:685614/1‑298 (MQ=255)
                                                          tAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTCGTTCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGCGTGACCTCCGGCCATCCCGTCAATCCCCTCCCCGGTGCAAAGGTCCTTCCCGGTTGAAcc                                                                                                                                  >  2:608960/1‑298 (MQ=25)
                                                          tAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTGGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAcc                                                                                                                                  <  2:813288/298‑1 (MQ=255)
                                                            aTTTTCTTTTTATAAATTGTGAAAAATTCCGGGGTGTTAACTCATTTTCATTTCCGGTTGTTGCGTGGCGCTAGGCTTCTGGTTCTTATTTATTTCGAAAGGGTTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGAGATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCGCGGTGCAAAGGTCCTTCCCGGTGAAACCGa                                                                                                                                <  2:945371/298‑1 (MQ=18)
                                                              ttGCCTTTTCTAAATTCTTAAAAATTCCGGGGTGTTAACTCTTTTTCATTTCTGGTTGTTGACTGGCGCTAGGGTTCTGGTGCTTATTTATTTACAAAGGATTTTACGATGAGCGGAAAGCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGTGCAGCCGGGGTACGTATTGGTGGCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACa                                                                                                                              <  2:209689/298‑1 (MQ=18)
                                                                   tttttgacatttgaaaaaaGTCCGGGTGGTTAATTCATTGTCTGTTCCTGGTGGTTACTGGCGCGAGGCTTCTGGGTCTTATTGATTTCGAAAGGATTTTACGGTGAGCTGGAAGCCGGCGGCGCGTCAGGGCGACATGACGAAGGATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGcc                                                                                                                         <  2:69046/284‑1 (MQ=17)
                                                                          aaTTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccg                                                                                                                  >  1:871804/1‑298 (MQ=255)
                                                                                  aaaaTTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTg                                                                                                          >  1:241932/1‑298 (MQ=255)
                                                                                    aaTTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGcc                                                                                                        >  1:799942/1‑298 (MQ=255)
                                                                                     agtCCGGGGTGTTACCTTATTGTATTTTCCTTTGGTTGACTGGCGGTAGGCTTCGTGTTCGTATTTATTTAAAAAGGGGTTTACGATGAGCGGAGACGCTGCGGCGGGTCAGGGCGACATGACGCAGGATGGCGGTAGCATTGTCCAGGGTGCGGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                                                                                                       <  2:859396/296‑1 (MQ=14)
                                                                                       tCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGtt                                                                                                     <  1:270174/298‑1 (MQ=255)
                                                                                         cGGGGTTTTACTTTATTTTCTTTTCGTGTTGTTGCATGGGGCTAGGCTTTGGTTTCTTATTTATTTCAAAAGGAGTGTAGGGTGAGCGGAAAACCGGGGGCGCGGCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGGCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCa                                                                                                   <  2:578892/298‑1 (MQ=14)
                                                                                           gggTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATc                                                                                                 <  1:942914/298‑1 (MQ=255)
                                                                                             gtgtTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCct                                                                                               >  1:343113/1‑298 (MQ=255)
                                                                                                ttAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGAAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGTGTACGTATTGGTGCCCCCACCGGCGCGGCCTGTCCGCTGTGCCCCGGCGGGGTGACGTCCGGCCATCCCGACCATCCCCTGTCCGACGCCACGGTCCTCCCCGCTGAGACCCCCACCCCCCTCCCCGGCCCCCACCCGTCCAGCctccc                                                                                            >  2:532270/1‑296 (MQ=18)
                                                                                                          aTCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTAc                                                                                  <  1:981540/298‑1 (MQ=255)
                                                                                                                             tGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                               >  1:869378/1‑298 (MQ=255)
                                                                                                                                 gcTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGcc                                                           >  1:897219/1‑298 (MQ=255)
                                                                                                                                        ttCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGTGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCg                                                    <  1:194520/298‑1 (MQ=255)
                                                                                                                                        ttCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCCGAGTGACGCCCGGCCATCCGGTCAATCCCCTGCTCGGTCCAAAGGTCCTTCCCGCTGCACCCCACCTCGCCCTGCCCCGCCCGCCGCCGCTCACCCTCCCCCGCCCCCACAGCAGTTCCCCGGCAAAACCGCCCGCACCg                                                    >  2:809947/1‑298 (MQ=14)
                                                                                                                                          cTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGt                                                  <  1:608960/298‑1 (MQ=255)
                                                                                                                                           tgattattttttatttaAAATGGGTTGCACGGTGAGCGGAAGACCGGGGGCGCGGCAGGGCGACATGAGGCAGTATGGCGGTAGCAGTGTCCAGGGTTCAGCCGGGGTACGTGTTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTg                                                 <  2:118239/289‑1 (MQ=11)
                                                                                                                                                  tatttattTCAACAGGATTTTACGAGGGGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGcc                                          <  1:235500/298‑1 (MQ=255)
                                                                                                                                                   atttatttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGCCGCAGTCTGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTCCGTATTGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGACACTCCCCTGCTCGGTGCACAGCTGCTTCCCGGGGCAACCGACATCGCCCTGCCCGGCCCGCTGCCGCTCGTCCTCCCCCCCCCCCACACCAGTTCCCGGCCCACAACCCCCccgccgcgcgggcggca                                         >  2:504994/1‑287 (MQ=11)
                                                                                                                                                        tttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCCCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCCCCCGCCCCTACAGCCGTTCCCGGCCACAAACGCCCCCCCCCGTGGGGCGCCCCggac                                    >  2:502773/1‑296 (MQ=11)
                                                                                                                                                         ttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAACCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGccc                                   >  1:857399/1‑298 (MQ=255)
                                                                                                                                                          tCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcccc                                  >  1:687793/1‑298 (MQ=37)
                                                                                                                                                                  aTTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGaa                          >  1:435924/1‑298 (MQ=18)
                                                                                                                                                                   ttttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGCTCAGCCGCGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTCCCCCCGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGACCCCGACATCGCCCCCCCCCCCCCCCTGCCGTCCATCCCCCCCCGCACCTACAGCACTCACCCCACAGAAACCCCCGCCCCGGCGCGGCGCCCCGCCCCCGGCTgcgac                         >  2:252914/1‑294 (MQ=11)
                                                                                                                                                                   ttttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGaaa                         >  1:293425/1‑298 (MQ=18)
                                                                                                                                                                    tttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGAAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGaaaa                        <  1:442029/298‑1 (MQ=17)
                                                                                                                                                                      tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCGGGGCGTTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCGGGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCACCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATg                      <  1:779609/298‑1 (MQ=11)
                                                                                                                                                                         gATGGGCGCGAAACCGGCGGGGCGTAGGGGGGGCATTAGTCATTGTGGGGGGAGCATTGTCCAGGGTTCGGCCGGGGTACGTATTGGTGCCCCCGCCGGCGTGGGCTGGTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCAGCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCt                   <  2:566002/298‑1 (MQ=11)
                                                                                                                                                                          aTGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGTCTCGGCCCCGGCTGGAAAATGCCTg                  <  1:280720/298‑1 (MQ=14)
                                                                                                                                                                           tGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCCCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCCCCCGCCCCTCCACCCGCTACCGGACAAAACAGCCCGCGCCGGCGGGGGCCCCCGGCCCCGGCCGCAAAACGCCCGc                 >  2:870033/1‑298 (MQ=11)
                                                                                                                                                                            gaccgcaaaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCg                <  1:32095/292‑1 (MQ=9)
                                                                                                                                                                             aGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg               >  1:282495/1‑298 (MQ=14)
                                                                                                                                                                              gCGGAAAACGGGCGGCGCGTAAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGa              <  1:183117/298‑1 (MQ=9)
                                                                                                                                                                                 gAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat           >  1:429265/1‑298 (MQ=14)
                                                                                                                                                                                     aCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGc       >  1:780042/1‑298 (MQ=14)
                                                                                                                                                                                      ccGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCt      >  1:886533/1‑298 (MQ=14)
                                                                                                                                                                                      ccGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCt      <  1:745228/298‑1 (MQ=14)
                                                                                                                                                                                       cggcggCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCtt     <  1:1896/298‑1 (MQ=14)
                                                                                                                                                                                         gcggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGccc                                                                                                                        <  1:1185/181‑1 (MQ=18)
                                                                                                                                                                                         gcggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGccc                                                                                                                        >  2:1185/1‑181 (MQ=18)
                                                                                                                                                                                          gggcgCGTCAGTGCGGGAGGACGCAGTTTGGCGGTAGCATTGGCCTGGGTTCAGGCGGGGTACGTATTGGTGCCCCCACCGGCGTGGTCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa  <  2:221078/297‑1 (MQ=11)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                            
GCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTGTTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACA  >  CP006584/4304829‑4305312

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