Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,558,907 T→A K207M (AAG→ATG)  iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,558,9070TA96.0% 211.2 / NA 78K207M (AAG→ATG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base A (36/37);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558668‑4559152
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
tGCCGTCGGAGCGGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGGCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                                                                                                                                                                                                             <  1:178846/250‑1 (MQ=255)
      cGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTGCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                                                                                                                                                                                                       <  2:19532/250‑1 (MQ=255)
                   gatgcgagacAGCGGACCACCCGCATGAACGGCGTGGGAAAAAGCTTCGGCAGGGGCGACCGGATCGTCAATGATGCTGGTGATAGCCATGGTGCCGCCTTTGCGCGCGAACTACCGCGCCTGCTCAAACAACGGTACGT‑TCCGGTTAACGAGGGTCGG‑CAGCAGCTTGCAGCTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCAGTTTACTATCGCCCAtgtagtatacgg                                                                                                                                                                                                                            >  2:383252/10‑239 (MQ=255)
                   gACGCGTGCCAGAGGACTCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTGTCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                                                                                                                                                                                                         <  2:368204/251‑1 (MQ=255)
                       cgTGCCAGCGGAATACCCGACTGACCGGCGCGGGCACTTCCTTCGGCAGGACCGACCGGCACGTCAATGCTCCTGGCGATATCGATGGTGACGCCTTTGCGCGCGACCTCCAGCGCCTGCTCACACCAAGGTACGT‑TGCGGTTACCGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCACTTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                                                                                                                                                                                                     >  2:302633/1‑251 (MQ=255)
                          gCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                                                                                                                                                                                                                  <  1:310310/251‑1 (MQ=255)
                          gCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                                                                                                                                                                                                                  >  1:439855/1‑251 (MQ=255)
                             gggggaatggccGCCTGAGCGTCGCGGGCAATACCTTCGGCATGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCGTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACTT‑TGCGGTTAACGTGGGTGGG‑CAGCAGCTTGCTGATCGGCACATCGCAGGTTTCCAGCAGGTCATAAATTGGCTGTAGCGCCCTTTTACTGTCGCCCTTGTGGAACACGGTGACGCCAGGttt                                                                                                                                                                                                               <  1:205071/241‑1 (MQ=255)
                               cggaacgcccGCCTGAACGGCGCGGGCAATCCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTTGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                                                                                                                                                                                                             <  2:372572/244‑1 (MQ=255)
                                         gccggAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGAGGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa                                                                                                                                                                                                   <  2:109654/247‑1 (MQ=255)
                                           cTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc                                                                                                                                                                                                 >  2:124374/1‑251 (MQ=255)
                                                cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCa                                                                                                                                                                                            <  1:125913/251‑1 (MQ=255)
                                                 ggcgcgctcAATACCTTCGGCAGGGGCGACTGGTTCGTCACTGCTGCTGGCGATCTAGATGGCGCCGCCTTTGCGCGCGAACTCCACCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCCGATCGGCACATCGCAGTTTTCCCGCAGGTCATAAATACGCTGTCACGCCATTTTACTCTCGCCCATGTGGAACACGGTGCCGCTAGGTTTACCGCCGAGCAAACCGCCa                                                                                                                                                                                            >  2:168576/9‑250 (MQ=255)
                                                 ggCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGACAATGCTGCTGGTGATATCGATGGAGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCa                                                                                                                                                                                            >  1:202708/1‑250 (MQ=255)
                                                       ggCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCggg                                                                                                                                                                                     <  1:423814/251‑1 (MQ=255)
                                                       ggCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGAAGCCAGGTTTACCGCCGAGCAAACCGCCAACGCggg                                                                                                                                                                                     >  2:129590/1‑251 (MQ=255)
                                                           atacgttgggccggagCGACCGGTTCGTCAATGCTGCTGGTGATCTCTAGGGTGCCGCCTTTGCGCGCGACCTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAGCACGTTGACGCAAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTc                                                                                                                                                                                 <  2:86526/236‑1 (MQ=255)
                                                             cccctCGGCAGGGGCGCCCGGGCCGTGAAGGCTGCTGGTGATATCGAGGGTGCCGCCTGTGCGCGCGAACCCCCGGGCCTGCTCAAACAACGGTACGTGTG‑GGTTAACGTGGGTCGG‑CAGCAGCTTGCTGACCGGCACATCGCAGTTTCCCAGCAGGTCAAAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGATCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                                <  2:365245/246‑1 (MQ=255)
                                                                     cAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCat                                                                                                                                                                        >  1:26469/1‑250 (MQ=255)
                                                                         ggCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCCACGCGGGATTCTGCCGCCCTATTgg                                                                                                                                                                   >  2:376301/1‑251 (MQ=255)
                                                                           cGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                 >  1:148695/1‑251 (MQ=255)
                                                                            gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                 >  1:332323/1‑250 (MQ=255)
                                                                            gACCGGATCGTCAATGCTGCTGGTGATATCGCAGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCCTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGGTTACCGCCGAGCAAACCGCCACCGCGGGATTCGGCCGCCatctgggcc                                                                                                                                                                 >  2:132997/1‑243 (MQ=255)
                                                                              ccGGTTCGTCAATGCTGCTGGTGATATCGAGGGTGCCGCTTTTGCGCGCGAACTCCAGCGCCTGCTCAACCAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTAATAAATAGGCTGTACCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAgg                                                                                                                                                              <  1:310892/251‑1 (MQ=255)
                                                                                ggTTCGTCAATGCTGCTGTGGATATCGATGGGGCCGCCTTTGCGCGCGAACGCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACAGCGCAGTTTTCCAGCAGGTCATAAATAGGCTTTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTg                                                                                                                                                            <  1:267773/251‑1 (MQ=255)
                                                                                    cGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATa                                                                                                                                                         >  1:180076/1‑250 (MQ=255)
                                                                                       cAATGCTGCTGGTGATATCGATGTTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑AAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCAAGGTGATaaaa                                                                                                                                                      <  1:13920/250‑2 (MQ=255)
                                                                                           gctgctGTTGATATCGATGTTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAGGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCg                                                                                                                                                 <  1:416189/251‑1 (MQ=255)
                                                                                           gctgcgggggatggcagtggggcacccTTTGCGCGCAACCCCCAGCTCCTGCTCAAACACCGTTACGT‑TGCGGTTAACG‑GGGTGGGGCACCACTTTGCTGATGGGAACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGCCGCCAATGTGGAACACGGGGACGCCAGTTTTCCCGCGGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTcc                                                                                                                                                  <  2:385019/225‑1 (MQ=255)
                                                                                            ctgctgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGTTTCTGCCGCCATATTGGCCCGGTGATAAACGTCCgg                                                                                                                                                >  2:151852/1‑251 (MQ=255)
                                                                                               ggggTGATATCGAGGGGGCCGTCTTTGCGCGCGAACTCCAGCGCCTTCTCAAACAACGGTAGGT‑GGCGGTAAAGGGGGGCCGG‑GAGAAGCTTGCTGGTCGGCACATCGCAGTTTTCCAGCAGGTCAAAAATAGGCGGTAACGCCATTTGACTGTCGCCCGTGTGGAACCCGTGGACGCCAGGTTTCCCGCCGAGAAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGc                                                                                                                                             <  2:101952/249‑1 (MQ=255)
                                                                                               ctgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGc                                                                                                                                             >  1:30762/1‑251 (MQ=255)
                                                                                                tgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCg                                                                                                                                            >  1:398657/1‑251 (MQ=255)
                                                                                                  gggaaaggtagggTGCCGCCTTTGCGGGCGAACCCCAGCGCCTGCGCAAACACCGTGACTT‑TGCGGTTAACGTGGGGCGG‑CATCAGCGTGTTGTTCGGCACTTCGCAGTTTTCCAGCAGGTCAGAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGACCACGGTGACGCAAGGTTGCCCGCCGAGCAAACCGCCAACGCGTGCTCCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGc                                                                                                                                          <  2:450094/240‑1 (MQ=255)
                                                                                                   ggaatacGATGGGGCCTCATTTGCGGGCGAACTCCAGGGCCGGCTAAAACAGCGGTAGGT‑TGCGGTAACCGGGGGTCGG‑AAGCAGCTTGCGGATCGGCACATCGCAGTTTTCCAGAAGGTCAAAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACAGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCa                                                                                                                                         <  2:91936/245‑1 (MQ=255)
                                                                                                    gATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAGCGCCATTTTACGGTCGCCCATGTGAAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAg                                                                                                                                        <  1:459829/251‑1 (MQ=255)
                                                                                                      tatCGATGATGACGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACCACGGTACGT‑TGCGGCTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAACTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCCGGTTTACCGCCGAGCAACCCGCCAACGCGGGCTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTCCGGCAGa                                                                                                                                       >  2:180311/1‑250 (MQ=255)
                                                                                                        tCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACg                                                                                                                                    >  1:20136/1‑251 (MQ=255)
                                                                                                         cGATGGTGCCGCCTTTGCGCGCGAACGCCGGCGCCGGGTCAAACAACGGTACGT‑TGCGGTTAACGGGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTAAAAAAGAGGCTGTAACGCCATTTTACGGTCGCCCATGGGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCAAAGGCGGGATTCTGCCGCCATATTGGACAGGTGATAAACGTCCGGTGCGGCAGAACg                                                                                                                                    <  1:105619/250‑1 (MQ=255)
                                                                                                            gggTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGGGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGa                                                                                                                                <  1:48922/250‑1 (MQ=255)
                                                                                                             ggTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAAGGTGGGTCGG‑CAGCATCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAAGAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAt                                                                                                                               <  1:60043/251‑1 (MQ=255)
                                                                                                               tGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGCCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                              >  1:300206/1‑250 (MQ=255)
                                                                                                                 ccgccTTTGCGAACGCACTCCAGAGCCTGCTCAAACAACGGTACAT‑TGCGGTTAAGGTGGGTCGG‑CCGACGCTTACGGATCGGAACATCGCAGTTTTCCAGCAGGTCATAAATAGACTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTAAAGCAAGGTTTACCACTGATAACACCGCCAACGCGGGTTTATGCAGCCATATTGACCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                              >  2:61928/1‑248 (MQ=255)
                                                                                                                       ttGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGA‑TGCGATACACGTGGGTCGG‑CAGCCGCTTGCTGATCGGCACATTGAAGTTTACCAGCAGGTCATAAATAGGCAGTCACGACATTTTACTGTCGCCCATGTGGAACCCGGAGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGAGATCAACGTCCGGGGCGGCAGAACGGTGCACAGAGATGgc                                                                                                                     >  1:294266/1‑251 (MQ=255)
                                                                                                                                  cTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCCTATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGc                                                                                                          >  2:405339/1‑251 (MQ=255)
                                                                                                                                    ccatcgccggcacaaaaaccggcaCGT‑TGCGGTGAACGTGGGGCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGGCATAAAGAGGCTGTAACGCCATTTTACTGTCCCCCAGGGGGAACACGGTGACGCAAGTTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCa                                                                                                        <  2:344452/228‑1 (MQ=255)
                                                                                                                                    ccAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGccc                                                                                                        >  2:247708/1‑250 (MQ=255)
                                                                                                                                       gcgcCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                     >  1:168681/1‑251 (MQ=255)
                                                                                                                                       gcgcCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                     <  1:129590/251‑1 (MQ=255)
                                                                                                                                                cAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CCGCAGCTTGCTGATCGGCACATCGCCGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcg                                                                                            >  2:176611/1‑251 (MQ=255)
                                                                                                                                                     aacGGTACGA‑TGAGGTTACCGTGGGTCGG‑AAGCTGCTTGCTGATCGGCAAATCGCAGTTTTCCAGAAGGTCATCAATAGGCTGTACAGCCATTTTTCTGTAGCCCATGTGGAACACGGTGACGCCAGGTTTACCGACGAGCCAACCGCCAACGCGGGCTTCTGCCGCCCTATTGGCCAGGTGATAAACGTACGGAGCGGCAGAACGGGGATCAGAGATGGCGAATTTAACGCCCATCATACGATCGGTAAt                                                                                       >  2:216125/1‑251 (MQ=255)
                                                                                                                                                      acGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACCTCGCAGTTTTCCAGCAGGACATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGGGATCAGAGATGGCACATTTCAAGCCAATCACACGATCGATAAt                                                                                       >  1:431269/1‑250 (MQ=255)
                                                                                                                                                      acGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGGGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATc                                                                                      >  2:314191/1‑251 (MQ=255)
                                                                                                                                                       cGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTTCTGATCGGCACATCGCAGTTTTCAAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCCTGTGGAGCACGGTGACGCCAGGGTTACCGCCGAGCAAACCGCCAACTCGGGTGTCTGCCGCCATATTGGCCAGGCGATATGCGTCCGGTACGGCGGCACGGTGATCATAGCTGGCGCAGTTCACGCCAATCACACGATCGATATTCg                                                                                     >  2:87080/1‑251 (MQ=255)
                                                                                                                                                         gTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGCACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGcc                                                                                   >  2:232702/1‑251 (MQ=255)
                                                                                                                                                          caCGT‑GGCGGTTAGCGGGGTTCGG‑CAGAGGCGTGCGGAGCGGCACCTCGCAGTTTCCAAGCAGTAAATAAAGAGGCGGAACCGCCATTTTACGGTCGCCCGTGTGGAACAGGGTGACGCCAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCa                                                                                  <  2:384817/250‑1 (MQ=255)
                                                                                                                                                            cGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTGGCTGATCGGCCCTGCTCAGTTTTCCGGCAGGTCATAAATAGGCTGTAACGCCATTTGACGGTCGCCAATGTGAAACACGTGGACGCCAGTTTGACCGCCGAGCAACCCGCCAAGGCGGGATTCTGCCGCCAGATTGTGCAGGTGTGAAAGGTCCGGGGCGGCAGAACGGTGATCAGAGAGGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACg                                                                                <  1:292944/251‑1 (MQ=255)
                                                                                                                                                             ggtggCGGTAAACGTGGGTGGG‑CGGCGGTTGGCGGATCGGCACATCGCAGTTTTCCAGCGGGTCATAAAGAGGCTGTACCCCCATTTTACGGTCGCCCATGGGGAACACGGTGGCGCCAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACg                                                                                <  2:148695/247‑1 (MQ=255)
                                                                                                                                                                tGCGGTTACCGGGGGCCGG‑CAGCAGCTTGCGGATCGGCCCACAGCAGTTTTCCCGCCGGTCATAAATAGGCTGTAACGCCATTTACCTGCCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACGGCCAACGCGGGATTCGGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTc                                                                              <  1:325598/250‑1 (MQ=255)
                                                                                                                                                                   ggTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                          <  1:159927/251‑1 (MQ=255)
                                                                                                                                                                             ggTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAc                                                                <  1:405339/251‑1 (MQ=255)
                                                                                                                                                                                        agcTTGCTGATCGGCCCATCGCAGTTTCCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGt                                                      <  1:27840/251‑1 (MQ=255)
                                                                                                                                                                                          cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGc                                                    >  1:335737/1‑251 (MQ=255)
                                                                                                                                                                                           ttGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCg                                                   >  1:168577/1‑251 (MQ=255)
                                                                                                                                                                                             gCTGATCGGCACATCGCAGTTTCCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCggg                                                 <  1:151852/251‑1 (MQ=255)
                                                                                                                                                                                                gATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGGGCGGCCGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACCCGATCGCTAATCGCCACGTGTTTTTCCACGGAACCCGTAATGGTGCGGGa                                                >  2:364521/1‑249 (MQ=255)
                                                                                                                                                                                                  ccGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTACCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAAGCGCCACGTCTTTTTCCACTGAACCCGGAATGGTGCGGGAAggg                                            <  1:180311/250‑1 (MQ=255)
                                                                                                                                                                                                            gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAg                                  <  1:124374/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGCGGGTCATAAATAGGCTGTAACGCCTTTTTCCTGTCGCCCATGTGGAACACGTTGACGCCGGGTTTACCTCCGAGCAAACCGCAAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                          <  1:116671/251‑1 (MQ=255)
                                                                                                                                                                                                                      agcagGTTATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTGTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAg                        <  1:243066/251‑1 (MQ=255)
                                                                                                                                                                                                                             cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAg                 <  1:247708/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                >  2:109194/1‑251 (MQ=255)
                                                                                                                                                                                                                               tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGAGCGGGCAGGGACACGATAAGCGCCGGTCAGCATCCTGgc               >  2:369297/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gCTGTAACGCCATTTTACTGTCGCCCATGTGGCACACGGTGACGCCAGGTTTACAGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATACTCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGCAGGGCCATGATAAGCGCAGGTCAGCATCCAGGCGCTGagg        >  2:177413/1‑249 (MQ=255)
                                                                                                                                                                                                                                      gCTGTAACGCCATTTTACTGTCGCCCATGTGGAACCCGGTGCCGCCAGTTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATg        <  1:394138/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCTGTAACCCCTTTTTACTGGCCCCCTGGTGGACCCCGGGTACGCCAGGTTTACCGCCGAGCAAACCGCCAACGGGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGTACATGATAAGCGCCGTTCATCATCCAGGCGCTGATg        <  2:439855/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aCGCCTTTTTCCTTTCGCCCAGGTAGAACACGGTAGCGCCAGGTTTACCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCTGTGCGGCAGAACGGTGATCAGAGATGTCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCATGCGGTGATGCcttct  <  2:233167/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGCCATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttct  >  2:130896/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
TGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGT‑TGCGGTTAACGTGGGTCGG‑CAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558668‑4559152

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