Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 499,060 T→G D305A (GAC→GCC)  aes ← acetylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913499,0600TG98.9% 332.7 / ‑5.3 93D305A (GAC→GCC) aesacetylesterase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (46/46);  total (47/46)
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.89e-01

GCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTAC  >  NC_000913/498813‑499305
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTcggcgg                                                                                                                                                                                                                                                    <  1:141766/251‑1 (MQ=255)
gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTcggcgg                                                                                                                                                                                                                                                    <  1:556353/251‑1 (MQ=255)
    ggAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGt                                                                                                                                                                                                                                                <  1:30571/251‑1 (MQ=255)
       gaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGgttt                                                                                                                                                                                                                                             >  1:416486/1‑251 (MQ=255)
             ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCCCGCCGGGACACATCGAAATTATGGGTAATCTTGTTTCCGCGTATCGCTAAAGCGGAGCGGGAAAGAAATGAGCGCCGTCGCTAAAAGCCTCGGCGGGGGGTTtcatcc                                                                                                                                                                                                                                        >  2:381254/1‑249 (MQ=255)
              gCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCatcatc                                                                                                                                                                                                                                      <  2:362482/251‑1 (MQ=255)
              gCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCatcatc                                                                                                                                                                                                                                      <  2:178034/251‑1 (MQ=255)
                      cATATTCAGGTGATGTGTCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagt                                                                                                                                                                                                                              <  2:71962‑M1/251‑5 (MQ=255)
                      cATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTTAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATTATCCGTGttgg                                                                                                                                                                                                                              >  2:345668‑M1/1‑247 (MQ=255)
                        tatTCAGGTGATATGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATTGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAAAGCCTCGGCGGCGGGTTTCAACATCCGTGtagtcg                                                                                                                                                                                                                            >  1:497199‑M1/1‑245 (MQ=255)
                            cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCGCGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGAGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtagg                                                                                                                                                                                                                        <  2:36583‑M1/251‑11 (MQ=255)
                            cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGCTTTCATCATCCGTGtagtcgtagg                                                                                                                                                                                                                        >  1:93207‑M1/1‑241 (MQ=255)
                             aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggc                                                                                                                                                                                                                       <  1:64337‑M1/251‑12 (MQ=255)
                                  aTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattag                                                                                                                                                                                                                  <  1:530732‑M1/251‑17 (MQ=255)
                                          gggCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgt                                                                                                                                                                                                          <  2:531413‑M1/251‑25 (MQ=255)
                                                 gctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCCCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgat                                                                                                                                                                                                   >  1:238947‑M1/1‑220 (MQ=255)
                                                   tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTCCGGCGCTTAATGCCACGCCGGGAAATATAGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGGTGAGCGGGAAAGAACTGAGCGGCGTCGGGAAGAGGCTCGGGGGGGGGTTTCATCATCCGTGgagtcgtagggattagacaattgtgcctgatgg                                                                                                                                                                                                 >  1:77298‑M1/1‑218 (MQ=255)
                                                    gcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATAACGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAAATGAGCGCCGTGGCGAAAAGCCCCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcatttgacatttgtgcctcatggc                                                                                                                                                                                                >  1:304913‑M1/1‑217 (MQ=255)
                                                       gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacg                                                                                                                                                                                             <  1:489552‑M1/251‑38 (MQ=255)
                                                           gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgctt                                                                                                                                                                                          <  1:693157‑M1/250‑41 (MQ=255)
                                                               ggAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgac                                                                                                                                                                                       <  2:674385‑M1/249‑44 (MQ=255)
                                                                  gaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTTCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgctgcgacgcttgacgcgtc                                                                                                                                                                                  >  1:147617‑M1/1‑203 (MQ=255)
                                                                   aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtag                                                                                                                                                                                                                         <  2:454075‑M1/211‑10 (MQ=255)
                                                                   aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtaggcgtag                                                                                                                                                                                                                         >  1:454075‑M1/1‑202 (MQ=255)
                                                                                   ccGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaacc                                                                                                                                                                 <  2:655466‑M1/251‑66 (MQ=255)
                                                                                      aGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggt                                                                                                                                                              <  2:223650‑M1/251‑69 (MQ=255)
                                                                                         aaaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatc                                                                                                                                                                             <  2:203563‑M1/233‑54 (MQ=255)
                                                                                         aaaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcggcgctttacgcgtcttatcagatcggaagaacacacc                                                                                                                                                           >  1:203563‑M1/1‑180 (MQ=255)
                                                                                               gTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAAAGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTTATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcctgacgcgtcttttcaggcctaccaccggggccccatcc                                                                                                                                                     >  1:585543‑M1/1‑174 (MQ=255)
                                                                                                        tCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaatt                                                                                                                                              <  1:192739‑M1/249‑85 (MQ=255)
                                                                                                        tCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaatt                                                                                                                                              <  1:186579‑M1/249‑85 (MQ=255)
                                                                                                         ccTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATTCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggt                                                                                                                                           <  2:77298‑M1/251‑88 (MQ=255)
                                                                                                         ccTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggt                                                                                                                                           <  1:381254‑M1/251‑88 (MQ=255)
                                                                                                           tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGGTAATCTTGTTGCCGCGCATCGCTAAAGCTGAGCGGTAAAGAACTGAGAGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAAGCGTGCAGCCGGCCTGGGGAGAGTTTGGACCCACAGGGCTACTGATTCGCCCCCTACGCCTGGGAAAGCAGACgg                                                                                                                                         >  1:613488/1‑251 (MQ=255)
                                                                                                                gCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttttgcctgatgcgacgcttgacgcgtcttatcaggcctacaacccgtgccgcctccggcaattggtgcacaa                                                                                                                                     >  1:541724‑M1/1‑157 (MQ=255)
                                                                                                                   ggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcc                                                                                                                                 <  1:595376‑M1/251‑98 (MQ=255)
                                                                                                                      gggAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggc                                                                                                                                                  <  2:31711‑M1/231‑81 (MQ=255)
                                                                                                                      gggAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggc                                                                                                                                                  >  1:31711‑M1/1‑151 (MQ=255)
                                                                                                                          aaaaaTATGAATAAATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGgagtcgtaggcgttagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcgattgggtcacaatgcctgaggcg                                                                                                                          >  2:667346‑M1/1‑147 (MQ=255)
                                                                                                                           aaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagt                                                                                                                                                                                                                              >  2:208095‑M1/1‑146 (MQ=255)
                                                                                                                           aaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagt                                                                                                                                                                                                                              <  1:208095‑M1/150‑5 (MQ=255)
                                                                                                                            naatATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacctttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                        >  2:327267‑M1/2‑145 (MQ=255)
                                                                                                                            aaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                        >  2:66076‑M1/1‑145 (MQ=255)
                                                                                                                             aaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatg                                                                                                                       >  1:5188‑M1/1‑144 (MQ=255)
                                                                                                                              atatGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgc                                                                                                                      >  1:602234‑M1/1‑143 (MQ=255)
                                                                                                                                atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                                    >  1:618867‑M1/1‑141 (MQ=255)
                                                                                                                                atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                                    <  1:345668‑M1/251‑111 (MQ=255)
                                                                                                                                      atatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgc                                                                                                               <  2:224175‑M1/250‑116 (MQ=255)
                                                                                                                                        atTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatc                                                                                                            >  1:571607‑M1/1‑133 (MQ=255)
                                                                                                                                          tCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgattgacgcatctt                                                                                                          >  2:104579‑M1/1‑131 (MQ=255)
                                                                                                                                           ccGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcctgacgcatctta                                                                                                         <  1:496087‑M1/251‑122 (MQ=255)
                                                                                                                                            cGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttat                                                                                                        >  1:30169‑M1/1‑129 (MQ=255)
                                                                                                                                             ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatc                                                                                                       >  1:516497‑M1/1‑128 (MQ=255)
                                                                                                                                               cgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatc                                                                                                       >  2:298521‑M1/1‑126 (MQ=255)
                                                                                                                                                gcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcagg                                                                                                    >  1:121578‑M1/1‑125 (MQ=255)
                                                                                                                                                gcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcagg                                                                                                    <  2:416486‑M1/251‑127 (MQ=255)
                                                                                                                                                 cTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggc                                                                                                   >  2:179149‑M1/1‑124 (MQ=255)
                                                                                                                                                 cTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggc                                                                                                   >  2:555970‑M1/1‑124 (MQ=255)
                                                                                                                                                 cTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgc                                                                                                                      <  2:424960‑M1/232‑109 (MQ=255)
                                                                                                                                                 cTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgc                                                                                                                      >  1:424960‑M1/1‑124 (MQ=255)
                                                                                                                                                  ttAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcc                                                                                                  >  2:409504‑M1/1‑123 (MQ=255)
                                                                                                                                                       gCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctaca                                                                                              >  1:699971‑M1/1‑118 (MQ=255)
                                                                                                                                                          aCGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgg                                                                                          <  1:649838‑M1/251‑137 (MQ=255)
                                                                                                                                                             ccGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggta                                                                                       <  2:497199‑M1/251‑140 (MQ=255)
                                                                                                                                                                   cATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggat                                                                                 <  1:513776‑M1/251‑146 (MQ=255)
                                                                                                                                                                    atatCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatc                                                                                <  2:640709‑M1/251‑147 (MQ=255)
                                                                                                                                                                          aaaTGATGGCTAGTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtagg                                                                          >  1:489919‑M1/1‑99 (MQ=255)
                                                                                                                                                                                    tAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataagg                                                                >  1:206363‑M1/1‑89 (MQ=255)
                                                                                                                                                                                          tgttgCCTCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgt                                                             <  2:699971‑M1/248‑166 (MQ=255)
                                                                                                                                                                                          tgttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgttt                                                           >  1:301323‑M1/1‑83 (MQ=255)
                                                                                                                                                                                          tgttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgt                                                             <  2:281319‑M1/248‑166 (MQ=255)
                                                                                                                                                                                           gttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttac                                                         <  1:179149‑M1/251‑170 (MQ=255)
                                                                                                                                                                                           gttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGCTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttac                                                         <  2:93207‑M1/251‑170 (MQ=255)
                                                                                                                                                                                               ccGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccg                                                     >  1:423314‑M1/1‑78 (MQ=255)
                                                                                                                                                                                               ccGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccg                                                     <  2:541724‑M1/251‑174 (MQ=255)
                                                                                                                                                                                                 gcgTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgc                                                    >  1:470946‑M1/1‑76 (MQ=255)
                                                                                                                                                                                                     aTCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGGGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtgtacgccgcatcc                                                <  2:602409‑M1/250‑179 (MQ=255)
                                                                                                                                                                                                     aTCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtagtccggataaggcgtttacgccgcatccg                                               <  1:298521‑M1/251‑180 (MQ=255)
                                                                                                                                                                                                       cGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggc                                             >  1:520143‑M1/1‑70 (MQ=255)
                                                                                                                                                                                                       cGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggc                                             <  2:549490‑M1/251‑182 (MQ=255)
                                                                                                                                                                                                           aaaGCTGAGCGGCAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggca                                            >  1:490519‑M1/1‑66 (MQ=255)
                                                                                                                                                                                                             aGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaat                                       >  1:52679‑M1/1‑64 (MQ=255)
                                                                                                                                                                                                             aGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaat                                       <  2:516497‑M1/251‑188 (MQ=255)
                                                                                                                                                                                                                 gAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatag                                     >  2:402435‑M1/1‑60 (MQ=255)
                                                                                                                                                                                                                    cGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcac                                <  2:301323‑M1/251‑195 (MQ=255)
                                                                                                                                                                                                                         aaGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagtt                           >  2:120969‑M1/1‑52 (MQ=255)
                                                                                                                                                                                                                                gAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccg                      <  1:360576‑M1/249‑205 (MQ=255)
                                                                                                                                                                                                                                         cgcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctg           <  1:409504‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                                                          gcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgc          <  1:327008‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                                          gcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaagagagcaccagttaaccgaacttactctgc          <  1:555970‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                                             aaGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggactacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcgataggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcc       >  1:172456‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                                agCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgccca     >  1:317753‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                                  ccTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaacggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcccaaat  >  2:459238‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                                  ccTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcc                                                                                                  <  1:73321‑M1/155‑129 (MQ=255)
                                                                                                                                                                                                                                                  ccTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcc                                                                                                  >  2:73321‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTAC  >  NC_000913/498813‑499305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: