Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 597,052 T→G 5.0% intergenic (+79/‑79) cusC → / → cusF copper/silver efflux system, outer membrane component/periplasmic copper‑ and silver‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913597,0520TG 5.0% 209.9 / 2.1 68intergenic (+79/‑79)cusC/cusFcopper/silver efflux system, outer membrane component/periplasmic copper‑ and silver‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base G (2/3);  total (32/36)
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.98e-01

GCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCT  >  NC_000913/596756‑597347
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                       
gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTg                                                                                                                                                                                                                                                                                                     >  2:585203/1‑301 (MQ=255)
  ttGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTggg                                                                                                                                                                                                                                                                                                   <  2:632996/300‑1 (MQ=255)
  cTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTggg                                                                                                                                                                                                                                                                                                   <  1:55619/301‑1 (MQ=255)
    taaCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGa                                                                                                                                                                                                                                                                                                 <  2:479249/300‑1 (MQ=255)
        gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGGAGCTTACATAAGTAAGTGACTggggtgagc                                                                                                                                                                                                                                                                                             <  2:380935/301‑6 (MQ=255)
              aTTAGCGCCCAGGAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCt                                                                                                                                                                                                                                                                                       <  1:528498/301‑1 (MQ=255)
              aTTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCt                                                                                                                                                                                                                                                                                       <  2:108388/301‑1 (MQ=255)
                 aGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGt                                                                                                                                                                                                                                                                                    >  2:620690/1‑301 (MQ=255)
                    gcgccgccgcGTTATCTGGCATCGCTTCAAATTTCTTTGCATCGGGCGCGGGCAGTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTGTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGc                                                                                                                                                                                                                                                                                 <  1:483339/293‑1 (MQ=255)
                     cccAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGcc                                                                                                                                                                                                                                                                                >  2:545733/1‑301 (MQ=255)
                         gcagcGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTcc                                                                                                                                                                                                                                                                            <  1:21283/301‑1 (MQ=255)
                            gcGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGc                                                                                                                                                                                                                                                                         <  2:574983/301‑1 (MQ=255)
                            gcGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGc                                                                                                                                                                                                                                                                         <  2:154789/301‑1 (MQ=255)
                               ttATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAAc                                                                                                                                                                                                                                                                      >  2:92028/1‑301 (MQ=255)
                                  tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCg                                                                                                                                                                                                                                                                   <  2:24038/301‑1 (MQ=255)
                                      gcttCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAatt                                                                                                                                                                                                                                                               <  1:91904/298‑1 (MQ=255)
                                         tCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGACAGCTTAACGAAATTTCTTTGTATACCGCACTTGTTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCCGTCACTTACTTATGTAAGCCCCTGGGGATTCACTCGCTTGTCCCCTTCCTGCCAccagcactatt                                                                                                                                                                                                                                                            >  1:399693/1‑291 (MQ=255)
                                                            ataccctttaccatgccagtggcgctaataacctgtggttgtgcttcgctcatggtttcatgatgatgttcgttagcctgggcattaaagccaataacggtaaacagactgaacattgcgacttgcagtgcttttttcatttattctctcctggagttaaaagtcATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACt                                                                                                                                                                                                                                         >  1:54543‑M1/166‑301 (MQ=255)
                                                            ataccctttaccacgccagtggcgctaataacctgtggttgtgcttcgctcatggtttcatgatgatgttcgttagcctgggcattaaagccaataacggtaaacagactgaacattgcgacttgcagtgcttttttcatttattctctcctggagttaaaagtcATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAaca                                                                                                                                                                                                                                         <  2:54611‑M1/136‑2 (MQ=255)
                                                                         atgatatCAGCACGGGGCAGTAAGTTATTTGGAAGTGCTGGATGCCGAGCTTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACGTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATa                                                                                                                                                                                                                            <  1:137008/298‑1 (MQ=255)
                                                                              ggggctctaataagttgtggttgtgcttcgctcatggtttcatgatgatgttcgttagcctgggcattaaagccaataacggtaaacagactgaacattgcgacttgcagtgcttttttcatttattctctcctggagttaaaagtcATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGa                                                                                                                                                                                                                       >  1:57204‑M1/148‑301 (MQ=255)
                                                                              atCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGa                                                                                                                                                                                                                       >  2:510874/1‑301 (MQ=255)
                                                                               tCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGaa                                                                                                                                                                                                                      >  2:523524/1‑301 (MQ=255)
                                                                                 aGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGaaaa                                                                                                                                                                                                                    <  2:297044/301‑1 (MQ=255)
                                                                                  gCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAAAATAAATGaaaaa                                                                                                                                                                                                                   >  1:7556/1‑301 (MQ=255)
                                                                                          cAGTAAGTTATCTGGAAGTGCTTGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGc                                                                                                                                                                                                           <  2:456034/301‑1 (MQ=255)
                                                                                                                    gCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGAACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTAcc                                                                                                                                                                                 <  1:511608/301‑1 (MQ=255)
                                                                                                                              tttttttTGCAACCCTATAAATTTTACTTGATCTTATTTATGCCCGTCAGTGTAACGATATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTGCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCGCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCt                                                                                                                                                                       <  1:545060/300‑1 (MQ=37)
                                                                                                                              ctttccctacacgtcgctcttcctatcTTGTTTTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCATTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCt                                                                                                                                                                       <  1:379008/275‑1 (MQ=25)
                                                                                                                              cTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCt                                                                                                                                                                       >  2:86629/1‑301 (MQ=255)
                                                                                                                                        aaCCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATACACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAAAATAAATGAAAAAAGAACTGCAAGTCGCAATGTTAAGTCTGTTAACAGTTATAGTCTTAAATGcaca                                                                                                                                                             >  2:503269/1‑301 (MQ=255)
                                                                                                                                         aCCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCCAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCCCTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGCCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATCCCCAg                                                                                                                                                            >  1:110455/1‑301 (MQ=255)
                                                                                                                                         aCCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCaa                                                                                                                                                            >  2:445990/1‑300 (MQ=255)
                                                                                                                                            cGACAAACTTTTCTTGTTCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCt                                                                                                                                                         <  1:582795/301‑1 (MQ=255)
                                                                                                                                                          tGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCtagatcggaagagcacacgtctgaactcc                                                                                                                                           >  2:379504/1‑273 (MQ=255)
                                                                                                                                                          ctctttccctacacgacgctcttccgatctaaTTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAAcatcatca                                                                                                                                           <  1:464203/271‑1 (MQ=255)
                                                                                                                                                            aTCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCCGGAGAGAATAAATGAAAACAGCACTGCAAGTCGAAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCAAGGCAAACGAAcatcatcatt                                                                                                                                         >  1:475051/1‑300 (MQ=255)
                                                                                                                                                                aagttATGCCCGTCAGTTTACAGAAATTTTTTGGTATACCGAACTTGGTGGCGGTTGAAAGCAATATACTCTTAATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTAACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAAGGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAAc                                                                                                                                     <  1:619832/298‑1 (MQ=37)
                                                                                                                                                                  aTTTTGCCCGTCAGGTTAACGAAAGTTTTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCa                                                                                                                                   <  1:445426/301‑1 (MQ=255)
                                                                                                                                                                    tATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATg                                                                                                                                 <  1:522852/301‑1 (MQ=255)
                                                                                                                                                                     tttcccGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGa                                                                                                                                <  2:7556/298‑1 (MQ=255)
                                                                                                                                                                     aTGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATACTCATGCACCCATGa                                                                                                                                >  1:573300/1‑301 (MQ=255)
                                                                                                                                                                       gCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGCAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCACAGGCTAACGAACAACATAATAAAACAATGAGc                                                                                                                              >  1:347669/1‑301 (MQ=255)
                                                                                                                                                                                      gatgttttgtggggtgcggcatgtggtggGGGTTGGGAGGAAGATACGCTGCATATGTCAGGGTGCAGGTGCTGCGGCTGCGTGCGCTCACCCCGGTGACTTACGGATGTAAGCTCCTGGGGATTCACGCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTGTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAg                                                                                                               <  1:584479/278‑1 (MQ=25)
                                                                                                                                                                                        aaTTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCAAGatcggaagagcacacgtctgaactccaa                                                                                                             >  2:464767/1‑274 (MQ=255)
                                                                                                                                                                                           ttCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGttat                                                                                                          >  2:350103/1‑301 (MQ=255)
                                                                                                                                                                                            tCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTACGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCAAGAGCAAAACAAAACCACAGGttaat                                                                                                         >  1:597230/1‑299 (MQ=255)
                                                                                                                                                                                             ttttGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGttatta                                                                                                        <  1:349640/300‑1 (MQ=255)
                                                                                                                                                                                                          aCTTGGTGGTGGTTGGCTGCACTTCATTTGTCATACTTTAAGTTTCATGTGCTGCTGCTGCGTTCGCTCACCCCAGTCACTTTCTTAGTGAAGTGCCTGGGGATTTACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACGTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGt                                                                                           <  1:502597/301‑1 (MQ=37)
                                                                                                                                                                                                             ttcagacgtgtgctcttccgatctTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGt                                                                                        >  2:606642‑M1/25‑301 (MQ=255)
                                                                                                                                                                                                                        gACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATca                                                                             >  2:610922/1‑300 (MQ=255)
                                                                                                                                                                                                                          cAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGaa                                                                           >  2:77126/1‑300 (MQ=255)
                                                                                                                                                                                                                                  TACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAg                                                                   >  2:400207‑M1/1‑301 (MQ=255)
                                                                                                                                                                                                                                   ACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAaga                                                                  <  2:470510‑M1/301‑2 (MQ=255)
                                                                                                                                                                                                                                     TCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGgttattagcgccactggcgtggtagatcggaagagcgacgtgtaggg                                                                <  1:605881‑M1/301‑24 (MQ=255)
                                                                                                                                                                                                                                           tACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCaaaaaaaa                                                          >  2:134800/1‑300 (MQ=255)
                                                                                                                                                                                                                                           aaCTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAAt                                                          <  1:77036/300‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGAGGTAAAGGGTATCGAACTGGAAAGAAAAAAAATCAccc                                                     >  1:358853/1‑300 (MQ=255)
                                                                                                                                                                                                                                                   gTTGCATGTGCTTCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCAccatcc                                                  <  1:86505/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tttgcgtTCGCTCACCCCATTCACTTACTTATGTAAGCTCCTGGGGATTCCCTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATtg                                   <  1:469909/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGAAtgatgc                               >  2:16683/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           gCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGa                          >  2:35580/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                  ccaggagcttacataAGTAAGTGACTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGcc                   >  2:280587/8‑301 (MQ=255)
                                                                                                                                                                                                                                                                                  ccaggagcttacataAGTAAGTGACTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGcc                   >  2:552306/8‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        tttccgatctGTAATCTCTGGGGGATTCACTCGCTTTTCGCCGTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGAt             <  1:342040/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             tttaTGTAAGTTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCa        <  2:359337/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  gTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGgaaaaagatcggaaa   >  2:342503/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                                   tAAGCTTCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCt  <  2:574015/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                       
GCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTTTATTTGCAACCCGACAAACTTTACTTGATCTGAATTATGCCCGTCAGGTTAACGAAATTTCTTTGTATACCGCACTTGGTGGCGGTTGACAGCAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTAAAGGGTATCGATCTGGAAAGCAAAAAAATCACCATCCATCACGATCCGATTGCTGCCGTGAACTGGCCGGAGATGACCATGCGCT  >  NC_000913/596756‑597347

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