Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,996 T→G D70A (GAT→GCT)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,9960TG99.1% 393.4 / NA 109D70A (GAT→GCT) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base G (61/46);  minor base A (1/0);  total (62/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 109 total aligned reads displayed.

ATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCT  >  NC_000913/4193749‑4194240
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
aTGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGCTCGCCAGCCGCCAGTAAGCGt                                                                                                                                                                                                                                                        >  1:358223/1‑251 (MQ=255)
  gcctgcccaggttttCCAGCCCCTGCGGGGCGGATGGCAGCTGTTGGGTGTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCGGGGAGACGTCTATTTCCATACCGTCATGTGTCGAACCGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTCCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTGGATCGCAAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                       <  2:464377/237‑1 (MQ=255)
    cGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGAGCGGATGGAATTGCGGTCGGTGGCTTGGAAATCTTCCTGCCCC‑AAAAGGGCGCCATACGCCTGCTGCTTGGTAGCCAGCCGCCCGCAAGCGTTGAt                                                                                                                                                                                                                                                    >  2:309118/1‑251 (MQ=255)
            cTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAAAAATc                                                                                                                                                                                                                                            >  2:529265/1‑251 (MQ=255)
                         tGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcg                                                                                                                                                                                                                               >  1:295616/1‑251 (MQ=255)
                          gCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcg                                                                                                                                                                                                                               <  1:460964/250‑1 (MQ=255)
                                    aggcttttGTTTGGTTTGGCTCGGGACCACATGTCCCCGCGCGATATAAGGGGGGCGTGCTCCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAGCGCCCAGCCGCAGGCCGG‑CCGCGCGGATGGCATGGAGATCGGTGTCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATCCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAgcg                                                                                                                                                                                                                   <  2:397964/244‑1 (MQ=255)
                                       cATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGCTAAACATTCGCGCGTTATAGCGG‑CGGACCCGCGCa                                                                                                                                                                                                                 >  2:781173/1‑250 (MQ=255)
                                        aTCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAAt                                                                                                                                                                                                               >  1:173425/1‑251 (MQ=255)
                                        aTCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGAGTCAAAACGCCCAGCCGCAAGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAgcgacat                                                                                                                                                                                                               >  1:452133/1‑247 (MQ=255)
                                          cTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTg                                                                                                                                                                                                             >  1:623370/1‑251 (MQ=255)
                                           ggtttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTgc                                                                                                                                                                                                            <  1:767854/250‑1 (MQ=255)
                                            gtttggtttgCGTCGGGAACCCATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTgcc                                                                                                                                                                                                           <  1:492769/251‑1 (MQ=255)
                                            gttgggtttgCGTCGGGAACACCTGTCCCAGATCGATAGAAGAGGGGCGTGCGGCCAGCGCGACGTCGATTTCCATAGCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGAGGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTgcc                                                                                                                                                                                                           <  2:500441/247‑1 (MQ=255)
                                                     gcttCGGGAACACTTGTCCCAGCGCGATATAAGAGGGGCGTGCTGCAAGCGCGACGTCAATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGAc                                                                                                                                                                                                   <  1:296303/247‑1 (MQ=255)
                                                       gTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATc                                                                                                                                                                                                <  1:539137/251‑1 (MQ=255)
                                                         cGGGACCACATGTCCCAGCGCGATATAAGAGGGGCGTGCTCCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAGCGCCCAGCCGCGGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATCCGCCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGT‑CGGCCCAGCGCAATTGCCGCCACGACATCgg                                                                                                                                                                                              <  1:186608/251‑1 (MQ=255)
                                                          gggAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGc                                                                                                                                                                                             >  1:529053/1‑251 (MQ=255)
                                                          gggAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGc                                                                                                                                                                                             >  1:639284/1‑251 (MQ=255)
                                                                   tGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGGGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCCGCGCCAATGCCCCCACACCCACGGC‑TTCCAcccc                                                                                                                                                                                   >  2:493389/1‑249 (MQ=255)
                                                                    gccccAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACctc                                                                                                                                                                                   <  1:636724/248‑1 (MQ=255)
                                                                                 ataAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGAGCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGAt                                                                                                                                                                     <  1:281565/251‑1 (MQ=255)
                                                                                   aagcggggCGTGCTGAAACCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGCTCGGTGGCTTGCAAATCTTCCTGCCCC‑ACATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGACCCATTGCCGACCCGCCCTAGCC‑TTCCAACTCTTCATCGCGCCgtcct                                                                                                                                                                   >  2:495094/5‑247 (MQ=255)
                                                                                   aaGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACCACATCAGC‑TTCCCCCCCTTCCTCGCCCCGGTCt                                                                                                                                                                   >  2:543707/1‑251 (MQ=255)
                                                                                      agGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCttt                                                                                                                                                                 >  1:702080/1‑250 (MQ=255)
                                                                                              gctgcCAGCGGGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTa                                                                                                                                                         <  1:691049/250‑1 (MQ=255)
                                                                                              gctgcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGt                                                                                                                                                          >  1:70864/1‑249 (MQ=255)
                                                                                                 gcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGCCATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGCTGCGTAGCTg                                                                                                                                                     >  1:160628/1‑251 (MQ=255)
                                                                                                   cAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTTCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGAAGCTggg                                                                                                                                                   >  1:759267/1‑250 (MQ=255)
                                                                                                    aGCGCGACGACGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCACCAGCCGCCAGTCAGCGTTGATAAACAATCTCGCGTTATAACGGACGACCC‑GCGGACGTGCCGCCACGCCCTTGGG‑TACCACCTCCTCATCGGCCCGATCtttagtgccaatccag                                                                                                                                                    >  2:210908/1‑236 (MQ=255)
                                                                                                       gcgACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAgt                                                                                                                                              <  1:84223/251‑1 (MQ=255)
                                                                                                         gACGTCGATTCCCATATCGTCATTTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAgtac                                                                                                                                            <  1:115868/251‑1 (MQ=255)
                                                                                                          gcGTCTATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAgtacg                                                                                                                                           <  1:404516/250‑1 (MQ=255)
                                                                                                              tgATTTCCATAACGTGATGTGTCGAACCGCCCAGCCGCAGGCACG‑CCGCGCGGAGGGCATTGAGACCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAgtacgtac                                                                                                                                        <  1:359051/249‑1 (MQ=255)
                                                                                                                       tatCGTCATGTGACGAAACGCCCCGCCGCAGGCCTG‑CCGCGCGGATGACATTGATATCGGTGGCTTGCAAATCTTCCTGCCCC‑AGATGGACACCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAACGTTGATAAACCATCGCGCGTTATCGCGG‑CGGCACAGCGCAATTGACGCCACGACCACGACGTTCCAC‑TCGTCACCGACCCGATCTTTGATGAGTCGCAGG‑GGAGTACGGAAGCCCGCATcc                                                                                                                              >  2:61141/1‑251 (MQ=255)
                                                                                                                        atCGTCATGTGCCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATcc                                                                                                                              <  1:574311/250‑1 (MQ=255)
                                                                                                                                gggcggaaccccccAGCCGCAG‑CCTGTCCGGGCGGAGGGCATGGAGATGGGTGGCTTGAAACTCTCCCGCCCCC‑AAATGGACGCCATACGCCTGAGGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACg                                                                                                                      <  2:554185/239‑1 (MQ=255)
                                                                                                                                  gtCGAACCGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGGGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTc                                                                                                                   <  1:519383/251‑1 (MQ=255)
                                                                                                                                    cGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCg                                                                                                                  <  1:529265/250‑1 (MQ=255)
                                                                                                                                          ccccAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑t                                                                                                          <  1:297581/250‑1 (MQ=255)
                                                                                                                                            ccAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCCGCGCAATTGCCGCCACGACCTCGGC‑TTCCACCTCTTCATCGCTCCGATCTTTGATGCGTTGCTGG‑AAAGTCCGAACGCCTTCATCCAACAGACGATCGATCCaa                                                                                                            >  2:671535/1‑247 (MQ=255)
                                                                                                                                              gcccccAGGCCTG‑CCCCCCGGATGCCATTGAGAGCGGTGGCTGGCAA‑TCCTCCCGCCCCCAAATGGACGCCTAACGCCTGATTCTGGTCCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTAc                                                                                                      <  2:174132/246‑1 (MQ=255)
                                                                                                                                              aGCCGCAGGCATG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑AGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGAGCTCTGATGCGTAGATGG‑AGAGAACGGACCCCTGCATCCAACAGACGTTCGATCCAC‑TGTa                                                                                                       >  1:492817/1‑250 (MQ=255)
                                                                                                                                                  gCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTAAGCCTGCATCCAACAGACGTTCGATCCAC‑TGTCCGCTg                                                                                                  >  1:550397/1‑251 (MQ=255)
                                                                                                                                                   cAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGAACTTTATTTCGTAGCTTG‑AGAGTACGGCCCCCCGCATCCCACAGACGTTCGATCCAC‑Tggaccctgg                                                                                                 >  2:68094/1‑243 (MQ=255)
                                                                                                                                                    aGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATCCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGGCATCGGC‑CTCCCCCTCTTCCTCCGGCCGCTCTTTGCTGCGTAGCTGG‑GGTGTACGTACGCCTGCATCCACCAGCCTGCCGATCCCC‑TGTACGCCGT‑c                                                                                               >  2:114797/1‑251 (MQ=255)
                                                                                                                                                    aGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCATCAGACGTTCGATCCAC‑TGTACGCTGT‑c                                                                                               >  2:540295/1‑251 (MQ=255)
                                                                                                                                                    aGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑c                                                                                               <  1:19213/251‑1 (MQ=255)
                                                                                                                                                     ggCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑cc                                                                                              >  1:816592/1‑251 (MQ=255)
                                                                                                                                                         tG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCCTCGCGCCGATCTTTGATGCGTAGCTGG‑GGAGGACGTACGCCTGCATCCAAAAGACGTTCGGTCCAC‑TGTAAGCTGT‑ccaccc                                                                                          >  2:228712/1‑250 (MQ=255)
                                                                                                                                                         tG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGGAGCTGGCAG‑GTACGTACGCCTGCCTCCAACAGACGTTCGATCCAC‑TGGACGCtagccacc                                                                                            >  2:641937/1‑243 (MQ=255)
                                                                                                                                                                 cgGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGt                                                                                    >  1:154611/1‑250 (MQ=255)
                                                                                                                                                                  ggATGGCATTGAGATCGGTGGCTTGCAACTCTTCCTGCCCC‑AACTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCACCGAAATAGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGGTCTGTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCAACCACCAAACGTTCGATACCC‑TGTAAGCAGA‑CCACCACCCGGTAc                                                                                  >  2:307078/1‑251 (MQ=255)
                                                                                                                                                                  ggATGGCATTGAGATCGGTGGCTTGCAAATCTTACTGACCA‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCCGTAAGCGTTGATAAACACTCGCGCGTTATAGCGG‑CGGCCCAGCGCCATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGCTCTTTGATGCTTAGCTGG‑AGAGTACGTACGCCCGCATCCCACAGACGTTCGATCCCC‑TGTACGCTGTACCAC‑ACCGGGTAc                                                                                  >  2:802250/1‑251 (MQ=255)
                                                                                                                                                                    aTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCCTCGCGCCGGTCTTTGCTGCGTAGCTGG‑AGAGTACGTACGCCTGCCTCCAAAAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACac                                                                                >  2:354426/1‑250 (MQ=255)
                                                                                                                                                                      ggCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGAAGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGAAATTGCCCCCACGACATCGGC‑TTCCACCTCTTCATAGCGCCGATCTTTGATGCGTAGCTGG‑AGGGTACGTACCACTGAATCCCCAAGAAGTTAGTTCCCA‑TGTACGCTGT‑CCACCACCGGGtcaagtc                                                                              >  2:482073/1‑245 (MQ=255)
                                                                                                                                                                        cATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCt                                                                            >  2:253550/1‑251 (MQ=255)
                                                                                                                                                                           tGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGtt                                                                                                                                                                                                                                     >  1:722525/1‑100 (MQ=255)
                                                                                                                                                                           tGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGa                                                                          >  1:443687/1‑250 (MQ=255)
                                                                                                                                                                           tGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGtt                                                                                                                                                                                                                                     <  2:722525/100‑1 (MQ=255)
                                                                                                                                                                            gatatggggggCTGGCAAACCTCCCCGCCCC‑AAATGGACGCCATACGCCTGTGGCTTTACCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAAc                                                                        <  2:39282/242‑1 (MQ=255)
                                                                                                                                                                            gagaTCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGCTCCAC‑TGTACGCTGT‑CCACACCCGGGAACCGTCCTGAAc                                                                        >  2:321428/1‑251 (MQ=255)
                                                                                                                                                                              aaTCGATGGCTTGCAACTCATCCTGCCCC‑AACTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATCACCACTCGCGCGTTATAACGG‑CGGCCCACCGCAATTGCCGCCACGCCCTCGGC‑TTCCACCAATTCATCACGCCGATCTTTTATGCGTATCTTG‑AGAGTACGTACGCCTGCATCCAACAGCCGTTCAATC‑ACATGTACGACGT‑CCACCACCGGGTACAGGCCTGAAcc                                                                       >  2:788915/2‑249 (MQ=255)
                                                                                                                                                                                ccGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGaaa                                                                    <  1:253550/250‑1 (MQ=255)
                                                                                                                                                                                 cGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGaaa                                                                    >  1:13180/1‑250 (MQ=255)
                                                                                                                                                                                 cGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGACGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTTCCGGCCCGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGTTCCAA‑TGTACGCTGT‑CCACCACCGGGtccgatcctgcaagaaa                                                                    >  2:495457/1‑234 (MQ=255)
                                                                                                                                                                                   gTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAACCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCAAGACATCCGC‑TTCCACCTCTTCCTCGGCCCGATCTTTTATGCGTAGCTGG‑GGAGTACGTTCGCCCTCATCCAACAGACGTTCGAATCAC‑TGGACCCCTT‑CCACCACCCGGGACAGTccctaaagaaa                                                                    >  2:214697/1‑239 (MQ=255)
                                                                                                                                                                                     ggCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCt                                                                                                                                                                   <  1:318585/154‑1 (MQ=255)
                                                                                                                                                                                     ggCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCt                                                                                                                                                                   >  2:318585/1‑154 (MQ=255)
                                                                                                                                                                                      gcgtgaaaatttcccTGCCCA‑AATTGGCCGCCATCCGCTGGATGTTTGGCCGCCAGCCGACAGTAAGCGTGGACAAACAATCGCGCGTTATAGCGG‑CGGCCAAGTGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTGGATGCGTAGCTGG‑AGAGTACGTACGCCTGCAGCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCCCCGGGTACAGGCCTGAACGAAAAGGTAc                                                              <  2:358615/238‑1 (MQ=255)
                                                                                                                                                                                         tGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCAAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGAGTACAGTCCTGAACGAAAAGGTACagg                                                           >  2:555728/1‑251 (MQ=255)
                                                                                                                                                                                          gCAAATCTCCCTCCCCC‑AAAGGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACagg                                                           <  1:627957/250‑1 (MQ=255)
                                                                                                                                                                                               ccTTCCTGCCCC‑AAAGGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGaaa                                                     <  1:556140/250‑1 (MQ=255)
                                                                                                                                                                                                      gCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGc                                              >  1:40539/1‑251 (MQ=255)
                                                                                                                                                                                                          c‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGAt                                          >  1:493952/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATAc                                        >  1:362830/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTGGACGCCATACGCCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATAc                                        <  1:811171/251‑1 (MQ=255)
                                                                                                                                                                                                                 gACGCCATACGCCGGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCa                                    <  1:321428/251‑1 (MQ=255)
                                                                                                                                                                                                                   cGCCATACGCCTGATGTTTGATCTCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACg                                  <  1:404711/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCTCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCAGCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCt                                <  1:257784/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAACGTAAAGGAGGAAAATCAGGCTgctaactcccgct                                >  2:795786/1‑239 (MQ=255)
                                                                                                                                                                                                                      cgtgcgcctaatgcctggccGCCATCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCTCCTCTTCAGCGCGCCTATCGTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCtt                               <  1:364387/232‑1 (MQ=255)
                                                                                                                                                                                                                      cATACGCCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCtt                               <  1:563187/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCa                                                                                                             >  1:122552/1‑174 (MQ=255)
                                                                                                                                                                                                                       aTACGCCGGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCa                                                                                                             <  2:122552/174‑1 (MQ=255)
                                                                                                                                                                                                                            ccTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCGGCCTCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCTTGAACGACAAGAGACAGGAGGAAAAGCAGGCTGAGAACTCACGcctcctcct                         >  2:227653/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGGTCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCtt                        >  1:604836/1‑250 (MQ=255)
                                                                                                                                                                                                                               gATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAACCAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTGCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTAc                      <  1:31870/251‑1 (MQ=255)
                                                                                                                                                                                                                                atctttgtgcgcaccccgccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCAGGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCTTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACAGCACGCTTCCTCCTTAc                      <  2:816592/235‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATCCCCCGCCGCCAGTAAGCGATGATAAACAATCGCGCGTTAGAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTGCATCGCGCCGATCTTGGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAg                  <  1:309118/249‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAtt           >  1:165086/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGAGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGGACAGTCATGAACGAAAAGGAAAAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAtt           >  2:670936/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAttt          >  2:299515/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTATGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTtc         >  1:390197/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ggccgccAGTAAGCGTTGATAAACAACCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCTCGCCATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTtct        <  2:529053/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCCTCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTtct        >  2:379345/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gccgccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTtctc       >  1:228280/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cgccAGTAAGCGTGGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCg     <  1:141202/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACCTCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGcc   >  1:331245/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGG‑CGTCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCct  >  1:300768/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
ATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTG‑CCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC‑AAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGG‑CGGCCCAGCGCAATTGCCGCCACGACATCGGC‑TTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGG‑AGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAC‑TGTACGCTGT‑CCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCT  >  NC_000913/4193749‑4194240

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: ‑