Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,585,613 C→T 3.1% E384K (GAA→AAA)  hsdR ← endonuclease R Type I restriction enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,6130CT 3.1% 513.0 / 2.5 149E384K (GAA→AAA) hsdRendonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base C (59/84);  new base T (3/2);  total (63/86)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 151 total aligned reads displayed.

ACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGT  >  NC_000913/4585326‑4585892
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                       
aCGGGTTTTACCGGTACCGGTCGCCATCGCCATCAGGATCTCTTGTAGCCCCTTGACGCTTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCAGTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTcg                                                                                                                                                                                                                                                                                        <  2:544243/289‑1 (MQ=255)
         accggtaccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCAAGCTCGCTCATGCCCGGGTTATCCGCAAACCACTGATTCTGTCTTTGCGGTTCGCTGCCGCGAATTTACAGCCGCTATTCCGGGCGGTTCCACTCGGGTACGGCTTTCGACATATTGCGGGTATCACCCCCGTCGCCATACTAGATGCCGCCTTTGGTTCTCATTGTTGCGCGGTATCTGCGcctggtg                                                                                                                                                                                                                                                                               >  2:405223/1‑284 (MQ=255)
               accggtCGCCATCGCCAGCAGGATCGCTGGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAg                                                                                                                                                                                                                                                                         <  1:250402/289‑1 (MQ=255)
                ccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGt                                                                                                                                                                                                                                                                        <  1:1282688/289‑1 (MQ=255)
                        catcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAg                                                                                                                                                                                                                                                                >  1:139660/1‑289 (MQ=255)
                         atcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg                                                                                                                                                                                                                                                               <  1:1221971/289‑1 (MQ=255)
                         atcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg                                                                                                                                                                                                                                                               <  2:422129/289‑1 (MQ=255)
                         atcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg                                                                                                                                                                                                                                                               <  1:488206/289‑1 (MQ=255)
                         atcccgacccggttctcTTTTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg                                                                                                                                                                                                                                                               <  2:1310054/276‑1 (MQ=255)
                                     ttctaTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAAc                                                                                                                                                                                                                                                    <  1:1313164/287‑1 (MQ=255)
                                          ttgttgCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTa                                                                                                                                                                                                                                                       >  1:180436/1‑280 (MQ=255)
                                          ttgttgCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTa                                                                                                                                                                                                                                                       <  2:180436/280‑1 (MQ=255)
                                           tgttgCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCgttgtt                                                                                                                                                                                                                                             >  1:1102256/1‑289 (MQ=255)
                                           tgttgCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCgttgtt                                                                                                                                                                                                                                             <  2:741285/289‑1 (MQ=255)
                                               gCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGcc                                                                                                                                                                                                                                         >  1:1087096/1‑289 (MQ=255)
                                                    ttGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTgg                                                                                                                                                                                                                                    >  2:785788/1‑289 (MQ=255)
                                                       aCGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGgtgt                                                                                                                                                                                                                                 >  1:416528/1‑289 (MQ=255)
                                                                ttttCCCCCGCGCGGCCGGCATCTTACTGATAATAACGCCGGCCGAGCTCGCTCATGCCAGGGTTATCGGCCAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCGCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCTATACCAGATGCCGCTGTTGGTGTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                                                                                                                                                                                        <  1:419137/289‑1 (MQ=255)
                                                                ttttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCCGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                                                                                                                                                                                        <  1:453582/289‑1 (MQ=255)
                                                                ttttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                                                                                                                                                                                        <  1:127711/289‑1 (MQ=255)
                                                                ttttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCACGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTGTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                                                                                                                                                                                        <  2:3040/289‑1 (MQ=255)
                                                                  ttCAACCGCGCGGACGGCATCTTCCTGCTAATAACGCAGGCCTAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTGTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTg                                                                                                                                                                                                                      <  1:632119/289‑1 (MQ=255)
                                                                        cgcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCa                                                                                                                                                                                                                <  1:668423/289‑1 (MQ=255)
                                                                         gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCTTGACAGCTGTTTTCAt                                                                                                                                                                                                               >  1:1129873/1‑289 (MQ=255)
                                                                                 ggCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCAAACTCTgg                                                                                                                                                                                                        <  1:8920/288‑1 (MQ=255)
                                                                                   cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCa                                                                                                                                                                                                      <  1:826961/288‑1 (MQ=255)
                                                                                   cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCa                                                                                                                                                                                                      <  1:107144/288‑1 (MQ=255)
                                                                                    aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACt                                                                                                                                                                                                    <  1:340596/289‑1 (MQ=255)
                                                                                          cTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTTGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCGTCa                                                                                                                                                                                              <  2:202953/289‑1 (MQ=255)
                                                                                               aataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCAc                                                                                                                                                                                         <  1:296645/289‑1 (MQ=255)
                                                                                                  aaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTc                                                                                                                                                                                      <  2:266901/289‑1 (MQ=255)
                                                                                                    cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAACCCACTGATTCTGGTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                    <  2:325834/289‑1 (MQ=255)
                                                                                                    cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                    <  1:1002173/289‑1 (MQ=255)
                                                                                                    cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                    <  2:1102256/289‑1 (MQ=255)
                                                                                                    cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                    <  2:252161/289‑1 (MQ=255)
                                                                                                     gccggCCCTGCTCGCTCATGCCAGGGTTATCGGGAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCTAGCATTTCCAGCAGCTCTTACGGGCGGTGCCACTCGGGTCAGGTTTTCGACATATTGCGGGTATCGCGCACGTCGCGATACCAGATGCCGCGTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAGACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                    <  2:324388/285‑1 (MQ=255)
                                                                                                      cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACCCTGGCACTGCTTCCTGCACTtcctcc                                                                                                                                                                                  >  2:709297/1‑289 (MQ=255)
                                                                                                       gggCCCAGCTCGCTCATGCCAGGGTTATCGGCAATCCACTGATTCTTTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCTCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                                                                                                 <  2:652828/288‑1 (MQ=255)
                                                                                                         gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCCGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGt                                                                                                                                                                               >  1:362898/1‑289 (MQ=255)
                                                                                                          cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                                                                                                 <  2:1087096/286‑1 (MQ=255)
                                                                                                              gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagta                                                                                                                                                                          >  1:834825/1‑289 (MQ=255)
                                                                                                               ctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCCTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCCAGTAGCCGCAGGGGATTTTAAACCGTTGGTTGCCGCTGGTGTCCTGCCCGCTGGTTTCATACTTTTGCACTGCCTCATTCACTTCCTTCGGTGagtag                                                                                                                                                                         >  2:565351/1‑289 (MQ=255)
                                                                                                                 cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCGGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTg                                                                                                                                                                       <  1:919566/289‑1 (MQ=255)
                                                                                                                 cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCACGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTg                                                                                                                                                                       <  1:1194990/289‑1 (MQ=255)
                                                                                                                   ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGc                                                                                                                                                                      >  2:211239/1‑288 (MQ=255)
                                                                                                                   ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGc                                                                                                                                                                      >  2:397423/1‑288 (MQ=255)
                                                                                                                          cAGGGTTATCGGCAAACCACGTATTCGGTGTTGGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGTTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaagca                                                                                                                                                               <  2:1050080/288‑1 (MQ=255)
                                                                                                                             ggTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTa                                                                                                                                                                                                                                                       >  1:1112496/1‑197 (MQ=255)
                                                                                                                                  tCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTcc                                                                                                                                                       >  2:974224/1‑288 (MQ=255)
                                                                                                                                       aaaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAgg                                                                                                                                                 >  1:1014196/1‑289 (MQ=255)
                                                                                                                                       aaaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAgg                                                                                                                                                 <  2:814703/289‑1 (MQ=255)
                                                                                                                                       aaaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAg                                                                                                                                                  <  1:785788/288‑1 (MQ=255)
                                                                                                                                        aaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCCCTGCTTCATGCACTTCATCCGGTGAGTAGTGCCCAAGCAAGGTTTCCCGCAGGa                                                                                                                                                >  2:641744/1‑289 (MQ=255)
                                                                                                                                         aCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCt                                                                                                                                                                                                 <  1:151692/239‑1 (MQ=255)
                                                                                                                                                  tCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATc                                                                                                                                      <  1:905511/289‑1 (MQ=255)
                                                                                                                                                    tGTTTTTGCGGTTCGCTGCCGAGCATTTCCGGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGa                                                                                                                                    <  1:440924/289‑1 (MQ=255)
                                                                                                                                                        tttGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGc                                                                                                                                                                                           <  1:187047/230‑1 (MQ=255)
                                                                                                                                                        tttGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGc                                                                                                                                                                                           >  2:187047/1‑230 (MQ=255)
                                                                                                                                                          ttcGGTTCGCTGCCCAGCATTTCCAGCAGCTCGTCCGGGCGGTGCCTCTCGGGTAAGGCATACGACATATTGCGGGTATCACGCACGTCGCGATACCAGTTGTCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCTTGTAGCAGAAGTGGATTTTAATCCGTTGTTTGCCGCTGGTGTCCTTCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAtt                                                                                                                              <  1:191154/287‑1 (MQ=255)
                                                                                                                                                          tGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTTTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAACCAGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCCGCCACTGCTTCATGCACTTCATCAGGGGAGTAGTGTCCCAGCAAGGTTTCCCGAAGGAAGCCATTATCGAAACAtt                                                                                                                              >  2:1232953/1‑289 (MQ=255)
                                                                                                                                                          tGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTGAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAtt                                                                                                                              <  2:139660/289‑1 (MQ=255)
                                                                                                                                                                tCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGGAGGGGATTTTAAACCGTTGGTTGCCGCGGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCCTCCGGTGAGGAGTTCTCAAGCAAGGTTTCCCGGCGGCAGCCATTATCGAAACATTTAATGt                                                                                                                        >  2:696437/1‑289 (MQ=255)
                                                                                                                                                                tCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                                                                                                 >  2:342487/1‑232 (MQ=255)
                                                                                                                                                                  gCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAg                                                                                                                      <  1:748914/289‑1 (MQ=255)
                                                                                                                                                                       cGAGCATTTCCAGCAACTCTTCCGGGCTGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATCCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTCAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                                                                                                  <  1:970085/288‑1 (MQ=255)
                                                                                                                                                                           cATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGa                                                                                                              <  2:964023/288‑1 (MQ=255)
                                                                                                                                                                                        tctTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTgg                                                                                                >  1:522053/1‑289 (MQ=255)
                                                                                                                                                                                            ccGGGCGGTGCCAATCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTTGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCggg                                                                                            <  1:575873/289‑1 (MQ=255)
                                                                                                                                                                                                  ggTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGc                                                                                                    >  1:343494/1‑275 (MQ=255)
                                                                                                                                                                                                  ggTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGc                                                                                                    <  2:343494/275‑1 (MQ=255)
                                                                                                                                                                                                    tGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAACCATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGACCGTCGa                                                                                    >  1:591436/1‑289 (MQ=255)
                                                                                                                                                                                                    tGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa                                                                                    >  2:122000/1‑289 (MQ=255)
                                                                                                                                                                                                    tGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa                                                                                    >  1:617264/1‑289 (MQ=255)
                                                                                                                                                                                                     gCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa                                                                                    <  1:696437/288‑1 (MQ=255)
                                                                                                                                                                                                       cACTCGGGTAAGGCTTTCGTTATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGata                                                                                  <  1:1232953/288‑1 (MQ=255)
                                                                                                                                                                                                                aaGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGtt                                                                        <  1:687155/289‑1 (MQ=255)
                                                                                                                                                                                                                        cGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCctct                                                                >  2:1119222/1‑289 (MQ=255)
                                                                                                                                                                                                                              atTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccacc                                                          >  1:179058/1‑289 (MQ=255)
                                                                                                                                                                                                                              atTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccac                                                           >  2:421537/1‑288 (MQ=255)
                                                                                                                                                                                                                              atTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCCCTaccac                                                           >  1:1202625/1‑288 (MQ=255)
                                                                                                                                                                                                                                tGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACcgc                                                        >  1:150297/1‑289 (MQ=255)
                                                                                                                                                                                                                                    ggTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGGGAGTAGTGCTCAAGCAAGGTTTCCCCCTGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTTCCACCCCGATg                                                    >  2:939993/1‑289 (MQ=255)
                                                                                                                                                                                                                                       aTCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGATGGTGTCCTGCCAGCTGGGTTCATAATCTGGCACTGCTTCATGCACGTCATCCGGTGCGCAGTGCTCAAGCAAGGTTTCCCTCAGGAAGCCATTATCGAAACATTTCCCGTAGCGATAAGACTCATTTAGCCTGGCGGGAACGTCGACATTGTTCAGTTTCTCCTCTCCCACAGCGTTGGGt                                                 >  2:898104/1‑289 (MQ=255)
                                                                                                                                                                                                                                        tCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGtt                                                >  1:785172/1‑289 (MQ=255)
                                                                                                                                                                                                                                           cgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTAACACCGCGATGGGTTTGa                                             >  1:739909/1‑289 (MQ=255)
                                                                                                                                                                                                                                              acgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAgg                                           <  1:709297/288‑1 (MQ=255)
                                                                                                                                                                                                                                               cgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGcc                                         >  1:1128970/1‑289 (MQ=255)
                                                                                                                                                                                                                                                gTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCg                                        >  1:313974/1‑289 (MQ=255)
                                                                                                                                                                                                                                                gTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCg                                        >  2:105401/1‑289 (MQ=255)
                                                                                                                                                                                                                                                  cgcgATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                                                                                                                                                                                        >  1:736225/1‑111 (MQ=255)
                                                                                                                                                                                                                                                    cgctACCGGTTGCCGCTTTTGGGCTTCATGTTTTGGCGTTTTTCGCGCCCGTTGGTACGGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGTTTTCCCGCAAGAAGCCATTATCGAAACATTTACTTTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACa                                     <  2:1129873/285‑1 (MQ=255)
                                                                                                                                                                                                                                                     gATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACaaa                                   <  1:32589/289‑1 (MQ=255)
                                                                                                                                                                                                                                                       tACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGGGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAAc                                  <  1:335319/288‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                                                                                                  >  2:506205/1‑200 (MQ=255)
                                                                                                                                                                                                                                                               gCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                                                                                                  <  1:506205/200‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                                                                                                  >  2:533637/1‑196 (MQ=255)
                                                                                                                                                                                                                                                                   cTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                                                                                                  <  1:533637/196‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      cgcgGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa                                                                                    <  1:237690/207‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        cgGTATTTGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCg                                        >  1:27734/1‑249 (MQ=255)
                                                                                                                                                                                                                                                                                        cgGTATTTGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCg                                        <  2:27734/249‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                       
ACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGT  >  NC_000913/4585326‑4585892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: