Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 137,378 C→T T99M (ACG→ATG)  cueO → multicopper oxidase CueO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913137,3780CT97.4% 271.7 / NA 80T99M (ACG→ATG) cueOmulticopper oxidase CueO
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (36/40);  minor base A (2/0);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATG  >  NC_000913/137134‑137623
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCGAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGc                                                                                                                                                                                                                                                   >  1:691126/1‑250 (MQ=255)
gCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGGGAAAATACAGCGCGGCAAAGCGGGAAAGGTTGATATCTAAAACCCACTGACGGGAGAGACAAAGTTgg                                                                                                                                                                                                                                                   >  2:141987/1‑249 (MQ=255)
        gCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACg                                                                                                                                                                                                                                         <  1:467058/251‑1 (MQ=255)
           gtgGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCGAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGc                                                                                                                                                                                                                                      <  2:691126/251‑1 (MQ=255)
             ggAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCt                                                                                                                                                                                                                                     <  1:150895/250‑1 (MQ=255)
                 ccGCGCAGTATTTGCGGCAGAACGCCCAGCGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAg                                                                                                                                                                                                                                <  1:474858/251‑1 (MQ=255)
                    cgcAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTAc                                                                                                                                                                                                                             <  2:210452/251‑1 (MQ=255)
                    cgcAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGGGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAAATGACGGGAGAGACAATGTTGCACTGGCCAGGGCTGGGAGTAc                                                                                                                                                                                                                             >  2:160851/1‑251 (MQ=255)
                          aTTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGt                                                                                                                                                                                                                        <  1:264924/250‑1 (MQ=255)
                               cGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGGGAAGTc                                                                                                                                                                                                                  >  1:246718/1‑251 (MQ=255)
                                  cAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGAc                                                                                                                                                                                                               <  1:141987/251‑1 (MQ=255)
                                    gAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTTGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAAAGGCAATCTGCTGGGGCCGGCGGGGAAAATAAAGCGCGGGAAAACGGCAACGGTTGGTATCTACAACCACATGAACGGAGGGACAACTTTGCACTGGCAACGGGTGGCAGTACCGGGAGAAggcggc                                                                                                                                                                                                               >  2:119490/1‑244 (MQ=255)
                                        gCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGACATTACAGAGCGGCAAAGCGGGCACCGTTGATATCCAAAACCAACTGACGGAAAAGACAATGTTGCACTGGGACGGGCTTGGAGTTCCggggggaggccgccgcggg                                                                                                                                                                                                         >  2:477846/1‑235 (MQ=255)
                                         cccAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGcc                                                                                                                                                                                                        <  2:449552/251‑1 (MQ=255)
                                          ccAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGccc                                                                                                                                                                                                       <  2:246718/251‑1 (MQ=255)
                                                            gATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGa                                                                                                                                                                                                <  2:88078/240‑1 (MQ=255)
                                                            gATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGcc                                                                                                                                                                                     <  1:486570/251‑1 (MQ=255)
                                                            gATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGa                                                                                                                                                                                                >  1:88078/1‑240 (MQ=255)
                                                                 gCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGt                                                                                                                                                                                 <  1:135134/250‑1 (MQ=255)
                                                                 gCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGt                                                                                                                                                                                 <  2:142657/250‑1 (MQ=255)
                                                                     aCGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCaa                                                                                                                                                                            <  1:525445/251‑1 (MQ=255)
                                                                                          aTTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGGCCGCAGGGGATTATTCCGCCAGGTGGCAAGCGCTCGGGGacggtggacgt                                                                                                                                                       >  1:407461/1‑251 (MQ=255)
                                                                                            tCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGacgttgaacgtt                                                                                                                                                      <  1:432268/250‑1 (MQ=255)
                                                                                               gTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATc                                                                                                                                                  <  1:151896/251‑1 (MQ=255)
                                                                                                  aaCTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAAc                                                                                                                                               <  1:18814/251‑1 (MQ=255)
                                                                                                      aTTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAAATATTCCGCCCGGTGGCAAGCGCTCGGGGACGTTGAACGTTGATCAACCTGc                                                                                                                                           >  1:694183/1‑251 (MQ=255)
                                                                                                              nggCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGACGGAAGAGACAAAGTTGCACTGGGACGGGCTGGGAGTAACGGGGGAAGTCGACGGCGGCCCGCAGGGGATTATTCCGCCCGGTGGCAAGCGCCCGGGGACGGTGAACGTTGATCAACCTGACGCTAccc                                                                                                                                   >  2:314058/2‑250 (MQ=255)
                                                                                                                gCCAGTCCACCTTTGGCGGGTAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTc                                                                                                                                                                      <  1:78910/214‑1 (MQ=255)
                                                                                                                      ccACCTTGGGCGGGAAAACTACAACTACCTGGGGCTATAACGGGAATCTGCGGGGGCGGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACACCCAACTGACGGAAGAGACAATGTGGCACTGGCACGGGCTGGAAGTCCCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGtt                                                                                                                            <  2:691076/250‑1 (MQ=255)
                                                                                                                        cccTTTGGCGGGAAAGCTGCAACTCCCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCa                                                                                                                         <  2:686553/250‑1 (MQ=255)
                                                                                                                               gggggAAAACTGCCACTCCCTGGGGCTATAACGGTAATCGTGTGGGGCCGGCGGTGAAATTACAGCCCGGGAAAGCGGTAACGTTTCATATCTACACCCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCGGGAAGTACCGGGTGAAGTCGACGGCGCCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGTTGACGTTGAACGTGGATCAGCCTGCCGCTACCTGCTTGTTCCATCCTCAt                                                                                                                  <  2:382249/249‑1 (MQ=255)
                                                                                                                               gCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGAcggcg                                                                                                                                                                                                           <  1:188469/162‑1 (MQ=255)
                                                                                                                               gCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAACGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGAcggcg                                                                                                                                                                                                           >  2:188469/1‑162 (MQ=255)
                                                                                                                                  ggAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGACATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGGACCGGGGGAAGTCGCCGGCGGGCCGCAGGGAATTATTCCGCCAGGTGGGAAGCGCTCGGGGACGTTGAAAGTTGATCAACCTGCCGCCACCTGCTGGTTCCCTCCGcctc                                                                                                                 >  2:372283/1‑246 (MQ=255)
                                                                                                                                                aCCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCCTCAGCACGGCAAAACCgg                                                                                                 >  1:505865/1‑251 (MQ=255)
                                                                                                                                                 ccTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCggg                                                                                                <  1:429539/251‑1 (MQ=255)
                                                                                                                                                  cgggggCTATACCGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTGCCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCACCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGc                                                                                               <  2:127120/249‑1 (MQ=255)
                                                                                                                                                   tGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGGGGCAAGCGCTCGGGGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCg                                                                                              >  1:2839/1‑251 (MQ=255)
                                                                                                                                                       gCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCCCGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGGGACGTTGAAAGGTTATCAACCGGCCGCTACCTGCTGGGTTCATCCGCATCAGCACGGCCAAACCggggggc                                                                                            >  2:44371/1‑245 (MQ=255)
                                                                                                                                                           taACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAggggg                                                                                      >  1:123712/1‑248 (MQ=255)
                                                                                                                                                           taACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAggggg                                                                                      >  1:648534/1‑248 (MQ=255)
                                                                                                                                                            aaCGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTgg                                                                                      <  1:335079/250‑1 (MQ=255)
                                                                                                                                                                 cAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATgg                                                                                >  1:443655/1‑251 (MQ=255)
                                                                                                                                                                  aaTCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAACGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCCGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGGCAGGTGGCGATggg                                                                               >  1:78805/1‑251 (MQ=255)
                                                                                                                                                                     ctgctgGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGGAACGGGTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCCCTGGCCCGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCCAAACCGGGCGGCAGGGGGCGGTGGggcg                                                                            >  2:679015/1‑250 (MQ=255)
                                                                                                                                                                       tgctgGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGggctg                                                                           >  1:9461/1‑251 (MQ=255)
                                                                                                                                                                        gctgGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGGGAAGGCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGGGACGGTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGGCCGGTGGGGGTGGggcggg                                                                          >  2:140615/1‑248 (MQ=255)
                                                                                                                                                                             ggCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTg                                                                                                                                            <  2:17950/180‑1 (MQ=255)
                                                                                                                                                                             ggCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTg                                                                                                                                            >  1:17950/1‑180 (MQ=255)
                                                                                                                                                                                 ggcggTGAAATTACAGCGCGGCCAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGgctgggctgg                                                                 >  1:543564/1‑251 (MQ=255)
                                                                                                                                                                                  gcggTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCtggt                                                                <  1:263362/251‑1 (MQ=255)
                                                                                                                                                                                     gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCACCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTggtggt                                                             <  1:601531/251‑1 (MQ=255)
                                                                                                                                                                                     gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGTTGGGCTggtggt                                                             <  1:481637/251‑1 (MQ=255)
                                                                                                                                                                                                gcgcGGCAAAGCGGTAACGGTTGATATCTGCAACCAACTGACGGAAGAGACAAGTTTGCGCAGGCACGGGCTGGAAGTACCGGGTGAATTCGACGGGGGCCCGCAGGGAATTATCCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                  <  1:443576/251‑1 (MQ=255)
                                                                                                                                                                                                   cgGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGCGGTGATTGAAGATGACg                                               >  2:53184/1‑251 (MQ=255)
                                                                                                                                                                                                     gCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGTTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGGTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                             <  1:477846/251‑1 (MQ=255)
                                                                                                                                                                                                     gCAAAGCGGTAACGGTGGATATCAACAACCAACTGACGGAAGAGACAATGTTCCCCTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGGGGCAAGCGCTCGGGGACGTGGAACGTTGATCACCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                             <  2:558029/251‑1 (MQ=255)
                                                                                                                                                                                                      aaaaGCGGGAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCGGGAAGTCCCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCACGCGCTCGGTGACGTTGAAGTTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGGGGTGATTGAAGATGACgaga                                            <  2:505865/250‑1 (MQ=255)
                                                                                                                                                                                                            ggTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACct                                                                                                                                   >  2:285260/1‑158 (MQ=255)
                                                                                                                                                                                                            ggTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACct                                                                                                                                   <  1:285260/158‑1 (MQ=255)
                                                                                                                                                                                                            ggTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGACGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTg                                       >  1:338324/1‑250 (MQ=255)
                                                                                                                                                                                                                aCGGTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAtt                                  >  1:376516/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGGTTGATATCTACAACCAACTGGCGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAttn                                 <  2:310495/251‑2 (MQ=255)
                                                                                                                                                                                                                   gTTGATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAt                               >  2:442651/1‑251 (MQ=255)
                                                                                                                                                                                                                      gATATCTACAACCAACTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgct                            >  2:140287/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCAACCAACTGACGGAAGAGAAAATGTTGCACTGGCACGGGCTGGAAGTACCGGTTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCTTGAAATTAATGCTGCCaaaa                     <  1:119490/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaccaacTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                <  1:15914/222‑1 (MQ=255)
                                                                                                                                                                                                                               aaccaacTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                >  2:15914/1‑222 (MQ=255)
                                                                                                                                                                                                                                accaacTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAAc                    >  1:583784/1‑249 (MQ=255)
                                                                                                                                                                                                                                 ccaacTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGt                 <  1:222804/251‑1 (MQ=255)
                                                                                                                                                                                                                                  caacTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCAACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACgg                                                                                                          >  2:192714/1‑161 (MQ=255)
                                                                                                                                                                                                                                   aacTGACGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGCACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGCGATCCTGAAATTAATGCTGCCAAAACAGTgg               >  1:375669/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATc         <  1:53184/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cGGAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCCGCCAAAACAGTGGGGTATc         >  2:429263/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ngAAGAGACAATGTTACACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCAGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg        >  2:326604/2‑251 (MQ=255)
                                                                                                                                                                                                                                           ggAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGACGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg        <  2:338324/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAAGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTCGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg        >  1:48841/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaGAGACAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgat      >  1:146952/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gaCAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg  >  2:461085/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gaCAATGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg  <  1:294853/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATCTGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGACGGAAGAGACAACGTTGCACTGGCACGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATG  >  NC_000913/137134‑137623

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