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→GGCdeoCdeoxyribose‑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

AACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCT  >  NC_000913/4617823‑4618301
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
gctcgtccgagctAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGAGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACg                                                                                                                                                                                                                                      <  1:375903/239‑1 (MQ=255)
           gTAAAGTGTCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAag                                                                                                                                                                                                                           <  1:58816/251‑1 (MQ=255)
             aaaGTGGGTGTGAACGCGACGCCGGAAAGCGCGAGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGGTTAAAACCGGCGGGCGGCGTGCGTACTGCGGTAGATGAGCAGAAATATATCGCTATCGCAGATGAACTGTTCGGTGCGGACTGGGCAGATGCGTGTCACTTCCGCTTTGGTGCTTCCAGCCTGCGGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAagag                                                                                                                                                                                                                         <  2:193063/251‑1 (MQ=255)
               aGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGTCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAgcg                                                                                                                                                                                                                       <  1:219610/251‑1 (MQ=255)
                  ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGcc                                                                                                                                                                                                                     >  1:45776/1‑250 (MQ=255)
                  ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGGAGAAAAAACCGTTTGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGGGCTGAGTGGGCAGATGCGCGTCACTACCGCTTTGGGGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCca                                                                                                                                                                                                                    <  1:276554/251‑1 (MQ=255)
                                 gCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACtaagtaa                                                                                                                                                                                                     <  1:87709/251‑1 (MQ=255)
                                     gcaaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGAGCTTCCAGCCTGCTGGCAAGCCTGTTGACAGCGCTGGGCCACGGTGACGGTACGAACGCCAGCAGCTACTAAGTAAGATg                                                                                                                                                                                                 >  2:430494/3‑251 (MQ=255)
                                     gAAAGCGCGCGCATCATGATGGCAGGGATCCGTGATCTGGGCGTAGAAACAACCGTTGGGTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCCCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCTCGGTGCCGGTACGAGCACCCGCAGCTCCtaagaacaatg                                                                                                                                                                                                 >  2:388835/1‑244 (MQ=255)
                                     gAAAGCGCGCGCATCATGATGGAAGTGATCCTTGATATGGGCGTAGAAAAAACCGTTGTTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAGGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGTTCTGACTGGTCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAATTAAGa                                                                                                                                                                                                   <  2:557682/249‑1 (MQ=255)
                                       aaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCt                                                                                                                                                                                               <  1:523746/251‑1 (MQ=255)
                                                 atcatGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGa                                                                                                                                                                                     <  1:430494/251‑1 (MQ=255)
                                                      gatgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgc                                                                                                                                                                                <  1:71953/251‑1 (MQ=255)
                                                       atgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                                                                                                                                                               >  1:466643/1‑251 (MQ=255)
                                                         ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACggtga                                                                                                                                                                                                                                  >  1:483439/1‑195 (MQ=255)
                                                         ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTGGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCCCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTCCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACggtga                                                                                                                                                                                                                                  <  2:483439/198‑4 (MQ=255)
                                                                aTCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATc                                                                                                                                                                      >  1:158222/1‑251 (MQ=255)
                                                                   cGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAgg                                                                                                                                                                   >  1:445083/1‑251 (MQ=255)
                                                                    gTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGCAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAACGCTGGGCCACGGTGACGGTAAGCGCGCCAGCCGCTACTCAGGAAGATGCTTTACGCCTGATCCGCTGCGCTTATCAGGc                                                                                                                                                                  >  2:5349/1‑251 (MQ=255)
                                                                     ggATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGcc                                                                                                                                                                 <  1:16165/250‑1 (MQ=255)
                                                                      gATATGGGCGTAGAAAAAACAGTTGGATTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGAGGCGAGTCACTACCGCTTTGGTGCTTACAACCTGATGGCAAGCCTGTTGAATGCGCTGGTCAACGGTGACGGTAAGAGCGCAAGCAGCTACTAAGGAAGTTGCTTTACGGCTGATGAGCTGCGCTTATCacgctt                                                                                                                                                                >  2:116612/1‑246 (MQ=255)
                                                                             gCGTAGAAAAAACCGTGGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGAc                                                                                                                                                         <  1:388835/251‑1 (MQ=255)
                                                                                 aGAAAAAACCGTTGTTTTCAAACCGCCGGGAGGCGTGCGTACGGCGGAAGATGCGCAGAAATATCTCGCCAGTGCAGATGATCTGGTCGGTGCTGACTGGGCAGATACGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTGCTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAt                                                                                                                                                     <  2:16722/251‑1 (MQ=255)
                                                                                 aGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTa                                                                                                                                                      <  1:132326/250‑1 (MQ=255)
                                                                                    aaaaaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCt                                                                                                                                                   >  1:361568/1‑250 (MQ=255)
                                                                                       aaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAccc                                                                                                                                               >  1:649923/1‑251 (MQ=255)
                                                                                       aaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAccc                                                                                                                                               <  1:495002/251‑1 (MQ=255)
                                                                                        aCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTCCTGCGGAAGATGCGCAGAAATATCGCGCCATCGCAGATTAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCg                                                                                                                                              <  2:571436/251‑1 (MQ=255)
                                                                                           gTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAGCTACCCGTAg                                                                                                                                           <  2:649923/251‑1 (MQ=255)
                                                                                            ttGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGGTGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTCCCCGTAgg                                                                                                                                          >  1:416489/1‑251 (MQ=255)
                                                                                               gTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCg                                                                                                                                       <  1:213778/251‑1 (MQ=255)
                                                                                                       ccagcgggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCggtaaaggc                                                                                                                               >  2:581628/4‑244 (MQ=255)
                                                                                                            ggtcgtcgggcgtACTGCGTATGATGCGCAGAAACATCTCGCCACCGCAGATGAACTGTCAGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAATCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTATGCCGGATAAGGCGTAGa                                                                                                                          <  2:290454/242‑1 (MQ=255)
                                                                                                            gggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcg                                                                                                                                                                                 <  2:58516/196‑1 (MQ=255)
                                                                                                            gggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcg                                                                                                                                                                                 >  1:58516/1‑196 (MQ=255)
                                                                                                            gggcggAGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGTCAGATGCGCGGCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTTCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGa                                                                                                                          <  1:402376/251‑1 (MQ=255)
                                                                                                             ggcggcGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGAc                                                                                                                         <  1:599119/251‑1 (MQ=255)
                                                                                                                ggcgTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGCCTGGGCAGATGCGCGTCCCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTa                                                                                                                                                                        >  1:462021/1‑201 (MQ=255)
                                                                                                                    tgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCg                                                                                                                  >  2:612782/1‑251 (MQ=255)
                                                                                                                      cgtACTGCGGCAGCTGCGCAGCAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATACGCGTCACTACCGATTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGCTAAGGCGTAGACGCATCCGGc                                                                                                                >  2:443681/1‑251 (MQ=255)
                                                                                                                         aCTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCaa                                                                                                              >  1:127583/1‑250 (MQ=255)
                                                                                                                             cGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                                                                                                         <  1:542510/251‑1 (MQ=255)
                                                                                                                             cGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                                                                                                         <  1:150351/251‑1 (MQ=255)
                                                                                                                             cGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                                                                                                         <  1:589358/251‑1 (MQ=255)
                                                                                                                              ggAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTAGGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccg                                                                                                        <  1:309085/251‑1 (MQ=255)
                                                                                                                                          gAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCtt                                                                                            <  1:62497/251‑1 (MQ=255)
                                                                                                                                                tctcGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTAAGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATACGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTc                                                                                      >  1:248376/1‑251 (MQ=255)
                                                                                                                                                   cGCCATCGCAGATGGACTGTTCGGTGATGACTGGGCAGATGGGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTc                                                                                      <  1:124547/248‑1 (MQ=255)
                                                                                                                                                      catcgCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCgggcgg                                                                                >  2:155270/5‑248 (MQ=255)
                                                                                                                                                        gcgCAGATGACCGGTTCGGGGCTGACTGGGCAGATGCGCGTCACTACCGCTTGGGTGCGTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGtt                                                                              <  2:639995/249‑1 (MQ=255)
                                                                                                                                                           cAGATGAACTGTTCGGTGCTGACTGGGCAGCTGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACtaagt                                                                                                                                                                                                       >  2:319019/1‑127 (MQ=255)
                                                                                                                                                           cAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACtaagt                                                                                                                                                                                                       <  1:319019/127‑1 (MQ=255)
                                                                                                                                                             gATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCtt                                                                         <  1:183457/251‑1 (MQ=255)
                                                                                                                                                              aTGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTCCCGCTTTGGTGCTTCCAGCCTACTGGCACGACTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGAAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGAATCTACCCGTAGGCCGGGTAAGGAGTAGACGCATCCGGCAACAGCCGCCTCATACTCTTTTCATCGGGCGGTTACCTTg                                                                        >  1:242424/1‑251 (MQ=255)
                                                                                                                                                                 aatggtgcGGTGCTGACGGGGCAGATGCGGGGCCGTACCGCTGTGGAGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGACCGCTAGCAGTTACTAAGTAAGATGCTTTAAGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGttt                                                                     <  2:248376/244‑1 (MQ=255)
                                                                                                                                                                   catttcgggtcgtaccgggCAGAGACGCGTCGCTACCGCTTTGGTGCGTCCCGCCTGCTGGCAAGCCTGTTGAAGGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGAGGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTATGAGACGTATCTTCCCGTAGGCCGGATAAGGCGTAGACGCATTCGGCAAAAGCCGCCTCCTACTCTTTTCCTCGGGAGGTTACCTTGTTtct                                                                   <  2:326521/235‑1 (MQ=255)
                                                                                                                                                                     gTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCg                                                                 >  1:335582/1‑251 (MQ=255)
                                                                                                                                                                         ggTGCTGACTGGGCAGATGAGCGTCACGACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCGGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCgg                                                                                    <  1:639027/228‑1 (MQ=255)
                                                                                                                                                                         ggTGCTGACTGGGCAGAGGCGCGTCACTCACGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCAAGTTTACGGTCAGCGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACTAGACGTCTCTCCCCACAGTACGGATACGGCGTAGACGACTCCGGCAACAGCCCCCTCATACTCTTTTCCTCgt                                                                                    >  2:639027/1‑227 (MQ=255)
                                                                                                                                                                            gCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGa                                                          <  2:127583/251‑1 (MQ=255)
                                                                                                                                                                            gCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGa                                                          <  1:612782/251‑1 (MQ=255)
                                                                                                                                                                               gccTGGGCAGATGCGAGGCACCACCGCTGTGGTGCGTCCATCCCGCTGGCAAGCCTGTTGAAAGCGTTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTCAGTCAGATGCTTTTCGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGCAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATTCTCTTTTCCTCGGGAGGTTACCTTGTTACGCGCACAAGAAat                                                       <  2:275555/249‑1 (MQ=255)
                                                                                                                                                                               gACTTGGCAGGTGCGCGTCACCACTGATTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGaaa                                                        <  2:281475/250‑1 (MQ=255)
                                                                                                                                                                                    ggCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGCTAAGGCGTAGACGCATCCGGCAAACGCCGCATCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTc                                                  >  2:534487/1‑251 (MQ=255)
                                                                                                                                                                                     gCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCg                                                 >  1:361341/1‑251 (MQ=255)
                                                                                                                                                                                                  ccTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGAt                                    <  1:72020/250‑1 (MQ=255)
                                                                                                                                                                                                  aCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGAt                                    <  1:5349/251‑1 (MQ=255)
                                                                                                                                                                                                  aCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGAt                                    <  1:44935/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGACAGCGCTGGGCCACGGTGACGGTAAGAGCGCAAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCAGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGGAATTATTCGTAAAAAACGTGATGGTc                                >  1:403130/1‑251 (MQ=255)
                                                                                                                                                                                                           ttGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATg                             <  1:443681/249‑1 (MQ=255)
                                                                                                                                                                                                           ttGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGAAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATgc                            >  1:658054/1‑250 (MQ=255)
                                                                                                                                                                                                             ggtgcTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGAAGCTTTACGCCTGATGCGCTGAGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTACTACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCAACAGAAATTATTCGTAAAACACGTGATGGTCaggaact                         >  2:649906/4‑245 (MQ=255)
                                                                                                                                                                                                              gtgcTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTg                        >  1:463604/3‑251 (MQ=255)
                                                                                                                                                                                                                           gctgGCAAGCCTGTTGAAAGCGCAGGGCCACGGTGACGGTACGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACACGAAATTATTAGTAAAAAACGTGATGGTAATGCACTGAGCGATgaagaa            >  2:6262/1‑250 (MQ=255)
                                                                                                                                                                                                                              ggCAAGCCTGTTGAAAGCGCTGGGCCTCGGTGACGTTAAGAGCGCCAGCAGCTACTAATTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAAtt         <  1:624689/250‑1 (MQ=255)
                                                                                                                                                                                                                               gCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCACAAGCCGCCTCATACTCTTTTCCTCGGGAg                                                                                 >  2:1408/1‑177 (MQ=255)
                                                                                                                                                                                                                                 aaGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGtt     >  2:620202/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cgacgccggAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATgcgctgcgct                                                                                                                                                                          >  1:46334/9‑85 (MQ=255)
                                                                                                                                                                                                                                  cgacgccggAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATgcgctgcgct                                                                                                                                                                          <  2:46334/77‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCt  >  1:166121/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
AACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCT  >  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: