Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,678 G→A P18S (CCA→TCA)  dnaA ← chromosomal replication initiator protein DnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,6780GA98.9% 312.1 / NA 89P18S (CCA→TCA) dnaAchromosomal replication initiator protein DnaA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/47);  minor base C (1/0);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCA  >  NC_000913/3883432‑3883924
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gagCAGGCGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAAc                                                                                                                                                                                                                                                    >  2:367145/3‑251 (MQ=255)
   cAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCa                                                                                                                                                                                                                                                 <  1:381128/251‑1 (MQ=255)
    aGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCCGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCAt                                                                                                                                                                                                                                                >  2:103417/1‑251 (MQ=255)
           gcgACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCaa                                                                                                                                                                                                                                         <  1:246015/251‑1 (MQ=255)
               cGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCgg                                                                                                                                                                                                                                     <  2:569387/251‑1 (MQ=255)
                gTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCggg                                                                                                                                                                                                                                    >  2:38305/1‑251 (MQ=255)
                 ttGCTCTTCACTGCCGCTTGTTTCGTTTGCGTCTCCGGTTTGGTGCCGTCTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCATTCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTTTTATCGCTCAGTTCCGCCTTCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCggg                                                                                                                                                                                                                                    <  1:512663/250‑1 (MQ=255)
                         cACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACAct                                                                                                                                                                                                                           >  1:23123/1‑251 (MQ=255)
                            tGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACActgct                                                                                                                                                                                                                        >  1:117231/1‑251 (MQ=255)
                                    gtgGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGc                                                                                                                                                                                                                >  1:212588/1‑251 (MQ=255)
                                        cGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGaa                                                                                                                                                                                                             >  1:280909/1‑250 (MQ=255)
                                             gCGTCACCGGTTTGGTGCCGACATCAAAACGCAGCTGTGGGGCATCCGCTCCTCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTCCCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGa                                                                                                                                                                                                        <  1:86058/250‑1 (MQ=255)
                                              cGTCACCGGTTTGGTGCCGACTTCAAAACGCAACTGTGGGGCATCCGCTCCGCAGAAACTGGTGAGCAGTCCATTGATATTATTAAGGTACTTGGCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGaca                                                                                                                                                                                                      <  2:497161/251‑1 (MQ=255)
                                                 cctcgtgttgttGCCGACTTCAAAACGCAGCTGTTGGTCATCCGCTCCGCAGAAACTGGTTAGCGTACCATTTATATTATTAAGGTACGTGTCCCGTACCCAACCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGCAAGTGACAcgt                                                                                                                                                                                                   <  2:19966/240‑2 (MQ=255)
                                                    cGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACCcggcgg                                                                                                                                                                                                >  2:364300/1‑251 (MQ=255)
                                                        ttGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACt                                                                                                                                                                                             <  2:308941/250‑1 (MQ=255)
                                                          ggTGCCGACTTCAAAACGCGGCTTTTGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACGTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTTGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGTCGGACTcc                                                                                                                                                                                           <  2:352800/250‑1 (MQ=255)
                                                          ggTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTcc                                                                                                                                                                                           >  1:567904/1‑250 (MQ=255)
                                                             gcctaCTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGGCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTc                                                                                                                                                                                       <  2:59147/247‑1 (MQ=255)
                                                             gCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGCATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTc                                                                                                                                                                                       >  1:537707/1‑251 (MQ=255)
                                                                       aaaCGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGt                                                                                                                                                                             <  2:117088/251‑1 (MQ=255)
                                                                           gCAGCTGTGGGGCATCCTCTCCGCAGAAACGGGTTAGCATTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCGGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGAt                                                                                                                                                                         <  2:404489/251‑1 (MQ=255)
                                                                            cAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATa                                                                                                                                                                        <  1:263959/251‑1 (MQ=255)
                                                                                tgtgGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATg                                                                                                                                                                     >  1:157697/1‑250 (MQ=255)
                                                                                   ggggCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGc                                                                                                                                                                                                                >  1:388023/1‑204 (MQ=255)
                                                                                   ggggCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGc                                                                                                                                                                                                                <  2:388023/204‑1 (MQ=255)
                                                                                       cATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGc                                                                                                                                                             <  2:541467/251‑1 (MQ=255)
                                                                                        atccgatccgAAGAAACTGGTTAGCAGTCCATTGATATTAGTAAGGTATGTGTCCCGTACCCAATCGAGGACAAAGCGGTTTGGCGCTGACAGGGCCAGCGTGGTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATGCTGTGGCTGATAGCTCTTCCTGCAAGCGGGCAAGACACTGCTGCCAGAGCGAAAGTGACTCGGCGGGCTCCGCTCGAACGAAAGTCGAATATGATTAAGGCt                                                                                                                                                            <  2:380882/251‑1 (MQ=255)
                                                                                          ccgctccgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTTCTTTGCCCGTACCCAATCGAGGACAAAACGTTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTg                                                                                                                                                           <  1:399972/250‑1 (MQ=255)
                                                                                          ccgctccgaAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGa                                                                                                                                                          <  2:177244/251‑1 (MQ=255)
                                                                                                cgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTTCTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACAtt                                                                                                                                                    <  2:4291/251‑1 (MQ=255)
                                                                                                       aCTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAtt                                                                                                                                             <  2:157697/251‑1 (MQ=255)
                                                                                                       aCTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAtt                                                                                                                                             <  1:509953/251‑1 (MQ=255)
                                                                                                         tGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTc                                                                                                                                                                            <  2:388785/218‑1 (MQ=255)
                                                                                                         tGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTc                                                                                                                                                                            >  1:388785/1‑218 (MQ=255)
                                                                                                          ggTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCTTACAGGGCCAGCGGGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCTACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgtt                                                                                                                                          <  2:181304/251‑1 (MQ=255)
                                                                                                          ggTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGAGTGTTATCGCTCAGCTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGCTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGACCAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgtt                                                                                                                                          >  1:557873/1‑251 (MQ=255)
                                                                                                              aGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACGCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGac                                                                                                                                       <  2:353530/250‑1 (MQ=255)
                                                                                                                cAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGacgt                                                                                                                                     >  1:504611/1‑250 (MQ=255)
                                                                                                                 aGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGTATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCCCTCGCACAAAAGTCGATAATGACTAAGGCTGAAACCTTCATGATTGTTGacgtac                                                                                                                                   >  2:29342/1‑251 (MQ=255)
                                                                                                                          gATATTATTAAGGTACTTGTCCCGTCCCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTTCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGaaa                                                                                                                            <  2:280909/250‑1 (MQ=255)
                                                                                                                             attattACGGTACTTGTCCCGTACCCAATCGAGGTCAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGACGTATCCACATCCTGAATTCTGCGGCTGATAACTACTCCTGCAACCGGGCAAGACACGGCTGCCACAGCGAAAGTGACCCGGCGGACTCCACTCGACCAACAGTCGATAATGACGAAGGCTGAAAGATTCATGATTGTTGACGTACCGCGAAAAGAc                                                                                                                        >  2:4039/1‑251 (MQ=255)
                                                                                                                             attattAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAc                                                                                                                        >  1:56251/1‑251 (MQ=255)
                                                                                                                             attattAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCACTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGa                                                                                                                         <  1:347213/250‑1 (MQ=255)
                                                                                                                               tattaAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTCCGTCGAAAAGAcac                                                                                                                      >  2:49040/1‑249 (MQ=255)
                                                                                                                               tattaAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAccc                                                                                                                      >  2:133416/1‑251 (MQ=255)
                                                                                                                                  taAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTg                                                                                                                    <  2:193435/250‑1 (MQ=255)
                                                                                                                                           gagcccGTCCCCAATCGAGTACAACAGGGTTTCGCGCGTACAGGGCGAGCGTTGTATCGGTCAGTTACGCCTGCAATGGGCGTATCCACAAACTGAATTCTGTTGCTAATAACTCATCCGGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGACGGACGCCACTCGAACAAAAGGCGAAAATGACTAAGGCTGAAACAGTCAGGATTATTGACGGACGTCGAAAAGGCCCTGCATGAGGGTg                                                                                                          <  2:296629/248‑1 (MQ=255)
                                                                                                                                              cccGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACg                                                                                                       >  1:460069/1‑251 (MQ=255)
                                                                                                                                              cccGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACg                                                                                                       >  1:56795/1‑251 (MQ=255)
                                                                                                                                                  tACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACg                                                                                                   >  1:528516/1‑251 (MQ=255)
                                                                                                                                                    cccAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACgt                                                                                                  >  2:138506/1‑249 (MQ=255)
                                                                                                                                                    cccAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGaa                                                                                                 <  1:43202/251‑1 (MQ=255)
                                                                                                                                                      cAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGGTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAAcc                                                                                               <  1:544789/251‑1 (MQ=255)
                                                                                                                                                           gAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGt                                                                                          <  1:103417/251‑1 (MQ=255)
                                                                                                                                                               aCAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCCCGAACCGCTGTCTGc                                                                                      >  2:426271/1‑251 (MQ=255)
                                                                                                                                                                 aaaaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAgg                                                                                                                                                              <  1:443274/177‑1 (MQ=255)
                                                                                                                                                                 aaaaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAc                                                                                                                        >  2:14343/1‑215 (MQ=255)
                                                                                                                                                                 aaaaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAc                                                                                                                        <  1:14343/215‑1 (MQ=255)
                                                                                                                                                                 aaaaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTAGAACAAAAGTCGATAATGACTAAgg                                                                                                                                                              >  2:443274/1‑177 (MQ=255)
                                                                                                                                                                   aaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGtt                                                                                  <  1:38305/251‑1 (MQ=255)
                                                                                                                                                                    aCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGGGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGAATAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTta                                                                                 <  2:23123/251‑1 (MQ=255)
                                                                                                                                                                    aCGGATTGGAGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATAATGCATTCTGTGGCTGATAACTCATCCTGCACTCGGGCAAGAAACAGCTGCCAAAGCGAACGTGCCACGGAGGACTCCACTCGACCAAAAGTCGATAATGACTAAGGCTGAAACCTTCATGATAGTTGACGAACGTCGAAAAGAACCTGCTTGAGGGTGACGCACGGACCGCTTACTGCGGttt                                                                                 >  2:368791/1‑250 (MQ=255)
                                                                                                                                                                     cGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGCACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTtat                                                                                >  2:496638/1‑251 (MQ=255)
                                                                                                                                                                      ggTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTtat                                                                                >  1:346885/1‑250 (MQ=255)
                                                                                                                                                                        tttGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTatatt                                                                             >  2:180671/1‑250 (MQ=255)
                                                                                                                                                                           ggCGCGTACTGGGCCAGCGTGTTATCGCTCAGTTCAGCCTGCGCTGGGCGTATCCACATACTGTATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAACGGCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTGGGGGGTGACGCACGAAACGCTGACTGCGGTGATATGccc                                                                          >  2:73183/1‑251 (MQ=255)
                                                                                                                                                                               cgTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATc                                                                      >  2:310765/1‑251 (MQ=255)
                                                                                                                                                                                  aCAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAg                                                                   <  2:56795/251‑1 (MQ=255)
                                                                                                                                                                                   cAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                  <  1:313906/251‑1 (MQ=255)
                                                                                                                                                                                     gggCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                <  2:59478/251‑1 (MQ=255)
                                                                                                                                                                                            cGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGACCGATCAAGATCCTGCaa                                                          >  2:487755/1‑250 (MQ=255)
                                                                                                                                                                                              tgtTATCGCTCAGTTCCGCCTGCAATGGGCGTATCTACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGAACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAAACGCTGTCTGCGGTTATATGCCCGATCAGGATCCTGCaaacc                                                       >  2:145251/1‑249 (MQ=255)
                                                                                                                                                                                              tgtTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGAGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       >  1:230682/1‑251 (MQ=255)
                                                                                                                                                                                                        cAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCACCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGAc                                              <  1:280109/250‑1 (MQ=255)
                                                                                                                                                                                                              cGCCTGCAATGGGCGTAGCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGAt                                       <  2:346885/251‑1 (MQ=255)
                                                                                                                                                                                                                           cGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACt                          >  1:526786/1‑251 (MQ=255)
                                                                                                                                                                                                                             tATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGCGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgc                        >  1:378464/1‑251 (MQ=255)
                                                                                                                                                                                                                                   caTACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAaga                  <  2:56251/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aTACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAagag                 <  1:57792/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tACTGAATTCTGTGGCTGATAACTCATCCTGCAATAGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGGACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAACGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGCACAAGATCATGCAAAACGATCGGGACCGCGGATCAAAGCCTAAACTGCGCAagaa                 >  2:578304/1‑249 (MQ=255)
                                                                                                                                                                                                                                      aCTGAATTCTGTGGCTGATAACGCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCGGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAGCGGGACCGCGGATCATAGCCTAAACTGCGCAAgaga                <  1:40281/250‑1 (MQ=255)
                                                                                                                                                                                                                                          aaTTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc           <  1:398234/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                           >  2:514146/1‑231 (MQ=255)
                                                                                                                                                                                                                                              cTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                           <  1:514146/231‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtgGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGACCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGttt      >  1:551402/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tgtgGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGtt       <  2:130702/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gtgGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtc     >  1:428288/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   gCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTca  <  2:567904/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GTGCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCA  >  NC_000913/3883432‑3883924

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