Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,134 G→A G139G (GGC→GGT)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,1340GA95.4% 322.6 / ‑5.9 111G139G (GGC→GGT) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base G (2/0);  major base A (59/45);  minor base C (1/4);  total (62/49)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

ACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAA  >  NC_000913/65888‑66375
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
acgctgatatggtTTATTCAGCAGCGTTTGCTGCATATCCGTTTACTGCGGCGCTAACTGACGGCAGAATATCCCCTTATAATCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGTGAATGGACCAGGACGCCGGGCATTTTCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGt                                                                                                                                                                                                                                               <  1:474684/239‑1 (MQ=255)
aCGCAGATAGTGTTAATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGtt                                                                                                                                                                                                                                              <  2:201924/251‑1 (MQ=255)
  gCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTg                                                                                                                                                                                                                                             >  1:558663/1‑250 (MQ=255)
       tAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCt                                                                                                                                                                                                                                       >  1:657240/1‑251 (MQ=255)
         gtgtTTATTCAGCAGCGTTTGCTGCATATCAGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGCCCAGGGCGCCGGGCGTTTGCGCTGCATCGATACCCTGTTTTTCGATGGTTTCTACGATGACGTTACCTGTTTCCCACTCATATTCACCGTGGATTTCTgg                                                                                                                                                                                                                                     >  2:229970/1‑250 (MQ=255)
           gtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCg                                                                                                                                                                                                                                    <  1:529814/250‑1 (MQ=255)
             ttATTCAGCAGCGTTTGCTGCATATCCGGGAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACACGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGCAGTTAACGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCg                                                                                                                                                                                                                                 >  2:175925/1‑251 (MQ=255)
              tattcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGAAGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGACC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTAGAAGGTTTCTACGATGACGTTACCGGTTTACCACTCATATTCACCGTTGATTTCTGcggcgg                                                                                                                                                                                                                                >  1:982549/5‑251 (MQ=255)
               attcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCCCGCTGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGACTGGACAAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGACGGTTTCTACTATGACGTTCCCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCATCggg                                                                                                                                                                                                                               >  2:719547/4‑250 (MQ=255)
                 tcagcatCGTTTGCTGCATATCCGGGAACTGTGGCGCTAACTGACGGTATAATATCCCCATATAAGCGCCCTCTTCCAGCACGATGGCGTTATGCACCGCAGCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTTAGAATGGACCAGGACGCCGTGCATTTGCGCTGCATCGATGCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTTATTTcgtcggcggtc                                                                                                                                                                                                                              <  2:989645/249‑11 (MQ=255)
                  cagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAt                                                                                                                                                                                                                            >  1:197071/1‑251 (MQ=255)
                   agcagcGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAtt                                                                                                                                                                                                                           <  1:811938/251‑1 (MQ=255)
                   agcagcGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCCTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAtt                                                                                                                                                                                                                           <  1:480240/251‑1 (MQ=255)
                     cagcGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCACCATATAAGCGACCTCTTCAAGCAAGATGGCGTTATGCACCGCATCTTAGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCATGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAtttt                                                                                                                                                                                                                         >  1:938581/1‑251 (MQ=255)
                      atcGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATGTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTg                                                                                                                                                                                                                        <  1:764003/249‑1 (MQ=255)
                      agcGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCCGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTAATATTCACCGTTGATTTCTGCGTCGGTCATTTTg                                                                                                                                                                                                                        >  2:968409/1‑251 (MQ=255)
                         gttggctgcATATCCGGTAACTGCGGCGCTAACTGACGGAAGAATATCCCCATAGAAGCGTCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCgt                                                                                                                                                                                                                     <  2:411690/247‑2 (MQ=255)
                           ttGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGCTACCCTGTTTTTCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGt                                                                                                                                                                                                                   >  1:580587/1‑251 (MQ=255)
                                   tatCAGGTAACGGCGGCGCTAACTGCCGGCCGAATATCCCCATATAAGCGACCTCTTCCAGAACGATGGCGTTATGCACCGCATCTTCGACATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATAGACCAGGACGCCGGGCATTTGCGCAGCATCGATACCCTGTTTTTCGAAGGTTGCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAATTCTACTTCGGTCTTTTTGCGGGTGCGGGGa                                                                                                                                                                                                            >  2:677789/1‑250 (MQ=255)
                                    atCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAgggcat                                                                                                                                                                                                          >  1:182277/1‑248 (MQ=255)
                                    atCCGGTAACTGCGGCGCTAACTGACGGCAGAATATACCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCA‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGCTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAt                                                                                                                                                                                                          >  1:819531/1‑251 (MQ=255)
                                     tCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATCTAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGATACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCAGGTGCAGGGactg                                                                                                                                                                                                         >  2:48962/1‑248 (MQ=255)
                                       cGGTAACTGTGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGt                                                                                                                                                                                                       >  2:715684/1‑251 (MQ=255)
                                          tAACTACGGCGCTACCTGACGACAGAATATCCCAATATAAGCAACATCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCA‑ATGCACATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGCTACCCTGTTTTACGAAGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTTGTTTTCTGCGTCGGTCATTTTGCGGGCGCCGGGAATgggtcc                                                                                                                                                                                                    >  2:906768/1‑247 (MQ=255)
                                             cTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGt                                                                                                                                                                                                  >  1:600557/1‑250 (MQ=255)
                                                 ggcgCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGaa                                                                                                                                                                                              >  1:441045/1‑250 (MQ=255)
                                                  gcgcTAACTGACGGCAGAATAGCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAAt                                                                                                                                                                                            <  1:860572/251‑1 (MQ=255)
                                                       aaCTGACGGCAGAATATCCCCGTATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCg                                                                                                                                                                                       <  1:138280/251‑1 (MQ=255)
                                                       aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGCCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGACTGGTGCCGTAGAAATAGTCg                                                                                                                                                                                       >  2:715992/1‑251 (MQ=255)
                                                              ggCAGAATATCCCCATATAAGCGACCTCATCCAGCACGAGGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCAACGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTGTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtggg                                                                                                                                                                                <  2:818964/251‑1 (MQ=255)
                                                               gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggt                                                                                                                                                                               <  1:275549/251‑1 (MQ=255)
                                                                 aGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtgg                                                                                                                                                                             <  1:479404/251‑1 (MQ=255)
                                                                    atatCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCAACTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATAGCCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCCTATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCCGGGACTGGTGCCGTAGAAATAGACGGAgtgggtggatc                                                                                                                                                                          >  2:118082/1‑248 (MQ=255)
                                                                       tCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAACTGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTTCGAAGGTTTCTACGATGAAGTTACCGGTTTCCCCCTCATATTCACCGTTGCTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGACGGGTGCCGTAGAATTAGTCGGCGTGGGTGGTGCCgg                                                                                                                                                                       >  2:988727/1‑251 (MQ=255)
                                                                            atattaccGACCTCTTCCAGCACGGTGGCCTTATGCTCCGCATCCTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGGTAATTGTCCAGGGCGGCGGGCATTTGCGCTTGAACTATACCCTGTTTTTCCAAGGTTTCTACGATGACGTTACCTGTTTCCCAATTATATTCACCCTTTATTTCTGGGGCGGTCATTTTTCCGGTTCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCt                                                                                                                                                                  <  1:936847/244‑1 (MQ=255)
                                                                              ataAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATTGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTgg                                                                                                                                                                <  1:175819/251‑1 (MQ=255)
                                                                                aaGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTgg                                                                                                                                                                >  1:269596/1‑250 (MQ=255)
                                                                                 aGCGACCTCTTCAAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTGCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTTCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGact                                                                                                                                                             >  2:446990/1‑249 (MQ=255)
                                                                                  gCGACCTCTTCCAGCACGATAGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATc                                                                                                                                                            >  1:383630/1‑251 (MQ=255)
                                                                                   cGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGCCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGTTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGACTCg                                                                                                                                                           >  2:29751/1‑251 (MQ=255)
                                                                                    gACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCg                                                                                                                                                           >  1:685297/1‑250 (MQ=255)
                                                                                     aCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGAc                                                                                                                                                         >  1:698292/1‑251 (MQ=255)
                                                                                      ccTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACGCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                                                                                                                        <  2:61933/251‑1 (MQ=255)
                                                                                                   gATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTgcgc                                                                                                                                           >  1:70423/1‑251 (MQ=255)
                                                                                                    aTGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGTCGTGAGAATGGCCCAGGACGCCGGGCATTTGCGCTGCAGCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACGTTCCCGGTTTCCCACCCATATTCACCGTTGATTTCTGCGTCGGTCATTTGGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAAGCGACTGTCCCGCCTGCGcc                                                                                                                                          <  2:781459/251‑1 (MQ=255)
                                                                                                      gtcGTTATGCACCGCACCTTCGGCATTTTTGCCC‑AGGCAAATGGGCCGGGAGAATGGTCCAGGCCGCCGGGCAGTTGCGCTGCATCGATTCCCTGTTTTTCGAGGGTTTCTACGATGACGTTACCGGTGTCCCACGCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGTTTGCTGGAATCGACTGGCCCTCCTGCGCCCa                                                                                                                                        <  2:359155/249‑1 (MQ=255)
                                                                                                         gTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGCATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCCTATTCACCGTTGCTTTCTGCGTCGGTAATTTTGCGGCTGCAGGGAATGGTGCCGGACGAATAGTCGGCGTGGGAGGTGCCGGTTGCTGGCATCGGCAGGCCCGCCCGCGCCCAgct                                                                                                                                     >  2:841327/1‑249 (MQ=255)
                                                                                                         gTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccccgct                                                                                                                                     >  2:295607/1‑247 (MQ=255)
                                                                                                         gTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccc                                                                                                                                         >  1:392587/1‑247 (MQ=255)
                                                                                                            aTGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGACCGCCTGCGCCCAGAtag                                                                                                                                   <  1:61680/250‑3 (MQ=255)
                                                                                                                  cGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTTCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGGCTGGCCCGCCTGCGCCCAGATAgtggcgtg                                                                                                                            >  2:189431/1‑251 (MQ=255)
                                                                                                                  cGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtg                                                                                                                            >  1:753752/1‑251 (MQ=255)
                                                                                                                  cGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtg                                                                                                                            <  1:684701/251‑1 (MQ=255)
                                                                                                                  cGCATCTTCGGCATTTTTGCAC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggagtg                                                                                                                            >  2:437666/1‑251 (MQ=255)
                                                                                                                    cATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggc                                                                                                                          <  2:689619/251‑1 (MQ=255)
                                                                                                                          cGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGt                                                                                                                    >  1:529092/1‑251 (MQ=255)
                                                                                                                          cGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGt                                                                                                                    <  1:340564/251‑1 (MQ=255)
                                                                                                                           ggCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTAGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAAATAGTGGCGTGGGGCGAGTg                                                                                                                   >  2:838572/1‑251 (MQ=255)
                                                                                                                              aTTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAgggtgg                                                                                                                >  1:553181/1‑246 (MQ=255)
                                                                                                                                ttttGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGTCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgtg                                                                                                               <  2:761169/250‑5 (MQ=255)
                                                                                                                                  tttgcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGTgggtgtg                                                                                                             >  1:832586/5‑245 (MQ=255)
                                                                                                                                         aTGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATg                                                                                                      >  1:191033/1‑251 (MQ=255)
                                                                                                                                              aatgggCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccgcc                                                                                                 <  1:478406/248‑1 (MQ=255)
                                                                                                                                                 gagCCGTGAGAATGGACCAGAACGCCGGGCATTTGCGCTGCATAGATACCCTGTTTTTCGAAGGTTTATACGATGACGTTACCGGTTTCCCACTCATCTTCACCGTTGCTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGTGCCGGTTGCTGGAATAGACTGGCCCGCCTGCCCCCAGATAGTGGCGGGGCGCGAgggggggtgtgacatgccgccaat                                                                                              >  2:336813/3‑228 (MQ=255)
                                                                                                                                                  ggCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGGGTGTGTGCACAATGCCGCCAATg                                                                                             >  1:306856/1‑251 (MQ=255)
                                                                                                                                                     cGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGac                                                                                          <  1:934991/251‑2 (MQ=255)
                                                                                                                                                          gAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGa                                                                                      >  1:813487/1‑250 (MQ=255)
                                                                                                                                                             gggcccAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTGCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGc                                                                                  <  2:197071/247‑1 (MQ=255)
                                                                                                                                                                aCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTg                                                                               <  1:437666/251‑1 (MQ=255)
                                                                                                                                                                  caggacgcggggaagttGCGCTGCATCAATACCCTGTTTTTCGAAGGTATCTACGATGACGTTACCTGTTTCACACTCATATGCACCGTTGATTTCTGCGTCGGTCAGTTGGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGTTGGTGCCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCTTGGCGCGAGTGTGTGTGCACAATGCCTCCAATGGACGGGAACGCCTg                                                                               <  2:536217/234‑1 (MQ=255)
                                                                                                                                                                   aggaCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTAACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCAGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGAt                                                                             >  1:910173/4‑250 (MQ=255)
                                                                                                                                                                      gcgcagggCATTTGCGCTGCATCGATACCCTGGTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATGCCCCGTTGATTTCTGCGTCGGGCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATagag                                                                         <  2:235129/246‑1 (MQ=255)
                                                                                                                                                                       cGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGGGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTTGGGGGGCACAATGCCGCCAATGGACGGGACTGCCTGATATagc                                                                        >  1:208428/1‑251 (MQ=255)
                                                                                                                                                                          cgtgCATTTGCGCTGCATCGATACCCTGTTTTTCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGGGTGCAAAAGGCCGCCCATGGACGGGAATGcatttcagggcagt                                                                     >  2:609605/4‑238 (MQ=255)
                                                                                                                                                                              catttcggcttgtcgagtcccttGTTTTTCGAAGGTTTCTACGATGACGTTGCCGTTTTCCCACTCATATTCACCGTTGATTGCGGCGTCGGGCATGTTGCGGGTGCAGGGATTGGTGCCGTAGAAATAGTCGGCGTGGGTGGGTCCGGTGGCGGGAATCGATTGGCCCGCCGGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAAGGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGt                                                                 <  2:709655/229‑1 (MQ=255)
                                                                                                                                                                              aaTTTGCGCTGCATCGATACCCTGTTTTTCGAAGATTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGt                                                                 <  1:157359/250‑1 (MQ=255)
                                                                                                                                                                                         cATCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACGTTTCCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTGGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGTGGCGGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgta                                                      <  2:537055/251‑2 (MQ=255)
                                                                                                                                                                                              aTACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCggg                                                 >  1:135323/1‑250 (MQ=255)
                                                                                                                                                                                              aTACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCCATGGACGGGAATGCCTGATAGAGCCGTCGGTGAGTTGGCggatcggg                                                 >  1:67473/1‑244 (MQ=255)
                                                                                                                                                                                                aCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCggggg                                               >  1:315787/1‑248 (MQ=255)
                                                                                                                                                                                                  ccTGTTTTTGGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTAATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGCTGCTGGAATCGACGGGCCCGCCTGCGCCCAGGTAGAGGCGTGGGGCGAGGGGGGGTGCACAAAGACGCCAAAGGACGGGAGTGGATGATAGAGAAGGCGGGGAGTGGGCGGATCggagggg                                             >  2:514757/1‑249 (MQ=255)
                                                                                                                                                                                                     gTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAggg                                           <  1:709737/250‑1 (MQ=255)
                                                                                                                                                                                                        tttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCtt                                        <  1:760431/250‑1 (MQ=255)
                                                                                                                                                                                                        ttgcgaagttttCTACGATGACGTTTCCGGTTTCCCACTCATATTCACCGTTTATTTCTGCGTCGGTCATTTGTCGGGTGCGGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTTCCGGTTGCTGGAATCGACTGGCCCGTCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCtt                                        <  1:571396/241‑1 (MQ=255)
                                                                                                                                                                                                            gaaGGTTTCTACGATGACTTTACCGGTTTCCCACTCATAGTCACCGTTGATTTCTGCGTCGGTAATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCGGGAATCGACGGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCg                                   <  2:338692/250‑1 (MQ=255)
                                                                                                                                                                                                               ggTTTCTACGATGACGTTACCGGTTTCCCTCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCgcac                                <  1:661426/251‑4 (MQ=255)
                                                                                                                                                                                                                gTTTCTACGCTGTCGTGTCCGGTTTCCCACTCATATTCACCGTTGTTTGCTGCGGCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAAGAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGATTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCAcc                               <  2:187820/251‑1 (MQ=255)
                                                                                                                                                                                                                gTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGAGCACAATGGCGCCAATGGACGGGAATGCCTGATAGAGCCGGCGGAGAGTTGGAGTATCGGGGGAGGGCTTTTTCGCAcc                               >  1:515798/1‑251 (MQ=255)
                                                                                                                                                                                                                 tttCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCCGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCt                              >  1:45951/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttCTACGATAACGTTACCGGTTTCCCACTCATATTCACCGTTGTTTTCTGCGTCGGTAATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCATGCGCCCAGATAGCGGCGTGGGGCGAGTGGGTGAGCACAATGACGCCAATGGTGGGGGATGCCTGCTACAGCAGTCGGTGAGATGGCGTATATGAGGAGGGCTTTTTCgcccctt                             >  2:431300/1‑245 (MQ=255)
                                                                                                                                                                                                                   tCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGTTAGTGGCGTGGCGCGGGTGTGTGTGACCAATGCCGCCAATGGAAGGGAGTGCCTGGTCGGGCAGTCGGTGGGTTGGCGCATAGGAGGGGGGCGTTTTCgccccttc                            >  2:655068/1‑244 (MQ=255)
                                                                                                                                                                                                                    cTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGGTTTTTTGCACCTTc                            >  1:461869/1‑250 (MQ=255)
                                                                                                                                                                                                                         aTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGTTGGAATCGACTgg                                                                                                                                                      <  1:310424/123‑2 (MQ=255)
                                                                                                                                                                                                                                  ccgttgtCCCACTCATAATCCCCGTTGATTGCTGCGTCGGTCATTTGGGGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGTGCCGGTTGCGGGAATCGACTGGCCCGCCTGCGCCCAGAGAGTGGCGTGTCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGt             <  2:191033/245‑1 (MQ=255)
                                                                                                                                                                                                                                    gttttCCCACTCATATTCCCCTTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGttt           <  1:175925/249‑1 (MQ=255)
                                                                                                                                                                                                                                        tCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGa        <  1:427427/250‑1 (MQ=255)
                                                                                                                                                                                                                                        tCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGGGGGCTTTTTCGCAACTTCAACCACTTCACCGGTTTCGAt       >  1:627667/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cccccTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATg      <  1:646951/247‑1 (MQ=255)
                                                                                                                                                                                                                                          ccACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGGGGGGCGCGAGTGTGTGTGCCCAATGCCGCCCATGGGCGGGAAAGCCTGATAGGGCCGTCGGGGGGGTGGCGTATGGGAGGAGGGCGTTTTTGCCCCTTCAAACCCTTCACCggatttggtgc     >  2:534189/1‑242 (MQ=255)
                                                                                                                                                                                                                                               cATATTCACCGTGGAGTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTaa  <  1:127141/249‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAA  >  NC_000913/65888‑66375

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