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

GCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAACGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATG  >  NC_000913/137134‑137623
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCGAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGc                                                                                                                                                                                                                                                   >  1:691126/1‑250 (MQ=255)
gCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGGGAAAATACAGCGCGGCAAAGCGGGAAAGGTTGATATCTAAAACCCACTGA‑CGGGAGAGACAAAGTTgg                                                                                                                                                                                                                                                   >  2:141987/1‑249 (MQ=255)
        gCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑Cg                                                                                                                                                                                                                                         <  1:467058/251‑1 (MQ=255)
           gtgGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCGAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGc                                                                                                                                                                                                                                      <  2:691126/251‑1 (MQ=255)
             ggAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCt                                                                                                                                                                                                                                     <  1:150895/250‑1 (MQ=255)
                 ccGCGCAGTATTTGCGGCAGAACGCCCAGCGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAg                                                                                                                                                                                                                                <  1:474858/251‑1 (MQ=255)
                    cgcAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTAc                                                                                                                                                                                                                             <  2:210452/251‑1 (MQ=255)
                    cgcAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGGGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAAATGA‑CGGGAGAGACAATGTTGCACTGGCC‑AGGGCTGGGAGTAc                                                                                                                                                                                                                             >  2:160851/1‑251 (MQ=255)
                          aTTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGt                                                                                                                                                                                                                        <  1:264924/250‑1 (MQ=255)
                               cGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGGGAAGTc                                                                                                                                                                                                                  >  1:246718/1‑251 (MQ=255)
                                  cAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGAc                                                                                                                                                                                                               <  1:141987/251‑1 (MQ=255)
                                    gAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTTGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAAAGGCAATC‑TGCTGGGGCCGGCGGGGAAAATAAAGCGCGGGAAAACGGCAACGGTTGGTATCTACAACCACATGAACGGA‑GGGACAACTTTGCACTGGCAACGGG‑TGGCAGTACCGGGAGAAggcggc                                                                                                                                                                                                               >  2:119490/1‑244 (MQ=255)
                                        gCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGACATTACAGAGCGGCAAAGCGGGCACCGTTGATATCCAAAACCAACTGA‑CGGAAAAGACAATGTTGCACTGGGA‑CGGGCTTGGAGTTCCggggggaggccgccgcggg                                                                                                                                                                                                         >  2:477846/1‑235 (MQ=255)
                                         cccAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGcc                                                                                                                                                                                                        <  2:449552/251‑1 (MQ=255)
                                          ccAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGccc                                                                                                                                                                                                       <  2:246718/251‑1 (MQ=255)
                                                            gATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGa                                                                                                                                                                                                <  2:88078/240‑1 (MQ=255)
                                                            gATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGcc                                                                                                                                                                                     <  1:486570/251‑1 (MQ=255)
                                                            gATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGa                                                                                                                                                                                                >  1:88078/1‑240 (MQ=255)
                                                                 gCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGt                                                                                                                                                                                 <  1:135134/250‑1 (MQ=255)
                                                                 gCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGt                                                                                                                                                                                 <  2:142657/250‑1 (MQ=255)
                                                                     aCGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCaa                                                                                                                                                                            <  1:525445/251‑1 (MQ=255)
                                                                                          aTTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGGCCGCAGGGGATTATTCCGCCAGGTGGCAAGCGCTCGGGGacggtggacgt                                                                                                                                                       >  1:407461/1‑251 (MQ=255)
                                                                                            tCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGacgttgaacgtt                                                                                                                                                      <  1:432268/250‑1 (MQ=255)
                                                                                               gTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATc                                                                                                                                                  <  1:151896/251‑1 (MQ=255)
                                                                                                  aaCTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAAc                                                                                                                                               <  1:18814/251‑1 (MQ=255)
                                                                                                      aTTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAAATATTCCGCCCGGTGGCAAGCGCTCGGGGACGTTGAACGTTGATCAACCTGc                                                                                                                                           >  1:694183/1‑251 (MQ=255)
                                                                                                              nggCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGA‑CGGAAGAGACAAAGTTGCACTGGGA‑CGGGCTGGGAGTAACGGGGGAAGTCGACGGCGGCCCGCAGGGGATTATTCCGCCCGGTGGCAAGCGCCCGGGGACGGTGAACGTTGATCAACCTGACGCTAccc                                                                                                                                   >  2:314058/2‑250 (MQ=255)
                                                                                                                gCCAGTCCACCTTTGGCGGGTAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTc                                                                                                                                                                      <  1:78910/214‑1 (MQ=255)
                                                                                                                      ccACCTTGGGCGGGAAAACTACAACTACCTGGGGCTATAACGGGAATC‑TGCGGGGGCGGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACACCCAACTGA‑CGGAAGAGACAATGTGGCACTGGCA‑CGGGCTGGAAGTCCCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGtt                                                                                                                            <  2:691076/250‑1 (MQ=255)
                                                                                                                        cccTTTGGCGGGAAAGCTGCAACTCCCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCa                                                                                                                         <  2:686553/250‑1 (MQ=255)
                                                                                                                               gggggAAAACTGCCACTCCCTGGGGCTATAACGGTAATCGTG‑TGGGGCCGGCGGTGAAATTACAGCCCGGGAAAGCGGTAACGTTTCATATCTACACCCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCGGGAAGTACCGGGTGAAGTCGACGGCGCCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGTTGACGTTGAACGTGGATCAGCCTGCCGCTACCTGCTTGTTCCATCCTCAt                                                                                                                  <  2:382249/249‑1 (MQ=255)
                                                                                                                               gCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGAcggcg                                                                                                                                                                                                           <  1:188469/162‑1 (MQ=255)
                                                                                                                               gCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAACGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGAcggcg                                                                                                                                                                                                           >  2:188469/1‑162 (MQ=255)
                                                                                                                                  ggAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGACATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGGACCGGGGGAAGTCGCCGGCGGGCCGCAGGGAATTATTCCGCCAGGTGGGAAGCGCTCGGGGACGTTGAAAGTTGATCAACCTGCCGCCACCTGCTGGTTCCCTCCGcctc                                                                                                                 >  2:372283/1‑246 (MQ=255)
                                                                                                                                                aCCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCCTCAGCACGGCAAAACCgg                                                                                                 >  1:505865/1‑251 (MQ=255)
                                                                                                                                                 ccTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCggg                                                                                                <  1:429539/251‑1 (MQ=255)
                                                                                                                                                  cgggggCTATACCGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTGCCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCACCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGc                                                                                               <  2:127120/249‑1 (MQ=255)
                                                                                                                                                   tGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGGGGCAAGCGCTCGGGGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCg                                                                                              >  1:2839/1‑251 (MQ=255)
                                                                                                                                                       gCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCC‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGGGACGTTGAAAGGTTATCAACCGGCCGCTACCTGCTGGGTTCATCCGCATCAGCACGGCCAAACCggggggc                                                                                            >  2:44371/1‑245 (MQ=255)
                                                                                                                                                           taACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAggggg                                                                                      >  1:123712/1‑248 (MQ=255)
                                                                                                                                                           taACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGGTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAggggg                                                                                      >  1:648534/1‑248 (MQ=255)
                                                                                                                                                            aaCGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTgg                                                                                      <  1:335079/250‑1 (MQ=255)
                                                                                                                                                                 cAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATgg                                                                                >  1:443655/1‑251 (MQ=255)
                                                                                                                                                                  aaTC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAACGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCCGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGGCAGGTGGCGATggg                                                                               >  1:78805/1‑251 (MQ=255)
                                                                                                                                                                     c‑tgctgGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGGAACGGGTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCCCTGGCC‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCCAAACCGGGCGGCAGGGGGCGGTGGggcg                                                                            >  2:679015/1‑250 (MQ=255)
                                                                                                                                                                       tgctgGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGggctg                                                                           >  1:9461/1‑251 (MQ=255)
                                                                                                                                                                        gctgGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGGGAAGGCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGGGACGGTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGGCCGGTGGGGGTGGggcggg                                                                          >  2:140615/1‑248 (MQ=255)
                                                                                                                                                                             ggCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTg                                                                                                                                            <  2:17950/180‑1 (MQ=255)
                                                                                                                                                                             ggCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTg                                                                                                                                            >  1:17950/1‑180 (MQ=255)
                                                                                                                                                                                 ggcggTGAAATTACAGCGCGGCCAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGgctgggctgg                                                                 >  1:543564/1‑251 (MQ=255)
                                                                                                                                                                                  gcggTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCtggt                                                                <  1:263362/251‑1 (MQ=255)
                                                                                                                                                                                     gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCACCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTggtggt                                                             <  1:601531/251‑1 (MQ=255)
                                                                                                                                                                                     gTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGTTGGGCTggtggt                                                             <  1:481637/251‑1 (MQ=255)
                                                                                                                                                                                                gcgcGGCAAAGCGGTAACGGTTGATATCTGCAACCAACTGA‑CGGAAGAGACAAGTTTGCGCAGGCA‑CGGGCTGGAAGTACCGGGTGAATTCGACGGGGGCCCGCAGGGAATTATCCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATg                                                  <  1:443576/251‑1 (MQ=255)
                                                                                                                                                                                                   cgGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGCGGTGATTGAAGATGACg                                               >  2:53184/1‑251 (MQ=255)
                                                                                                                                                                                                     gCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGTTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGGTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                             <  1:477846/251‑1 (MQ=255)
                                                                                                                                                                                                     gCAAAGCGGTAACGGTGGATATCAACAACCAACTGA‑CGGAAGAGACAATGTTCCCCTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGGGGCAAGCGCTCGGGGACGTGGAACGTTGATCACCCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACgag                                             <  2:558029/251‑1 (MQ=255)
                                                                                                                                                                                                      aaaaGCGGGAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCGGGAAGTCCCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCACGCGCTCGGTGACGTTGAAGTTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGGGGTGATTGAAGATGACgaga                                            <  2:505865/250‑1 (MQ=255)
                                                                                                                                                                                                            ggTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACct                                                                                                                                   >  2:285260/1‑158 (MQ=255)
                                                                                                                                                                                                            ggTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACct                                                                                                                                   <  1:285260/158‑1 (MQ=255)
                                                                                                                                                                                                            ggTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGACGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTg                                       >  1:338324/1‑250 (MQ=255)
                                                                                                                                                                                                                aCGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAtt                                  >  1:376516/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGGTTGATATCTACAACCAACTGG‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAAttn                                 <  2:310495/251‑2 (MQ=255)
                                                                                                                                                                                                                   gTTGATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAt                               >  2:442651/1‑251 (MQ=255)
                                                                                                                                                                                                                      gATATCTACAACCAACTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAAtgct                            >  2:140287/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCAACCAACTGA‑CGGAAGAGAAAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGTTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCTTGAAATTAATGCTGCCaaaa                     <  1:119490/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaccaacTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                <  1:15914/222‑1 (MQ=255)
                                                                                                                                                                                                                               aaccaacTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGAc                                                >  2:15914/1‑222 (MQ=255)
                                                                                                                                                                                                                                accaacTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAAc                    >  1:583784/1‑249 (MQ=255)
                                                                                                                                                                                                                                 ccaacTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGt                 <  1:222804/251‑1 (MQ=255)
                                                                                                                                                                                                                                  caacTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCAACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACgg                                                                                                          >  2:192714/1‑161 (MQ=255)
                                                                                                                                                                                                                                   aacTGA‑CGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGCACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGCGATCCTGAAATTAATGCTGCCAAAACAGTgg               >  1:375669/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATc         <  1:53184/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cGGAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCCGCCAAAACAGTGGGGTATc         >  2:429263/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ngAAGAGACAATGTTACACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCAGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg        >  2:326604/2‑251 (MQ=255)
                                                                                                                                                                                                                                           ggAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGACGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg        <  2:338324/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAAGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTCGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCg        >  1:48841/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaGAGACAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCgat      >  1:146952/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gaCAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg  >  2:461085/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gaCAATGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGatgatg  <  1:294853/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCGGCAGAACGCCCAACGTTACCGATCCCTGATTTGCTCACGACCGATGCCCGTAATCGCATTCAGTTAACTATTGGCGCAGGCCAGTCCACCTTTGGCGGGAAAACTGCAACTACCTGGGGCTATAACGGCAATC‑TGCTGGGGCCGGCGGTGAAATTACAGCGCGGCAAAGCGGTAACGGTTGATATCTACAACCAACTGA‑CGGAAGAGACAACGTTGCACTGGCA‑CGGGCTGGAAGTACCGGGTGAAGTCGACGGCGGCCCGCAGGGAATTATTCCGCCAGGTGGCAAGCGCTCGGTGACGTTGAACGTTGATCAACCTGCCGCTACCTGCTGGTTCCATCCGCATCAGCACGGCAAAACCGGGCGACAGGTGGCGATGGGGCTGGCTGGGCTGGTGGTGATTGAAGATGACGAGATCCTGAAATTAATGCTGCCAAAACAGTGGGGTATCGATGATG  >  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: ‑