Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,299 T→C S321S (TCT→TCC)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,2990TC98.9% 286.9 / NA 90S321S (TCT→TCC) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (51/38);  total (52/38)
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.65e-01

ATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑ACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGT  >  NC_000913/1061‑1541
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
aTGCTAGGTTGTTGAAGTCGATGTTCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGGTCTTCTCCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATTCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCTGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATg                                                                                                                                                                                                                                         <  2:12359/251‑1 (MQ=255)
       gttgttGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGAGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAggg                                                                                                                                                                                                                                   >  1:273300/1‑250 (MQ=255)
         tgttgAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGAt                                                                                                                                                                                                                                 <  1:227933/250‑1 (MQ=255)
          gttgAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATg                                                                                                                                                                                                                                <  1:69372/250‑1 (MQ=255)
             gAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATACCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGtt                                                                                                                                                                                                                             >  1:100579/1‑250 (MQ=255)
                   gATGTCCTACCAGGAAGCGATGGCGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATgg                                                                                                                                                                                                                      >  1:149596/1‑251 (MQ=255)
                          tACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGgcgcg                                                                                                                                                                                                               >  2:243403/1‑251 (MQ=255)
                           aCCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGgcgcgc                                                                                                                                                                                                              >  2:121099/1‑251 (MQ=255)
                              aGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGtg                                                                                                                                                                                                           >  1:519283/1‑250 (MQ=255)
                                  aGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTg                                                                                                                                                                                                       >  1:231809/1‑251 (MQ=255)
                                      attgAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTGCCAGATCCCTTGCCTGATTAAAAAGTCCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGTTGCCAGCCGTTATGAAGACGAATTACCGGTCAAGGGCAGTTCCAATCTGAATAATATGGCAAGGTTCTGCGTTTCCTGCCCGGGTATTAAAGGGATGGGTGGCATGGCGgcgcgctggtttgcagc                                                                                                                                                                                                   <  2:412174/248‑12 (MQ=255)
                                      aTGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCCTTCACCCC‑ATCGCCCAGTTCCAGATCACTTGACTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACACTCATTGGTGCCAGCCGTGATGAAGACGACTTACCGGCCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGGGTGTTTGCAGc                                                                                                                                                                                                   >  1:296005/1‑251 (MQ=255)
                                       tGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TACAGGTACGATCATTGGTGCCAGCCGTGATGAAGACGAATTACAGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCg                                                                                                                                                                                                  >  1:75706/1‑251 (MQ=255)
                                        ggAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATCACACCC‑ATCGCCCAGTTCCAAACCCCTTGCCTGATTAAAAATAACGCCAATCCTCAAGCTACC‑GGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAACGGCATTTCCAATCTGAAATATATGGCAATGTTCAGCGTTTCCGGCCCGGGGCTGAAAGGGATGGTTGGCATGGCGTCGCGCGGCTTTGCAGCg                                                                                                                                                                                                  >  1:581603/1‑250 (MQ=255)
                                             tttCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGgtgta                                                                                                                                                                                             >  2:483405/1‑249 (MQ=255)
                                             tttCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTc                                                                                                                                                                                             >  2:357042/1‑250 (MQ=255)
                                              ttCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCa                                                                                                                                                                                            <  1:357042/250‑1 (MQ=255)
                                                    cTTCGGCGCTAAAGTTATTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGACAGGGATGGTTGGCATGGCGGCGCGGGTGTTTGCAGCGATGTCACGCGccc                                                                                                                                                                                     >  1:515480/1‑251 (MQ=255)
                                                      tCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTTCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGccc                                                                                                                                                                                     <  1:582031/250‑1 (MQ=255)
                                                         gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCGTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTTGTGCCAGCCGTGATCAAGTCGAATTACCGGTCAAGGGCATTTCCTATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGTTGATGAAATGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCTCGt                                                                                                                                                                                   <  1:402833/249‑1 (MQ=255)
                                                           gcTAAAGTTCTTCCCCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTc                                                                                                                                                                              <  1:591227/251‑1 (MQ=255)
                                                                 gTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCggggg                                                                                                                                                                        >  1:172979/1‑247 (MQ=255)
                                                                      tcACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATGCCGGAAATCCTCAAGC‑TCCAGTTACGCTCATTGTTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGGATCCGGCCCGGGGATGAGAGGGATGGTTGGCATGTCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCt                                                                                                                                                                    <  1:359262/249‑1 (MQ=255)
                                                                      tcACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGTTTGCAGCGCAGTCACAAGCCCGTATTTCCG‑GGGTGCttg                                                                                                                                                                  >  1:53861/1‑249 (MQ=255)
                                                                      tcACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTAAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCt                                                                                                                                                                    <  1:352360/249‑1 (MQ=255)
                                                                       cACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGCTTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTGTGCGGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGAt                                                                                                                                                                 >  2:189294/1‑251 (MQ=255)
                                                                        ccccccGCACCATTACCCCC‑ATCGCCCAGTGCCAGATCCCGTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATGACCGGTCAAGGGCATTTCCAATCTGAATAATATGTCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGTTTGTCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGAtt                                                                                                                                                                <  2:123938/250‑1 (MQ=255)
                                                                             cGCACCATTACCCCC‑ATCGCCCCGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCG‑TGGTGCTGATTACGc                                                                                                                                                            >  1:161638/1‑250 (MQ=255)
                                                                                aCCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGGTTACGCAatc                                                                                                                                                        >  1:384469/1‑251 (MQ=255)
                                                                                aCCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAatc                                                                                                                                                        <  1:243403/251‑1 (MQ=255)
                                                                                 ccATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGTTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAatca                                                                                                                                                       <  1:23007/251‑1 (MQ=255)
                                                                                  cATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAatcat                                                                                                                                                      >  1:65395/1‑251 (MQ=255)
                                                                                   aTTACCCCC‑AGCGCCCAGTTCCAGATCCATTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAatcatc                                                                                                                                                     <  1:608666/251‑1 (MQ=255)
                                                                                     tACCCCC‑CTCGCCCAGTTCCAGATCCCTTGCATGATTACAAATACCGGAAATCCTCAAGC‑TCCAGGTCCGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGGTGTCCCGCGCCCGCATTTCCG‑TGGTGCTGATTTCGCaagc                                                                                                                                                        >  2:102503/1‑244 (MQ=255)
                                                                                     tACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAatc                                                                                                                                                        <  1:102503/246‑1 (MQ=255)
                                                                                       tcccacagcGCCCAGTTCCATGTCCCTCGCCTGACTACAAAGTCCGGACAACCTCAAGC‑TCCAGCTGCGCTCTTTGGTGCCAGCCGTGATTCAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATTATATAGGAATGTACAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCTGAGCGCGTGTTTGCAGCGATGTGACGCGCCCGTATGTCCG‑TGGTGCGGATTCCGCAATCATCtt                                                                                                                                                   <  2:337812/242‑1 (MQ=255)
                                                                                        cccccacagCCCAGTTCCAGATCCCTTGCCGGATTAAAAATCCCTGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTTCCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGTATGAAAGGGATGGTTGGTATGGCGGCGCACGGGGTTGCAGCGATGTCTCGCGCCCGTTTTTGCGATGG‑GCTGATTACGctgtcagcttcc                                                                                                                                                 <  2:603549/243‑12 (MQ=255)
                                                                                                  ccAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCGTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCTGTCAAGGTCATTTCCAATCTGAATAATATGGCAATGTGCAGCTTTTCCGGCCCGGGGATTAAAGGGATGGTAGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGGATTCATCTTCCTAATACatta                                                                                                                                       <  2:35527/251‑4 (MQ=255)
                                                                                                  ccAGTACCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCGCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTGCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCTCGCGCCCGTATATCCG‑TGGGGCTGATTACGCAATCATCTTCCGAATACAgta                                                                                                                                       <  2:60219/251‑3 (MQ=255)
                                                                                                         cAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAACACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGttt                                                                                                                                >  2:414465/1‑251 (MQ=255)
                                                                                                          aGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCCGCGCTGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTc                                                                                                                               >  1:38696/1‑251 (MQ=255)
                                                                                                           gATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCt                                                                                                                              <  2:161638/251‑1 (MQ=255)
                                                                                                                 ttGCCTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTc                                                                                                                        >  1:561776/1‑251 (MQ=255)
                                                                                                                   gCCTGATTAAAAATACCGGAACTCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGGATACAGTATCAGTTTCTGCGTTccg                                                                                                                      >  2:110636/1‑250 (MQ=255)
                                                                                                                     cTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCACGCGCCCGGATTTCCG‑TGGTGCTGGTTACGCAATCATCTTCCGCATACAGTATCAGTTTCTGCGTTccgc                                                                                                                     >  2:472822/1‑248 (MQ=255)
                                                                                                                     cTGATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTccgc                                                                                                                     >  2:175262/1‑248 (MQ=255)
                                                                                                                       gATTAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGTCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCaaa                                                                                                                  <  1:323895/251‑1 (MQ=255)
                                                                                                                          tAAAAATACCGGAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGc                                                                                                                <  1:110636/250‑1 (MQ=255)
                                                                                                                           aaaaaTACCGGAAATCCCCAAGC‑TCCATGTTCGCGCATTGGTGCCAGCCGGGATGAAGACGAAGTACCGGGCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCGGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTGCGTTCCGCAAAGCg                                                                                                               <  2:381568/250‑1 (MQ=255)
                                                                                                                                 aCCGGAAATTCTCAAGC‑TCCAGGTACGCGCATTGGTGCCAGCCGTGAAGAAGACGCATTGCCGGTCAAGGGCATTTCCAATCTGAACACTAGGGCAATGTTCAGCGTTTCCGGCCCGGGGAAGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTTCG‑TGGTGCTGATTACGCAATCATCTTCCGAATCCAGTATCAGTTTCTGCGTCCCGCAAAGCGACtgtgt                                                                                                        <  2:533448/251‑1 (MQ=255)
                                                                                                                                   cgtaaaaccTCAAGC‑TCCAGGTACGCGCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAAGACAGTATCATTTTCTGCGTTCAGCAAAGCGACTGTGTGc                                                                                                      <  2:87690/244‑1 (MQ=255)
                                                                                                                                    ggAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCg                                                                                                     >  1:434441/1‑251 (MQ=255)
                                                                                                                                     gAAATCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCACTCTGAATAATATGGCACTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGAGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCgg                                                                                                    >  1:377384/1‑250 (MQ=255)
                                                                                                                                      aaaTCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAg                                                                                                   <  1:414465/251‑1 (MQ=255)
                                                                                                                                       aaTCCTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGTTGAAAGGGATGGTTGGCATGGCGCCGCGCGTGTTTGCAGAGATGTCACGCGCCCGTATTTCAG‑TGGTGCTGATTACGCAATCATCTTCCGAATCCAGTATCAGTTTCTGCGTTCCGCAAAGCGTCTGTGTGCGAGc                                                                                                  >  1:382787/1‑251 (MQ=255)
                                                                                                                                          ccTCAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTTATGAAGACGAAGTACCGGTCAAGGTCATTGCCAATCTGAATAAGAGGGCAATGTTCAGCTTTTCCGGCCCGGGGATGAAAGGGGTGGTTGGCATGGCGGCTCGCGTGTTTGCAGCGATGTCACGCGCCTGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACATTATCAGTTGCTGCGTTCCGCAAAGCGACTGTGTTCGAGCTGa                                                                                               <  2:296005/251‑1 (MQ=255)
                                                                                                                                            ccAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAAc                                                                                             <  2:129594/250‑1 (MQ=255)
                                                                                                                                             cAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCCGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACg                                                                                            >  2:132323/1‑251 (MQ=255)
                                                                                                                                             cAAGC‑TCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACg                                                                                            >  1:557567/1‑251 (MQ=255)
                                                                                                                                                    ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCaa                                                                                       <  1:253078/250‑1 (MQ=255)
                                                                                                                                                    ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGAGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                      >  1:396849/1‑251 (MQ=255)
                                                                                                                                                     cAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg                                                                                     >  1:592189/1‑251 (MQ=255)
                                                                                                                                                        gTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCACGGGCATTTCCAATCTGAATAATATGGCAATGTTCCGCGTTTCCGGCCTGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCCCGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGGATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGc                                                                                    >  1:619097/1‑249 (MQ=255)
                                                                                                                                                         gaCGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGGCAAGGGCATGTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAgg                                                                                 <  2:515480/250‑1 (MQ=255)
                                                                                                                                                           cGCTCATTGGAGCCAGCCGTGATGAAGACGAATTACCGGACAAGGGCATTTCCAATCTGGATAATATGGAAATGTTCAGCGTTTCCGGCACGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGAGTGTTTGCCGCGCTGTCACGGGCCCGTATTTCCG‑GGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCACAACGACTTGGCGCGAGCGGAACGGGCAATGCAggg                                                                                >  2:621644/1‑249 (MQ=255)
                                                                                                                                                                aTTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGtt                                                                          <  1:132323/251‑1 (MQ=255)
                                                                                                                                                                  tGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTc                                                                         <  1:542522/250‑1 (MQ=255)
                                                                                                                                                                   ggTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCt                                                                        <  1:189294/250‑1 (MQ=255)
                                                                                                                                                                      gacAGCCGTGATGAAGACTAATTACCTGTCAAGGGCAGTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAACGGGCAATGCAGGAATAGTTCTAcc                                                                     <  2:377384/248‑1 (MQ=255)
                                                                                                                                                                       ccAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATAGGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTACCTg                                                                   <  2:400791/251‑1 (MQ=255)
                                                                                                                                                                         aGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGAGGTTTTGTAGCGATGGCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTACGAATACCGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGCAGAGTTCTACCttgg                                                                 >  1:573853/1‑248 (MQ=255)
                                                                                                                                                                              tgatgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGGGTTTGCAGCGCTGTCACGCGCCCGTATTTCCG‑TGGGGCTGATTACGCAATCCTCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTg                                                             >  2:185186/1‑250 (MQ=255)
                                                                                                                                                                                   aaGACGAATTACCGGTCAAGGGCAGTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTACGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAaagaa                                                       <  2:401480/251‑1 (MQ=255)
                                                                                                                                                                                        gAATTACCGGTCAAGGGCATTTCAAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCCGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGACCTGGAAGAAGGCtt                                                  >  2:35643/1‑251 (MQ=255)
                                                                                                                                                                                         aaTTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCCATGCAGGAAGAGTTCTACCTGGGACTGAAAGAAGGCtt                                                  >  2:215259/1‑250 (MQ=255)
                                                                                                                                                                                           ttACCGGTCAAGGGCATTTCCACTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACt                                               >  1:496727/1‑251 (MQ=255)
                                                                                                                                                                                                ggTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGc                                          <  2:519283/251‑1 (MQ=255)
                                                                                                                                                                                                ggTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGc                                          >  1:311910/1‑251 (MQ=255)
                                                                                                                                                                                                    aaGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCt                                      >  2:234876/1‑251 (MQ=255)
                                                                                                                                                                                                      gggCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGGGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTGACTGGAGCCGCTgg                                    >  2:482308/1‑251 (MQ=255)
                                                                                                                                                                                                          aTTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGcggg                                >  1:271458/1‑248 (MQ=255)
                                                                                                                                                                                                                aaTCTGAATAATATGGCAATGTTCAGCGTTTCTGGCCCGGGGCTGACAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCCATCCTCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg                                                                                     >  2:112717/1‑192 (MQ=255)
                                                                                                                                                                                                                aaTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATg                                                                                     <  1:112717/192‑1 (MQ=255)
                                                                                                                                                                                                                   cTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGTCTTACTGGAGCCGCTGGCGGTGACGGGACg                       >  2:224192/1‑251 (MQ=255)
                                                                                                                                                                                                                       ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACGGAACGGCTg                   >  1:267049/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGGTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGCGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGCGGAAAGGCGGGCCCTTAtct         >  2:294235/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCg       >  1:19880/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGAATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGCAACCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGCAAGAAGGCTTAATGGAGACGCTGGCGGCGACCGGACGGCTGGcccatttctcggggg   >  1:572932/1‑238 (MQ=255)
                                                                                                                                                                                                                                          cGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGAGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCggtggt  >  1:576449/1‑249 (MQ=255)
                                                                                                                                                                                                                                          cGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGAGCGCGTGTTTGCAGCGATGTCACGCGCCCGTCTTTCCG‑TGGCGCTGATTACGCCATCATCTTCCGAATCCAGTATCCGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTCTCTCgggggt  >  1:562241/1‑249 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
ATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCC‑ATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGC‑ACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCG‑TGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGT  >  NC_000913/1061‑1541

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