Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,858,227 T→C 4.2% L59L (CTA→CTG)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,858,2270TC 4.2% 572.3 / 6.5 172L59L (CTA→CTG) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (87/75);  new base C (3/4);  total (91/81)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 172 total aligned reads displayed.

TTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATAACTGAT  >  NC_000913/3857943‑3858510
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                           
ttCGCGTGTAAAGCTTCGACGCTGTTTAGGCCATTGCGCAACAGGCTAATCGCTTTTTCCGTCAGTTCGCTAAGTTAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACAATCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTTGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGt                                                                                                                                                                                                                                                                                          <  1:371114/288‑1 (MQ=255)
     cTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGc                                                                                                                                                                                                                                                                                     >  1:1601149/1‑288 (MQ=255)
          aaTCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCGGTGc                                                                                                                                                                                                                                                                               >  1:1360755/1‑289 (MQ=255)
          aaTCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTgg                                                                                                                                                                                                                                                                               >  1:81587/1‑288 (MQ=255)
              ttCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGc                                                                                                                                                                                                                                                                           >  1:1507902/1‑289 (MQ=255)
               tCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCCTAGTGGAGATGTCTTTGCCTGGCGCGTCGTTGTTGGTCCACGGGGCGGACTGGTTCGC‑GTAGATACCGGTAAATAGTccctgttgcttataaa                                                                                                                                                                                                                                                                          >  2:1265201/1‑275 (MQ=255)
               tCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACGCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCg                                                                                                                                                                                                                                                                          >  2:406771/1‑289 (MQ=255)
                  aCGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTAc                                                                                                                                                                                                                                                                       >  1:429857/1‑289 (MQ=255)
                             gCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                                                                                                                                                                                 <  1:1186244/284‑1 (MQ=255)
                              ccATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGaa                                                                                                                                                                                                                                                            <  1:98969/288‑1 (MQ=255)
                                  tGCGCCATAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGCACAAACCGGTGAACAggg                                                                                                                                                                                                                                                       >  1:20232/1‑288 (MQ=255)
                                  tGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGgt                                                                                                                                                                                                                                                       >  1:964111/1‑289 (MQ=255)
                                       cacaGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATCCCGGTAACTAGTCCGGCCCGTGCAGGCGTCAAAACCGGTGAAAAGgt                                                                                                                                                                                                                                                       >  2:1186244/1‑284 (MQ=255)
                                       cacaGGCTAATCTATTTTTCCTTCCGTTCGCTAAGATAGTTCGCCCCCTCGAACCAGTAATCAGAGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGGCGAGATGCCAGTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCCGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAgg                                                                                                                                                                                                                                                  <  2:1440914/289‑1 (MQ=255)
                                        acaGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGc                                                                                                                                                                                                                                                 <  1:107645/289‑1 (MQ=255)
                                         caGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCCCGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGGGAACAGGTGTAGGCg                                                                                                                                                                                                                                                >  2:357445/1‑289 (MQ=255)
                                                           cGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGaa                                                                                                                                                                                                                                                            <  2:220925/259‑1 (MQ=255)
                                                             tCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGTTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTc                                                                                                                                                                                                                            <  2:712403/289‑1 (MQ=255)
                                                             tCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCtt                                                                                                                                                                                                                             >  2:770184/1‑288 (MQ=255)
                                                                      tAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCCCTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAg                                                                                                                                                                                                                   <  1:1111138/289‑1 (MQ=255)
                                                                        aggtatgtCGCCCCTTGGAACCAGTACTCGGCGTCCCCATCCGGCTGACCCTCGCCAGTGCCTAAATGGTCAGCACCATCGAGATGCCATTTGCCGATGTAACATGTGGGATAGCCGGCATCCTTAAAGTAGCGCCCCAGAGTGGAGAAGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGAc                                                                                                                                                                                                                 <  2:615792/282‑1 (MQ=255)
                                                                                  ccccATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTATTTACCGGTATCTA‑CGCTAACCAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACCGGTGTAgg                                                                                                                                                                                                                                                  >  1:316315/1‑246 (MQ=255)
                                                                                                 aaTCAACGTCCCAGTCCGGCGGACACACGCCAATGACGAAATAGTCATGACCGTCGAGATGCCATTTGCCGAGGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGGGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAg                                                                                                                                                                                           <  2:1296670/286‑1 (MQ=255)
                                                                                                     aGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTCGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGttt                                                                                                                                                                                     >  2:1146874/1‑288 (MQ=255)
                                                                                                     aGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGttt                                                                                                                                                                                     >  2:646411/1‑288 (MQ=255)
                                                                                                       cGTCACACTCCGGCGGACACGCGCCAGTGCAGAAATAGTAATGACCGTCGCGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCGCCATAGTGGAGATGTTTTTGACTGGCGCGACGTTGTTGGTCCACGGGACGGACTGGGTAGC‑GAAGATACCGGGACATAGTCCGGCCCCGGCAGGCGTACAAACAGGTGAAACGGTTTCGGCGGAATTCAAGAGAAACCCTTGCTCCGGCCGGCTACACATCATTTGCGTATTCGACGGCTTaa                                                                                                                                                                                   >  2:1501509/1‑287 (MQ=255)
                                                                                                        gTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTAcc                                                                                                                                                                                  <  1:1093158/288‑1 (MQ=255)
                                                                                                        gTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCg                                                                                                                                                                                 <  1:646411/289‑1 (MQ=255)
                                                                                                            cACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCtat                                                                                                                                                                             >  1:942088/1‑289 (MQ=255)
                                                                                                                 cggcggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCa                                                                                                                                                                         >  1:472858/1‑288 (MQ=255)
                                                                                                                 cggcggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGGGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTGTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCaa                                                                                                                                                                        <  2:795634/289‑1 (MQ=255)
                                                                                                                    cggACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAAcc                                                                                                                                                                      <  1:1170038/288‑1 (MQ=255)
                                                                                                                      gACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTACC‑GTAGATACCGGTAAATAGTCCGGCCCGTCCAGGCGTACAAACCGGTGAACAGGTTTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCCATATTTTGCTAATTCAGCGGTTTACCGCTATAGCAACCGa                                                                                                                                                                    >  2:1600590/1‑288 (MQ=255)
                                                                                                                               cAGTGCCGAAATCGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGGGGGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtggt                                                                                                                                                          <  1:674297/289‑1 (MQ=255)
                                                                                                                                   gCCGAAATATTCATGACCGGCGAGATGCCATTTGCCGATGGAACAGGTGTGATAGCCGGCATCCGTAAAGTAGCGCCCCATAGTGGAGATTGTTTTGCCTGGCGCGACGTTGTTGGGCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGc                                                                                                                                                       <  2:815366/288‑1 (MQ=255)
                                                                                                                                   gCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGc                                                                                                                                                       <  1:406771/288‑1 (MQ=255)
                                                                                                                                     cGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGc                                                                                                                                                                                                   <  1:90699/242‑1 (MQ=255)
                                                                                                                                     cGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGc                                                                                                                                                                                                   >  2:90699/1‑242 (MQ=255)
                                                                                                                                     cGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTg                                                                                                                                                    >  1:280871/1‑289 (MQ=255)
                                                                                                                                       aaaTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAgg                                                                                                                                                                                                                                                  >  1:1039577/1‑193 (MQ=255)
                                                                                                                                        aaTAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGt                                                                                                                                                 <  1:755986/289‑1 (MQ=255)
                                                                                                                                          tcgTCATGCCTGTCGGGATGCCATTTGCCGATGTAACGGTTAGGATAGCCGGCATCCTTAAAGTTGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTAt                                                                                                                                               <  2:1165768/287‑1 (MQ=255)
                                                                                                                                            gTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTCTCAATATTTTGCGTATTAAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGTTATCg                                                                                                                                             >  2:820433/1‑289 (MQ=255)
                                                                                                                                                     gTCGAGATGCCATGTGACAATTGAACAGGTGGGATAGCCGGCAGCCTTAAAGTAGGGCCCCATAGTGGAGATGTTTGTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACg                                                                                                                                    <  2:1394263/289‑1 (MQ=255)
                                                                                                                                                      tCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGa                                                                                                                                   >  2:470111/1‑289 (MQ=255)
                                                                                                                                                             gCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTCTCGGTCATGACGAACaacaa                                                                                                                            >  2:1335113/1‑285 (MQ=255)
                                                                                                                                                                     gcGATGTAACGGGTGTGATAGCCGGCATCCTTAAAGTGGCGCCCCATAGTGGAGATGTTTTTGCCTGGGGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTgggg                                                                                                                     <  1:583564/287‑1 (MQ=255)
                                                                                                                                                                      cGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCg                                                                                                                   <  1:460914/289‑1 (MQ=255)
                                                                                                                                                                         tGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTCTTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAAAAGAAAATTGGGGCggt                                                                                                                 >  1:1235859/1‑286 (MQ=255)
                                                                                                                                                                            aaCAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCAt                                                                                                             <  2:472858/289‑1 (MQ=255)
                                                                                                                                                                             aCAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGTGTATCGGTCATGCCGAACAGAAAATTGGGGCGTTTtata                                                                                                            >  2:1519576/1‑289 (MQ=255)
                                                                                                                                                                                  tgggATAGCCGGCATCCTTAACGTAGCGCCCCAGAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAgg                                                                                                          <  2:429857/283‑1 (MQ=255)
                                                                                                                                                                                        aGCCGGCATCCTTAAAGTAGCGCCCCATGGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                                                                 <  1:1335113/289‑1 (MQ=255)
                                                                                                                                                                                         gCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATGAAAGCGAATCCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACAATATTGGTGGCCTGGGTATCGGTCATGACGAACCGACAATTGGGGCGTTTAATAGGCTATCCTTc                                                                                                 >  2:310519/1‑288 (MQ=255)
                                                                                                                                                                                            ggCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGaa                                                                                                                                                                                                                                                            >  2:1615668/1‑130 (MQ=255)
                                                                                                                                                                                            ggCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGaa                                                                                                                                                                                                                                                            <  1:1615668/130‑1 (MQ=255)
                                                                                                                                                                                            ggCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCg                                                                                                                                             >  1:254701/1‑241 (MQ=255)
                                                                                                                                                                                              cATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTc                                                                                            <  2:808742/288‑1 (MQ=255)
                                                                                                                                                                                              cATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCg                                                                                           >  1:1015745/1‑289 (MQ=255)
                                                                                                                                                                                               aTCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCgg                                                                                          <  1:29689/289‑1 (MQ=255)
                                                                                                                                                                                                    taagttaccgcCCCTTGTTGGGGATGTTTTTGCTTGGCTCTACGTTGTGGGACCGCGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTa                                                                                     <  1:869301/281‑1 (MQ=255)
                                                                                                                                                                                                    gaaaGTAGCGCCCCAGAGTGGAGAGGTGTTGGCCTGGCGCGGCGTTGTTGGTCACCGGGCCGGGCTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTGGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTGGGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTa                                                                                     <  2:351330/288‑1 (MQ=255)
                                                                                                                                                                                                         acaggcCCCATAGTGGGGATGTTTTTTCCGTGGGGGACGTTGTTGGTCCACGGGACGGACTGGTTATC‑GGAGATACCGGTAAATACTCCGGACCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGt                                                                                <  2:887273/285‑1 (MQ=255)
                                                                                                                                                                                                          acagcCCCATAGTGGTGATGTTTTTGCCTGGGGTGACGTTGTTGGTCGTCGGGCCTGACTGGTTAGC‑GTAGATACCGGTAAATGGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCTTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGt                                                                                <  2:1296544/285‑1 (MQ=255)
                                                                                                                                                                                                           gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGCCGAACAGAAAATTTGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGtt                                                                               >  2:1050018/1‑288 (MQ=255)
                                                                                                                                                                                                           gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGtt                                                                               >  1:1476920/1‑288 (MQ=255)
                                                                                                                                                                                                               cccATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGc                                                                          <  1:80934/289‑1 (MQ=255)
                                                                                                                                                                                                                        gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTATTTACCGGTATCTA‑CGCTAACCAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAgg                                                                                                                                                                                                                                                  <  1:1511528/112‑1 (MQ=18)
                                                                                                                                                                                                                        gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTATTTACCGGTATCTA‑CGCTAACCAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAgg                                                                                                                                                                                                                                                  >  2:1511528/1‑112 (MQ=18)
                                                                                                                                                                                                                           aTGTTTTCGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGAAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTTATAGGCTATACTTCCCGTCGGGAGTACAGGTCACGGCAGGGCCGTGGCg                                                              >  1:23905/1‑289 (MQ=255)
                                                                                                                                                                                                                              tttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCa                                                                                                                                                                                                            >  2:236199/1‑144 (MQ=255)
                                                                                                                                                                                                                              tttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGAt                                                            >  1:223442/1‑288 (MQ=255)
                                                                                                                                                                                                                                     cTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAg                                                    <  2:1476920/289‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACACACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGGTGGGTCGTGACGATGCCGGCAGCTg                                                 >  1:597408/1‑289 (MQ=255)
                                                                                                                                                                                                                                               tgtacgcctgcACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCACGTGCCGGCGTACAAACCGGTGAACAGGTGTAGGAGGAATTAACGCGAATACCTTGCGCCGCCAGACTATCTATATTTTGCGTATTCAGCCGTTTACCGCTATCGCAACCGACCATATTCGTGGCCTGGGTATCGGTCATCCCGAACAGACAATTGGGGCGTTTTATAGGCTATCCTTCCCGCCGGGACTAAAGGTTACGGCTGGGTCGTGCCGCTCCCGCCGGCCGTCaaacc                                          >  2:269727/11‑287 (MQ=255)
                                                                                                                                                                                                                                                    tgGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAgg                                     <  1:352834/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                                                  <  1:955528/192‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                                                  >  2:955528/1‑192 (MQ=255)
                                                                                                                                                                                                                                                         agagcgacGGACTGGTTAAC‑GTAAATACCGGTTAATAGGCCGGGCCGTGAAGGGGTATAAACCGGTGAACAGGTTTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCATACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGAc                                <  2:1586460/282‑1 (MQ=255)
                                                                                                                                                                                                                                                           cGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTACAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACaa                              >  1:1545702/1‑289 (MQ=255)
                                                                                                                                                                                                                                                            tggCCGGACTGGTTAGC‑GTAGATACCGGTAAACAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAAt                             <  1:1144473/288‑1 (MQ=255)
                                                                                                                                                                                                                                                            gggCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAAt                             <  1:222838/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                cGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACt                         >  1:415091/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  gACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTAc                       >  1:1147387/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                        ttAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGGGGAAAGGTTACGGCTGGGCCGTGACGATGCCGCCGGCTGTCAAAAATCAGCTAGAAAATCCCTAcggctgt                 >  1:1586460/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                        ttAGC‑GGAGAACCCCGGAAAAAGTCCGGCGCGTGCAGGCGTACAAACCTGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACggcggc                 <  2:20232/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCaa               <  1:770184/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCa                                                                                                              <  1:995552/187‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATCGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAAGGCAAtc        >  1:1489600/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                  aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                                                                 >  2:311328/1‑200 (MQ=255)
                                                                                                                                                                                                                                                                                  aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCGGGTAGACAATCACTACGGCGGCAATGAAATa        <  2:942088/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATa        >  2:540635/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                  aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAAAAGGTGTCGGCGGAATTAAAGCGAATAACTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGAAACCGACCATATTGGTGGCCTGGGGATCGGTCATGAGGAACAGAAAATAGGAGCGTTACATAGGCTATCCTTCCAGCCGGGAGTAAAGGTTACGCTTCGGTCGTGACGATGCCGCCAGCTGTCAAAATTCTGCTATACCATCCCTCCGTCGGCCAAGAAATa        >  2:857946/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                       ggTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGGGTAGGAGGAATTAAAGCGAATACCTGCCGCGGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAACATATTGGTGGCCGGGGTATCGGTCATGACGAACAGAAAATTGGGGCGGTTAATAGGCTATCCTTCCCGTAGGGAGAAAAGGTTACGGCTGGGTCGTGACGATGCCGACAGCCGGCAAAAATCACGCAGACAATCACTAAGGCGGAAAACAAATAACTGa   >  2:480285/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                       ggTAAAAAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGTTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTTCCGTCGGGAGTAAAGGTTACGGCTGGGTCGCGACGATTCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTCCGGCGGCAATGaatacactga   >  2:835701/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                         tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATAACTGat  <  2:1561518/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATAACTGat  <  2:982743/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCAGGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATAACTGat  <  2:1015745/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tAAATAGTCCGACCCTTGCCGGCCTACAAACCGGTTAACAGGTTTAGGGGGAATTAAAGCGAATACCTTCCGCCGCCAAACTATCAATATTTTGCGTATTCTGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGTCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACGACGGCGGCAATGAAATAACTGat  <  2:1579532/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tAAAAAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTGGCGGATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATAACTGat  <  1:835701/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                           
TTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGACACTCGCCAGTGCCGAAATAGTCATGACCGTCGAGATGCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGC‑GTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACAATCACTACGGCGGCAATGAAATAACTGAT  >  NC_000913/3857943‑3858510

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