Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,069 T→C G249G (GGT→GGC)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,0690TC92.5% 220.6 / 1.3 80G249G (GGT→GGC) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (32/42);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCT  >  NC_000913/4617823‑4618301
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
gctcgtccgagctAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGAGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑Cg                                                                                                                                                                                                                                      <  1:375903/239‑1 (MQ=255)
           gTAAAGTGTCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAag                                                                                                                                                                                                                           <  1:58816/251‑1 (MQ=255)
             aaaGTGGGTGTGAACGCGACGCCGGAAAGCGCGAGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGGTTAAAACCGGCGGGCGGCGTGCGTACTGCGGTAGATGAGCAGAAATATATCGCTATCGCAGATGAACTGTTCGGTGCGGACTGGGCAGATGCGTGTCACTTCCGCTTTGGTGCTTCCAGCCTGCGGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAagag                                                                                                                                                                                                                         <  2:193063/251‑1 (MQ=255)
               aGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGTCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAgcg                                                                                                                                                                                                                       <  1:219610/251‑1 (MQ=255)
                  ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGcc                                                                                                                                                                                                                     >  1:45776/1‑250 (MQ=255)
                  ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGGAGAAAAAACCGTTTGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGGGCTGAGTGGGCAGATGCGCGTCACTACCGCTTTGGGGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCca                                                                                                                                                                                                                    <  1:276554/251‑1 (MQ=255)
                                 gCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACtaagtaa                                                                                                                                                                                                     <  1:87709/251‑1 (MQ=255)
                                     gcaaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGAGCTTCCAGCCTGCTGGCAAGCCTGTTGACAGCGCTGGGCCA‑CGGTGACGGTACGAACGCCAGCAGCTACTAAGTAAGATg                                                                                                                                                                                                 >  2:430494/3‑251 (MQ=255)
                                     gAAAGCGCGCGCATCATGATGGCAGGGATCCGTGATCTGGGCGTAGAAACAACCGTTGGGTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCCCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCT‑CGGTGCCGGTACGAGCACCCGCAGCTCCtaagaacaatg                                                                                                                                                                                                 >  2:388835/1‑244 (MQ=255)
                                     gAAAGCGCGCGCATCATGATGGAAGTGATCCTTGATATGGGCGTAGAAAAAACCGTTGTTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAGGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGTTCTGACTGGTCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAATTAAGa                                                                                                                                                                                                   <  2:557682/249‑1 (MQ=255)
                                       aaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCt                                                                                                                                                                                               <  1:523746/251‑1 (MQ=255)
                                                 atcatGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGa                                                                                                                                                                                     <  1:430494/251‑1 (MQ=255)
                                                      gatgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgc                                                                                                                                                                                <  1:71953/251‑1 (MQ=255)
                                                       atgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                                                                                                                                                               >  1:466643/1‑251 (MQ=255)
                                                         ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑Cggtga                                                                                                                                                                                                                                  >  1:483439/1‑195 (MQ=255)
                                                         ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTGGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCCCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTCCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑Cggtga                                                                                                                                                                                                                                  <  2:483439/198‑4 (MQ=255)
                                                                aTCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATc                                                                                                                                                                      >  1:158222/1‑251 (MQ=255)
                                                                   cGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAgg                                                                                                                                                                   >  1:445083/1‑251 (MQ=255)
                                                                    gTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGCAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAACGCTGGGCCA‑CGGTGACGGTAAGCGCGCCAGCCGCTACTCAGGAAGATGCTTTACGCCTGATCCGCTGCGCTTATCAGGc                                                                                                                                                                  >  2:5349/1‑251 (MQ=255)
                                                                     ggATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGcc                                                                                                                                                                 <  1:16165/250‑1 (MQ=255)
                                                                      gATATGGGCGTAGAAAAAACAGTTGGATTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGAGGCGAGTCACTACCGCTTTGGTGCTTACAACCTGATGGCAAGCCTGTTGAATGCGCTGG‑TCAACGGTGACGGTAAGAGCGCAAGCAGCTACTAAGGAAGTTGCTTTACGGCTGATGAGCTGCGCTTATCacgctt                                                                                                                                                                >  2:116612/1‑246 (MQ=255)
                                                                             gCGTAGAAAAAACCGTGGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGAc                                                                                                                                                         <  1:388835/251‑1 (MQ=255)
                                                                                 aGAAAAAACCGTTGTTTTCAAACCGCCGGGAGGCGTGCGTACGGCGGAAGATGCGCAGAAATATCTCGCCAGTGCAGATGATCTGGTCGGTGCTGACTGGGCAGATACGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTGCTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAt                                                                                                                                                     <  2:16722/251‑1 (MQ=255)
                                                                                 aGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTa                                                                                                                                                      <  1:132326/250‑1 (MQ=255)
                                                                                    aaaaaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCt                                                                                                                                                   >  1:361568/1‑250 (MQ=255)
                                                                                       aaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAccc                                                                                                                                               >  1:649923/1‑251 (MQ=255)
                                                                                       aaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAccc                                                                                                                                               <  1:495002/251‑1 (MQ=255)
                                                                                        aCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTCCTGCGGAAGATGCGCAGAAATATCGCGCCATCGCAGATTAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCg                                                                                                                                              <  2:571436/251‑1 (MQ=255)
                                                                                           gTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAGCTACCCGTAg                                                                                                                                           <  2:649923/251‑1 (MQ=255)
                                                                                            ttGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGGTGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTCCCCGTAgg                                                                                                                                          >  1:416489/1‑251 (MQ=255)
                                                                                               gTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCg                                                                                                                                       <  1:213778/251‑1 (MQ=255)
                                                                                                       ccagcgggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCggtaaaggc                                                                                                                               >  2:581628/4‑244 (MQ=255)
                                                                                                            ggtcgtcgggcgtACTGCGTATGATGCGCAGAAACATCTCGCCACCGCAGATGAACTGTCAGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAATCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTATGCCGGATAAGGCGTAGa                                                                                                                          <  2:290454/242‑1 (MQ=255)
                                                                                                            gggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcg                                                                                                                                                                                 <  2:58516/196‑1 (MQ=255)
                                                                                                            gggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcg                                                                                                                                                                                 >  1:58516/1‑196 (MQ=255)
                                                                                                            gggcggAGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGTCAGATGCGCGGCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTTCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGa                                                                                                                          <  1:402376/251‑1 (MQ=255)
                                                                                                             ggcggcGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGAc                                                                                                                         <  1:599119/251‑1 (MQ=255)
                                                                                                                ggcgTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGCCTGGGCAGATGCGCGTCCCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTa                                                                                                                                                                        >  1:462021/1‑201 (MQ=255)
                                                                                                                    tgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCg                                                                                                                  >  2:612782/1‑251 (MQ=255)
                                                                                                                      cgtACTGCGGCAGCTGCGCAGCAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATACGCGTCACTACCGATTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGCTAAGGCGTAGACGCATCCGGc                                                                                                                >  2:443681/1‑251 (MQ=255)
                                                                                                                         aCTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCaa                                                                                                              >  1:127583/1‑250 (MQ=255)
                                                                                                                             cGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                                                                                                         <  1:542510/251‑1 (MQ=255)
                                                                                                                             cGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                                                                                                         <  1:150351/251‑1 (MQ=255)
                                                                                                                             cGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                                                                                                         <  1:589358/251‑1 (MQ=255)
                                                                                                                              ggAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTAGGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccg                                                                                                        <  1:309085/251‑1 (MQ=255)
                                                                                                                                          gAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCtt                                                                                            <  1:62497/251‑1 (MQ=255)
                                                                                                                                                tctcGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTAAGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATACGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTc                                                                                      >  1:248376/1‑251 (MQ=255)
                                                                                                                                                   cGCCATCGCAGATGGACTGTTCGGTGATGACTGGGCAGATGGGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTc                                                                                      <  1:124547/248‑1 (MQ=255)
                                                                                                                                                      catcgCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCgggcgg                                                                                >  2:155270/5‑248 (MQ=255)
                                                                                                                                                        gcgCAGATGACCGGTTCGGGGCTGACTGGGCAGATGCGCGTCACTACCGCTTGGGTGCGTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGtt                                                                              <  2:639995/249‑1 (MQ=255)
                                                                                                                                                           cAGATGAACTGTTCGGTGCTGACTGGGCAGCTGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACtaagt                                                                                                                                                                                                       >  2:319019/1‑127 (MQ=255)
                                                                                                                                                           cAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACtaagt                                                                                                                                                                                                       <  1:319019/127‑1 (MQ=255)
                                                                                                                                                             gATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCtt                                                                         <  1:183457/251‑1 (MQ=255)
                                                                                                                                                              aTGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTCCCGCTTTGGTGCTTCCAGCCTACTGGCACGACTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGAAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGAATCTACCCGTAGGCCGGGTAAGGAGTAGACGCATCCGGCAACAGCCGCCTCATACTCTTTTCATCGGGCGGTTACCTTg                                                                        >  1:242424/1‑251 (MQ=255)
                                                                                                                                                                 aatggtgcGGTGCTGACGGGGCAGATGCGGGGCCGTACCGCTGTGGAGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGACCGCTAGCAGTTACTAAGTAAGATGCTTTAAGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGttt                                                                     <  2:248376/244‑1 (MQ=255)
                                                                                                                                                                   catttcgggtcgtaccgggCAGAGACGCGTCGCTACCGCTTTGGTGCGTCCCGCCTGCTGGCAAGCCTGTTGAAGGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGAGGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTATGAGACGTATCTTCCCGTAGGCCGGATAAGGCGTAGACGCATTCGGCAAAAGCCGCCTCCTACTCTTTTCCTCGGGAGGTTACCTTGTTtct                                                                   <  2:326521/235‑1 (MQ=255)
                                                                                                                                                                     gTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCg                                                                 >  1:335582/1‑251 (MQ=255)
                                                                                                                                                                         ggTGCTGACTGGGCAGATGAGCGTCACGACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCGGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCgg                                                                                    <  1:639027/228‑1 (MQ=255)
                                                                                                                                                                         ggTGCTGACTGGGCAGAGGCGCGTCACTCACGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑AGTTTACGGTCAGCGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACTAGACGTCTCTCCCCACAGTACGGATACGGCGTAGACGACTCCGGCAACAGCCCCCTCATACTCTTTTCCTCgt                                                                                    >  2:639027/1‑227 (MQ=255)
                                                                                                                                                                            gCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGa                                                          <  2:127583/251‑1 (MQ=255)
                                                                                                                                                                            gCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGa                                                          <  1:612782/251‑1 (MQ=255)
                                                                                                                                                                               gccTGGGCAGATGCGAGGCACCACCGCTGTGGTGCGTCCATCCCGCTGGCAAGCCTGTTGAAAGCGTTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTCAGTCAGATGCTTTTCGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGCAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATTCTCTTTTCCTCGGGAGGTTACCTTGTTACGCGCACAAGAAat                                                       <  2:275555/249‑1 (MQ=255)
                                                                                                                                                                               gACTTGGCAGGTGCGCGTCACCACTGATTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGaaa                                                        <  2:281475/250‑1 (MQ=255)
                                                                                                                                                                                    ggCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGCTAAGGCGTAGACGCATCCGGCAAACGCCGCATCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTc                                                  >  2:534487/1‑251 (MQ=255)
                                                                                                                                                                                     gCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCg                                                 >  1:361341/1‑251 (MQ=255)
                                                                                                                                                                                                  ccTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGAt                                    <  1:72020/250‑1 (MQ=255)
                                                                                                                                                                                                  aCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGAt                                    <  1:5349/251‑1 (MQ=255)
                                                                                                                                                                                                  aCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGAt                                    <  1:44935/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGACAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCAAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCAGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGGAATTATTCGTAAAAAACGTGATGGTc                                >  1:403130/1‑251 (MQ=255)
                                                                                                                                                                                                           ttGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATg                             <  1:443681/249‑1 (MQ=255)
                                                                                                                                                                                                           ttGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGAAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATgc                            >  1:658054/1‑250 (MQ=255)
                                                                                                                                                                                                             ggtgcTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGAAGCTTTACGCCTGATGCGCTGAGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTACTACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCAACAGAAATTATTCGTAAAACACGTGATGGTCaggaact                         >  2:649906/4‑245 (MQ=255)
                                                                                                                                                                                                              gtgcTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTg                        >  1:463604/3‑251 (MQ=255)
                                                                                                                                                                                                                           gctgGCAAGCCTGTTGAAAGCGCAGGGCCA‑CGGTGACGGTACGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACACGAAATTATTAGTAAAAAACGTGATGGTAATGCACTGAGCGATgaagaa            >  2:6262/1‑250 (MQ=255)
                                                                                                                                                                                                                              ggCAAGCCTGTTGAAAGCGCTGGGCCT‑CGGTGACGTTAAGAGCGCCAGCAGCTACTAATTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAAtt         <  1:624689/250‑1 (MQ=255)
                                                                                                                                                                                                                               gCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCACAAGCCGCCTCATACTCTTTTCCTCGGGAg                                                                                 >  2:1408/1‑177 (MQ=255)
                                                                                                                                                                                                                                 aaGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGtt     >  2:620202/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cgacgccggAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATgcgctgcgct                                                                                                                                                                          >  1:46334/9‑85 (MQ=255)
                                                                                                                                                                                                                                  cgacgccggAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATgcgctgcgct                                                                                                                                                                          <  2:46334/77‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCt  >  1:166121/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
AACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCT  >  NC_000913/4617823‑4618301

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