Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,488 G→A L61L (CTG→CTA mog → molybdopterin adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,4880GA98.7% 237.7 / NA 77L61L (CTG→CTAmogmolybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (53/23);  total (54/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAG  >  NC_000913/9243‑9726
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
tctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtagtggg                                                                                                                                                                                                                                           >  1:142216/1‑250 (MQ=255)
         cccTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGa                                                                                                                                                                                                                                   <  2:479322/250‑1 (MQ=255)
         cccTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGAGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGCAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATg                                                                                                                                                                                                                                    >  1:38293/1‑249 (MQ=255)
             tCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTg                                                                                                                                                                                                                               <  2:466669/250‑1 (MQ=255)
              cGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGcc                                                                                                                                                                                                                             >  1:328550/1‑251 (MQ=255)
                cAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCAt                                                                                                                                                                                                                           >  1:67089/1‑251 (MQ=255)
                 aGTGAATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTCGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGAAAACCCGCTTAATCCCCGATGAGCAGGCGCTCATCGAGCACACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATc                                                                                                                                                                                                                          >  2:161841/1‑251 (MQ=255)
                       tCATTCTGTTTAACGAGCCTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGc                                                                                                                                                                                                                    >  2:563832/1‑251 (MQ=255)
                            cTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTcacc                                                                                                                                                                                                               >  2:503500/1‑251 (MQ=255)
                             tGTTTAACGAGACTGTTTAAACGGAAACATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGACTGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCAGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAACTGAGTTGCCATCTGGTGCTcacc                                                                                                                                                                                                               >  2:480779/1‑250 (MQ=255)
                              gTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTcaccac                                                                                                                                                                                                             >  2:140029/1‑251 (MQ=255)
                              gTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCGcaccac                                                                                                                                                                                                             <  2:272288/251‑1 (MQ=255)
                                ttAACGAGCCTGATTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACgg                                                                                                                                                                                                           >  2:433443/1‑251 (MQ=255)
                                 tAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGCACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACggg                                                                                                                                                                                                          >  1:398220/1‑251 (MQ=255)
                                 tAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACggg                                                                                                                                                                                                          >  1:267396/1‑251 (MQ=255)
                                      agaCTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATg                                                                                                                                                                                                                                         <  2:201158/215‑1 (MQ=255)
                                      agaCTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATg                                                                                                                                                                                                                                         >  1:201158/1‑215 (MQ=255)
                                           gTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGc                                                                                                                                                                                                >  1:392146/1‑251 (MQ=255)
                                               aaaCGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGCGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCgg                                                                                                                                                                                            >  2:147045/1‑251 (MQ=255)
                                                  cGGAAAAATCTTGATGAATACTTTCCGTATTGGCTTAGTTTCCATCTCTGATCGAGCATCCAGCGGCGTTTATCAGGATAAAGGCATACCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGgcgc                                                                                                                                                                                         >  2:396140/1‑251 (MQ=255)
                                                  cGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATACCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGCTCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCAGgcgc                                                                                                                                                                                         >  2:212188/1‑251 (MQ=255)
                                                     aaaaaTCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGcgtc                                                                                                                                                                                      >  1:1314/1‑251 (MQ=255)
                                                           cTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACg                                                                                                                                                                                >  1:281360/1‑251 (MQ=255)
                                                                        ttACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGGCGTGACGTAACGCCCGATGc                                                                                                                                                                   >  2:443246/1‑251 (MQ=255)
                                                                         tACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGc                                                                                                                                                                   <  1:530750/250‑1 (MQ=255)
                                                                          aCGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCg                                                                                                                                                                  >  1:175993/1‑250 (MQ=255)
                                                                            gTATTGGCTTAGTTTCCATCTCTGATCGCGCCTCCAGCGGCGTTTATCAGGATAAAGACATCCCTGCGCTGGAAGCATGGCTGACATCGGCGCTAACCACGCCGTTTGCACTGGAAACCCGCTTAATCCCCGATGCGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCgc                                                                                                                                                                 >  2:621669/1‑248 (MQ=255)
                                                                             tATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGGCGTGACGTAACGCCCGATGCGACGc                                                                                                                                                              >  2:466989/1‑251 (MQ=255)
                                                                                 ggCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCAGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGTTGCGACGCTGGc                                                                                                                                                          >  1:56330/1‑251 (MQ=255)
                                                                                          tCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACAACGCCGTTTGAACTGGAAACCCGATTAATCCCCGATGAGCAGGCGATCATCGAGCAAAAGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCAGGCGCGTCGTGACGTAACGCCCGATGCGCCGCTGGCAGTAGCgg                                                                                                                                                  >  1:562912/1‑250 (MQ=255)
                                                                                            cATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATACAGGCCTCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCTAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGAcc                                                                                                                                               >  2:275935/1‑251 (MQ=255)
                                                                                             ctctcTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCAGTTTGAACTCGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTAGCCATCTGGTGCTCACAACGGGCGGAACTGGCCCAGCGCGTCGAGACGTAACGCCCGATGCGCCGCTGGCAGTAGCGGACcg                                                                                                                                              >  1:315155/2‑251 (MQ=255)
                                                                                                   gATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCgaga                                                                                                                                         <  2:128953/250‑1 (MQ=255)
                                                                                                      cgcgCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGc                                                                                                                                      <  1:506127/250‑1 (MQ=255)
                                                                                                         gcATCCAGCGTCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGTGGACCGCGAGATGCCTgg                                                                                                                                  <  1:394527/251‑1 (MQ=255)
                                                                                                          cATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGc                                                                                                                                 >  1:537396/1‑251 (MQ=255)
                                                                                                           aTCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGAAGCCTGGCt                                                                                                                                <  2:469361/251‑1 (MQ=255)
                                                                                                            tCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGATGACATCGGCGCTAACCACGCCGTTTGCACTGGAAGACCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAGTAGGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCtt                                                                                                                               >  1:582104/1‑251 (MQ=255)
                                                                                                             ccAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCTCTAACCACGCCTTTTTAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTACCGCCCGATGCGACGCTGGCAGTAGCGGACCGCTAGATTCCTggttt                                                                                                                               <  2:166917/250‑4 (MQ=255)
                                                                                                              cAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTg                                                                                                                             <  1:212188/251‑1 (MQ=255)
                                                                                                                gcggcgTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGt                                                                                                                           <  1:455423/251‑1 (MQ=255)
                                                                                                                   gcgTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAAATGGAAAACCGCTTACTCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGAAGTAACGCCCGATGCGACACTGGCAGTAGCGGACCGCGCGATGCCTGGCTTTGGTGa                                                                                                                         >  1:89175/1‑250 (MQ=255)
                                                                                                                           cAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATgc                                                                                                                 <  1:621669/250‑1 (MQ=255)
                                                                                                                            aGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCACCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATgcgc                                                                                                               >  1:56757/1‑251 (MQ=255)
                                                                                                                              gATAAAGGCATCCATGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACTCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCTGACCGCGGGATGCCTGGCTTTGGTGGACAGATGCGcc                                                                                                              >  2:288173/1‑250 (MQ=255)
                                                                                                                                tAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGAGAGATGCCTGGCTTTGGTGAACAGATgcggc                                                                                                              >  1:483297/1‑246 (MQ=255)
                                                                                                                                 aaaGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCGTTGGTGAACAGATGCGCCAGAt                                                                                                          <  2:67089/251‑1 (MQ=255)
                                                                                                                                    ggCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGAta                                                                                                         >  1:552699/1‑248 (MQ=255)
                                                                                                                                         cccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGc                                                                                                  <  1:70936/251‑1 (MQ=255)
                                                                                                                                          ccTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGACTGc                                                                                                  >  1:30822/1‑250 (MQ=255)
                                                                                                                                           cTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGACGTTTGAACTGGAAACCCGCTTAATCCCAGATGAGCCGGCGATCCTCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCAAGGGCGGAACTGGCCCGGCGCGTCGTGACGTCACGCCCGATGCGACGCTGGCAGTAGAGGACCGCGCGATGCCTGGCTTTGGTGAAGAGATGCGCCAGATAAGCCTGCAt                                                                                                >  1:276595/1‑251 (MQ=255)
                                                                                                                                                     agaaTGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGCTGCGCCAGATAAGCCTGCATTTTGTACCa                                                                                       >  2:433844/1‑250 (MQ=255)
                                                                                                                                                              gACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTTGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGAt                                                                               <  1:396140/249‑1 (MQ=255)
                                                                                                                                                                cATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGCTGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCtt                                                                           >  1:312779/1‑251 (MQ=255)
                                                                                                                                                                    ggCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAACCGTTGTGTGAGCTAGAGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTCGCCCGACGCGTCGTGACGTAACGCCCGATGCGCCGCTGGCAGTAGCGGACCGCAAGATGCCTGGCTTTGGTGACCAGATGCGCAAGATAAGCCTGAATTTTGTACCAACTGCGATCCTTTcgc                                                                       >  1:618164/1‑251 (MQ=255)
                                                                                                                                                                            ccACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGAAAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGtt                                                               >  2:542528/1‑250 (MQ=255)
                                                                                                                                                                                      tGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAgggtggtgtgatt                                                      <  1:466989/250‑12 (MQ=255)
                                                                                                                                                                                       gAACTGGAAACCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCg                                                    <  2:293509/250‑1 (MQ=255)
                                                                                                                                                                                           tGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCaa                                                 >  2:339280/1‑250 (MQ=255)
                                                                                                                                                                                              aaaCCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAg                                             <  1:241964/251‑1 (MQ=255)
                                                                                                                                                                                                aCCCGATTAATCCCCGATGAGCAAGCGATCATCGAGCACAAGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTACCGCCCGGTGCGACGCTGGCAGTAGCGGACAGCGAGATGCCGGGCTTTGGTGTACCGATGCGCCAGATAAACCTGCATTTGGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAgac                                           >  2:249565/1‑249 (MQ=255)
                                                                                                                                                                                                  ccGCTTAATACCCGCTGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCAGAGCGTCGTGACGTAACGCCAGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCCTTGGTGAACAGAGGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGttgctgtgattctcacaccggcgc                                         >  2:496726/1‑228 (MQ=255)
                                                                                                                                                                                                    gCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTg                                       >  2:520273/1‑251 (MQ=255)
                                                                                                                                                                                                     cTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTg                                       >  2:10752/1‑250 (MQ=255)
                                                                                                                                                                                                       tAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGAt                                     <  1:452844/250‑1 (MQ=255)
                                                                                                                                                                                                            cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGACATGAGTTGCCATCTGGTGCTCCCCACGGGCGGAACTGGACCGGCGCGCCGTGACGTAACGCCCGATGCGCCGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGCTAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTa                                >  1:592982/1‑250 (MQ=255)
                                                                                                                                                                                                              cGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACt                             <  1:140029/251‑1 (MQ=255)
                                                                                                                                                                                                              cGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACt                             >  2:192175/1‑251 (MQ=255)
                                                                                                                                                                                                                 tGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTAc                          >  1:100363/1‑251 (MQ=255)
                                                                                                                                                                                                                      aGGCGATCATAGAGAAACCGTTGTGTGAGATAGTGGATGAAATCAGTTGCCATATGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGCCGTCACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCCGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGt                     >  2:85740/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTc                    >  2:39808/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAg                  <  2:267396/251‑1 (MQ=255)
                                                                                                                                                                                                                          gATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTACCTTACCCGGTCAGc                 >  2:26291/1‑251 (MQ=255)
                                                                                                                                                                                                                           atcatcGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGcc                <  1:339280/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTc          >  1:257045/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAt       >  1:220056/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTaa    >  1:419594/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gtgtgAGCTAGTGGATGAAATGCGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAag  >  2:257675/1‑248 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
TCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAG  >  NC_000913/9243‑9726

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