Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,620,877 T→C 3.6% intergenic (+50/‑7) deoB → / → deoD phosphopentomutase/purine nucleoside phosphorylase 1, nicotinamide 1‑beta‑D‑riboside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,620,8770TC 3.6% 352.1 / 3.3 113intergenic (+50/‑7)deoB/deoDphosphopentomutase/purine nucleoside phosphorylase 1, nicotinamide 1‑beta‑D‑riboside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (57/50);  new base C (2/2);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-02

Only 100 of 116 total aligned reads displayed.

ACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCA  >  NC_000913/4620643‑4621106
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
aCCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATaaaa                                                                                                                                                                                                                                     >  2:904644/1‑239 (MQ=255)
          cGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGTTAAAACAATGGcga                                                                                                                                                                                                                            >  2:724958/1‑236 (MQ=255)
            aTCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAGCCGGGCTCACTGGTTCATCGTGAAACCTTCGCGGACATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGTTATGGAATATGGCAAAGCCAGTTTCTGTTGGATTTGGGCGGAGGGTGGGTTCCGCCTTTGTTATGTCACAAAAAGGaaaaaaaacaggtaaccc                                                                                                                                                                                                                         >  1:411478/1‑227 (MQ=255)
            aTCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCCAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGTTGGATTTGGCCGGAGCGTTGACTCCGCCTTTGTTATGTCCCAAAAAAGATAAAAAAATgggcacc                                                                                                                                                                                                                          >  1:43557/1‑233 (MQ=255)
              ccGGCCTGGCCCGGTATTTGCCACACGGGTGAACAAATTCCGGTACTGGTTTATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCTAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCGTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCca                                                                                                                                                                                                                       <  2:1024437/239‑1 (MQ=255)
              ccGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCca                                                                                                                                                                                                                       <  2:606068/239‑1 (MQ=255)
               cGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGGGAAACCTTCGCGGATGTCGGCCAGACTCTGGCAAAACATTTTGGTACTTCTGATATGGAATATGGCAAAACCATGTTCTGTTGGATTTGGTCGGTGCGTTGACTCCGCCTTTGTTATGTCAAAAAAAGGAAAAAACAAGGGGTACCCaac                                                                                                                                                                                                                      >  1:592105/1‑239 (MQ=255)
               cGACCTGGAACGGTACTGAACACCCGGGTGAACCCATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCca                                                                                                                                                                                                                       <  2:592105/238‑1 (MQ=255)
                         cGGTACTGACCCAGCGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATg                                                                                                                                                                                                            <  2:129250/238‑1 (MQ=255)
                           gTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGGGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCca                                                                                                                                                                                                                       >  1:869334/1‑226 (MQ=255)
                           gTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCca                                                                                                                                                                                                                       <  2:869334/226‑1 (MQ=255)
                                 gccaccgggcggaaccacttCCGGTACTGGGATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATgg                                                                                                                                                                                                   <  2:197307/221‑1 (MQ=255)
                                 aCCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGCATCCGCCTTTGTTATGTACCAAAAAGGATAAACCAATGGCTACCCACACCATTAATGAAGAAATgg                                                                                                                                                                                                   >  1:522493/1‑239 (MQ=255)
                                 aCCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACCAATTAATGCAGaaagtg                                                                                                                                                                                                   >  2:370333/1‑236 (MQ=255)
                                           tGAACCCATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAAt                                                                                                                                                                                                                                 <  2:540561/200‑1 (MQ=255)
                                           tGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGCTATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAAt                                                                                                                                                                                                                                 >  1:540561/1‑200 (MQ=255)
                                           tGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGAt                                                                                                                                                                                                                                         >  2:512485/1‑192 (MQ=255)
                                           tGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACAt                                                                                                                                                                                                                 <  1:975695/215‑1 (MQ=255)
                                           tGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACAt                                                                                                                                                                                                                 >  2:975695/1‑215 (MQ=255)
                                           tGAAAACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGAt                                                                                                                                                                                                                                         <  1:512485/192‑1 (MQ=255)
                                              acacATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTg                                                                                                                                                                                       <  2:690312/238‑1 (MQ=255)
                                                acaTTCCGGTACTGGGATATGGCCCCAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACg                                                                                                                                                                                    <  2:43557/239‑1 (MQ=255)
                                                        gTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAAAGGCTACCCCCCACATTAATGCAGAAAGGGGCGATTTCCCTGAAGAAGgtttt                                                                                                                                                                            >  1:370652/1‑238 (MQ=255)
                                                            tGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACCTTAATGCAGAAATGGGCGCTTTCGCTAACGTAGGTTTGATg                                                                                                                                                                         >  2:146069/1‑238 (MQ=255)
                                                            tGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGCACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCAAAAAAAGGATAAAAAAATGGTGACCCCCCACATTAAGTCAGAAATGGGCGAtttgcctacggtagattttgtgc                                                                                                                                                                        >  1:36414/1‑219 (MQ=255)
                                                              gTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCCCAAAAAGGTTAAAACAATGGCTACCCCCCACATTAATTCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATCCCa                                                                                                                                                                      >  1:707352/1‑239 (MQ=255)
                                                              gTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTGATGCCAg                                                                                                                                                                     >  1:906100/1‑239 (MQ=255)
                                                              gTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGGAGTTTTGATGCCa                                                                                                                                                                      >  2:789036/1‑239 (MQ=255)
                                                                        cGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCCAAATATTTTGGTACTTCTGATATGGAATATGGCCCAGCCCTGTTCTGTTGGCTTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACCAACAGGATAAAACATTGGCAACCCCCCCCTTTAATCCAGAACTGGCCGCTTTTCCCAACGTACTTTTTGTGCCCggtcccccc                                                                                                                                                             >  1:599553/1‑231 (MQ=255)
                                                                        cGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGAccc                                                                                                                                                             >  2:963358/1‑237 (MQ=255)
                                                                        cGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGc                                                                                                                                                           >  2:877760/1‑239 (MQ=255)
                                                                        cGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCAGCCTTTGTTATGTCCCAAAAAGGATAAAACAATGGCTACCCCACCCATTAATTCACAAATGGGCGATTTAGCAGACGTAGTTTTTATGCCCGGCGcccccc                                                                                                                                                           >  1:534579/1‑237 (MQ=255)
                                                                            aGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTGATGCCAGGCGACCCGCtgcgt                                                                                                                                                      <  2:906100/239‑1 (MQ=255)
                                                                              taaacccgggccacttgggttcacGTGAAACCTTCGCGGATAGCGGCCAGACTCTGGCAAAATTTTTTGGTTCTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgt                                                                                                                                                      <  1:485881/215‑1 (MQ=255)
                                                                                aaaCCGGGTTCCTTGGGGTATCGTGAAACCTTCGCGTATATCGGCCAGACTTGGGCAAAATTTTTTGGTTATTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgtgcg                                                                                                                                                   <  1:789036/239‑1 (MQ=255)
                                                                                aaaCCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgtgcg                                                                                                                                                   <  2:370652/239‑1 (MQ=255)
                                                                                aaaCCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgtgcg                                                                                                                                                   <  2:36414/239‑1 (MQ=255)
                                                                                       gCTCACTGGGTAATTGTGAAACCTTCGGCGATATCGGCCGATCTATGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGTGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGtat                                                                                                                                             <  2:910879/238‑1 (MQ=255)
                                                                                        cTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCCAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGAGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCCCACATTAATGCGGAAATGGGGGATTTCGCCGACGTAGGTTTGAAGCCGGGGGACCCCCGGGGGGCGAAggaaaa                                                                                                                                           >  1:549671/1‑234 (MQ=255)
                                                                                        cTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATCAAACAATGGCTACCCCCCACATTAATGCAGAAATGGGCGATTTAGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTatat                                                                                                                                           >  2:485881/1‑239 (MQ=255)
                                                                                                     cGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCACCACATTAATGCAGAAATGGGCGGTTTCGCTGACGTAGTTTTGATTCCAGTGGAACCGCTGCGTGCGACGCAAATTGCTGCAACtt                                                                                                                                >  2:709909/1‑237 (MQ=255)
                                                                                                          aaCCTTCGGTGTTGTTGGACAGACCTTGGAAAAATTTTTTGGTTATTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAAGGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGa                                                                                                                         <  1:370333/239‑1 (MQ=255)
                                                                                                             cTTCGCGGATATCGGCCAGACTCTGGCAAAATTTTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGa                                                                                                                      <  2:828591/239‑1 (MQ=255)
                                                                                                             cTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGa                                                                                                                      <  1:1013894/239‑1 (MQ=255)
                                                                                                             cTTCGCAGATATCGGCCAGACTCTGGCAAAATATTTGGGTAATGCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGGTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGa                                                                                                                      <  2:599553/239‑1 (MQ=255)
                                                                                                                      tatCGGCCAGACTCTGGCAAAATTTTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCgtg                                                                                                              <  2:337583/238‑1 (MQ=255)
                                                                                                                       atccgacAGACGCTGGCAAAATATTTTGGTACTTCTGATCTGGAATATGGCAAAGACATGTTCTGATCGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTACTGCACAAATGGCCGATTTCCCTCCCGTACTTTTAATGCAGGCTGACCCCGGACCTGCCACGTATATTGTTCAAAATTTCCTTTaaaaaacccgcgag                                                                                                            >  1:987638/7‑227 (MQ=255)
                                                                                                                       atCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGTCGACCCGCTGCGGGCGAAGAATATTGCTGAAACTTTCCTTGAAGATGCCCgtgaa                                                                                                            >  2:487960/1‑239 (MQ=255)
                                                                                                                             cAGTCTTTGGCAAAATATTTTTGTTCTTCTGATATGGAATATGGCAAAGCCTTTTTTTGATGGATTTGGGCGGAGCGTTTACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACACCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGaac                                                                                                      <  1:69549/239‑1 (MQ=255)
                                                                                                                             cAGACTCTGGCAAAATATTTTGGTTCTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGaac                                                                                                      <  2:258873/239‑1 (MQ=255)
                                                                                                                               gACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACCAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTCCGTGCGAAGAATATTGCTGAAACTTTCCTTGAACATGCCCGTGAAGTGaaca                                                                                                     >  2:229600/1‑238 (MQ=255)
                                                                                                                                    tGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTc                                                                                               >  2:101969/1‑239 (MQ=255)
                                                                                                                                          aaTATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAATGGAACAACGTTCGCGGTa                                                                                         >  2:606923/1‑239 (MQ=255)
                                                                                                                                                      cTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGACCAACGTTCGCGGAATGCTGGGCTTc                                                                              >  2:650180/1‑238 (MQ=255)
                                                                                                                                                       ttgtgtttgtgattttGGCGAAGTCATTGTCTGATGAAGTTGGGGGGTGGTTTGACCTCGCCTTTGTTAGTTCACAAATAGGATAAAACAATGGCTCACCCACACAGTAATTCAGAAATGTGCGATGTCGCTGACGTAGTTTTGATGCCATGCGACCCGCTGCGTGTGAAGTATATTGCTGAACCTTTCCTTGAAGATGCTCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCAc                                                                            <  1:20529/224‑1 (MQ=255)
                                                                                                                                                             tatGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGt                                                                        >  2:636867/1‑237 (MQ=255)
                                                                                                                                                                    tatGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTATCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAAAGTTCCCGGTATGAGGGGCTTCACggggaatgaaaaa                                                               >  1:835520/1‑229 (MQ=255)
                                                                                                                                                                      tGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGc                                                                                                                                                                                                                              <  2:585456/80‑1 (MQ=255)
                                                                                                                                                                      tGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGc                                                                                                                                                                                                                              >  1:585456/1‑80 (MQ=255)
                                                                                                                                                                       ggCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGCCGACCCGCTGCGTGCGCAGGATATTGCTGAAACTTTCCTTGAAAATGCCCGTGAAGTGAACAACGTTGGCGGTATGCTGGGCTTCCCCGGTACTTACAAAgg                                                             >  2:173057/1‑238 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCg                                                                                                                                                                                          <  2:776932/109‑1 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCg                                                                                                                                                                                          >  1:776932/1‑109 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGCACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCaa                                                       >  2:813023/1‑239 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGCACAACGTTCCGGGTCTACTGGGCTTCACCGGTACTTACAACAGCCGCaa                                                       >  1:48261/1‑239 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTcgcg                                                                                            <  1:1021737/202‑1 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTcgcg                                                                                            >  2:1021737/1‑202 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACATTCCCGGAATCCGGGCCTTAACCCGAACTTACAACAGCCGCaa                                                       >  1:214580/1‑239 (MQ=255)
                                                                                                                                                                            aGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGGTAACAAAGTTCGCGGTATGCTGGGCTTCACCGAGACTTACAAAGGCCCCaa                                                       >  1:616144/1‑239 (MQ=255)
                                                                                                                                                                             ttcccttttttgtGGGTTTGGGGGTGAACGTTGACCCCGCCTTTGTTGTGGCAAATAAGGGATAAAGCAATGGCTACCCCGCCAATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTGGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTGCCTTGAAGTTGCACGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCCCCGGTACTTACAAAGGCCGCaa                                                       <  1:173057/227‑1 (MQ=255)
                                                                                                                                                                                        tggggtttGGGCGGGGCTTTGCCTCCGTCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAAt                                           <  1:487960/234‑1 (MQ=255)
                                                                                                                                                                                        aTGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCg                                                                                                                                                                                          >  2:358067/1‑97 (MQ=255)
                                                                                                                                                                                        aTGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCg                                                                                                                                                                                          <  1:358067/97‑1 (MQ=255)
                                                                                                                                                                                             tttGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCGGTGAATGGAACAACGTTCGCGGTATGCTGGGCTTCACCGGGACTAACAAAACCCGCAAAATTTCCGTAATGGGtg                                      >  1:516832/1‑238 (MQ=255)
                                                                                                                                                                                             tttGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGt                                       <  2:574530/238‑1 (MQ=255)
                                                                                                                                                                                             tttGGGCGGAGCGTTGACCCCGCCTTTGGTAGGTCACAAAAAGGATAAAACAATGGCTACCCCGCACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGGGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGt                                       <  1:606923/238‑1 (MQ=255)
                                                                                                                                                                                             tttGGGGGGGCGTTGACCCCGCCTTTGTTATGTCAAAAAAAGGATAAAACAATGGCTCCCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTc                                      <  1:101969/238‑1 (MQ=255)
                                                                                                                                                                                              ttGGGCGGAGCTTTGACTCCGCGTTTGTTATGTCACAAAAAGTATAAACCAATGGCTACCCCACCCATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCa                                     <  1:955661/239‑1 (MQ=255)
                                                                                                                                                                                                   cGGAGCGTTTACTCCGGCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACTCATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGt                                 <  2:522493/238‑1 (MQ=255)
                                                                                                                                                                                                       aGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATgg                             >  1:297164/1‑239 (MQ=255)
                                                                                                                                                                                                       aGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGACGACCCACTTAGTGCGAAGAATATTGCTGCAACTTTCCTTGAAGCTGCCCGGGAAATAAACAACGTTCGCGGCATTCTGGTCTTTAACGGTACTACAAAAGACCACAAAAATTAGCTAAAGGTCCACGGTagtg                             >  1:76872/1‑236 (MQ=255)
                                                                                                                                                                                                       aGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGACAAAAc                                                                                                                                                                                                                                    >  1:87211/1‑42 (MQ=255)
                                                                                                                                                                                                       aGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGACAAAAc                                                                                                                                                                                                                                    <  2:87211/42‑1 (MQ=255)
                                                                                                                                                                                                            ttccTCCGCTTTTGTTATGTCACAAAAAGGATAAAACAAGGGCTTCCCCACAAATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATc                        <  1:963358/236‑1 (MQ=255)
                                                                                                                                                                                                            tgtccccctcTTTTTTATTACACAAAAAGGTAAAAACAAGGGCTACCTCCCACATTAATGCAGAAATGGGCGATTTCGATGACGTAGGTTTGACGCCAGGCGACCCGCTTCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGTCCGCACAATTTCCGTAATGGGTCACGGTATGGGTATc                        <  1:904644/230‑1 (MQ=255)
                                                                                                                                                                                                            tggcTCCGCCTTTGTTATGGCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGGAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATc                        <  2:13978/236‑1 (MQ=255)
                                                                                                                                                                                                                 ccGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCCCCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGt                    >  2:252778/1‑238 (MQ=255)
                                                                                                                                                                                                                   gCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgtgcg                                                                                                                                                   >  2:718146/1‑109 (MQ=255)
                                                                                                                                                                                                                   gCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACt                                                                                                                                 >  1:160989/1‑127 (MQ=255)
                                                                                                                                                                                                                   gCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAATAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGATCACGGGATGGGTATCCCGt                    >  2:401437/1‑236 (MQ=255)
                                                                                                                                                                                                                   gCCTTTGTTATGTCACAAAAAGGACAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCg                                                                                                                                                            <  2:564084/100‑1 (MQ=255)
                                                                                                                                                                                                                     cTTTGTTATGGCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGc               <  2:616144/239‑1 (MQ=255)
                                                                                                                                                                                                                        tGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGCCCTGCTcc            >  2:607444/1‑239 (MQ=255)
                                                                                                                                                                                                                          ttATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa           >  2:757261/1‑238 (MQ=255)
                                                                                                                                                                                                                             tgtacaaaaaaGGATAAAACAATGGCTACCCCACACAGTAATGCAGAAATGGGCGATTTCGCTGATGGAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTa       <  2:986958/233‑1 (MQ=255)
                                                                                                                                                                                                                              gttactaaaaTGATTAATCAATGGCTCCCCCACACATGAATGGATGAATGGGCGATTGGGCGGGAGTAGTTTTGATGCCAGGTGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCGTGAAGATGCCCGTGAAGGGAACAACGTTCGGGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTa       <  1:49395/232‑1 (MQ=255)
                                                                                                                                                                                                                               tCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGt                                       <  2:260475/205‑1 (MQ=255)
                                                                                                                                                                                                                               tCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGt                                       >  1:260475/1‑205 (MQ=255)
                                                                                                                                                                                                                                  caAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGGTCCATCTccacca  >  2:175499/1‑239 (MQ=255)
                                                                                                                                                                                                                                  caAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATACCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCa  <  2:707352/239‑1 (MQ=255)
                                                                                                                                                                                                                                         gATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATc                        <  2:46296/210‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
ACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCA  >  NC_000913/4620643‑4621106

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