Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 581,640 A→T 5.0% pseudogene (287/309 nt) aaaD → pseudogene, DLP12 prophage, tail fiber assembly protein family,Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913581,6400AT 5.0% 283.3 / 3.1 79pseudogene (287/309 nt)aaaDpseudogene, DLP12 prophage, tail fiber assembly protein family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (44/31);  new base T (2/2);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGATGCAGG  >  NC_000913/581364‑581927
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                               
cGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa                                                                                                                                                                                                                                                                                     <  1:144639/289‑1 (MQ=255)
 gAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCTTAAGCAGCGCAac                                                                                                                                                                                                                                                                                    >  1:86457/1‑289 (MQ=255)
 gAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAac                                                                                                                                                                                                                                                                                    >  2:455938/1‑289 (MQ=255)
                 aGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTg                                                                                                                                                                                                                                                                    >  2:479840/1‑289 (MQ=255)
                 aGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTg                                                                                                                                                                                                                                                                    >  2:427388/1‑289 (MQ=255)
                     tGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGCGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGTTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACCCCCTTATTCCAGCTGGCtt                                                                                                                                                                                                                                                                >  2:282667/1‑289 (MQ=255)
                     tGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAac                                                                                                                                                                                                                                                                                    >  1:362940/1‑269 (MQ=255)
                     tGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCt                                                                                                                                                                                                                                                                 >  2:59694/1‑288 (MQ=255)
                     tGCGGCATTTTGTCCGCGCCGGGCTTCGCCCACTGTTCAGGCCGGAGCCACAGACCGTCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAACGCGATGGGCAGCGACTACATCCGTGTGGTGAATGTGGTGAAGTCTGCCAGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACCTTATAGACCATAACCAGCTCCACACCTTTACTCACGCTGGCtt                                                                                                                                                                                                                                                                >  2:235177/1‑289 (MQ=255)
                     tGCGGCACTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAac                                                                                                                                                                                                                                                                                    <  2:362940/269‑1 (MQ=255)
                          cATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTgg                                                                                                                                                                                                                                                           >  1:596685/1‑289 (MQ=255)
                           aatttGTCCGCGCCGGGCTTCCCTCACTGTTCAGGCCCGAGGCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTgg                                                                                                                                                                                                                                                           <  2:135584/286‑1 (MQ=255)
                           aTTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTgg                                                                                                                                                                                                                                                           <  2:650452/288‑1 (MQ=255)
                            ttttGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCt                                                                                                                                                                                                                                                         <  2:25268/289‑1 (MQ=255)
                            ttttGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCt                                                                                                                                                                                                                                                         <  1:85890/289‑1 (MQ=255)
                                                                     cacaGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGa                                                                                                                                                                                                                <  1:59694/289‑1 (MQ=255)
                                                                         gACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGatcat                                                                                                                                                                                                            <  2:70284/289‑1 (MQ=255)
                                                                           ccgccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGatcat                                                                                                                                                                                                            <  1:479840/287‑1 (MQ=255)
                                                                              ccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCgggg                                                                                                                                                                                                       >  2:439879/1‑289 (MQ=255)
                                                                                        ggCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCCTCGGGGTAAAACGGtt                                                                                                                                                                                             >  2:119140/1‑289 (MQ=255)
                                                                                         gCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTaaaa                                                                                                                                                                                                  >  2:444516/1‑283 (MQ=255)
                                                                                         gCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTaaaa                                                                                                                                                                                                  <  1:444516/283‑1 (MQ=255)
                                                                                              tGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCGTCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGAc                                                                                                                                                                                       >  1:452677/1‑289 (MQ=255)
                                                                                                          tctcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAggg                                                                                                                                                                           >  1:219444/1‑289 (MQ=255)
                                                                                                             cccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAc                                                                                                                                                                        <  2:87237/289‑1 (MQ=255)
                                                                                                                gAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCTTAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAcgcg                                                                                                                                                                     <  2:86457/289‑1 (MQ=255)
                                                                                                                gAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAcgcg                                                                                                                                                                     <  2:693548/289‑1 (MQ=255)
                                                                                                                      aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAACCGGTTTATGACGTAGCGTCCGGGGACGCGttattt                                                                                                                                                               >  1:300955/1‑289 (MQ=255)
                                                                                                                      aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTAATGACGTAGCGTCAGGGGACGCGgtattt                                                                                                                                                               >  1:617198/1‑289 (MQ=255)
                                                                                                                                   gAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGCGACTGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTTGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTCAACTCg                                                                                                                                                  >  1:477751/1‑289 (MQ=255)
                                                                                                                                    aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACt                                                                                                                                                    <  1:439879/286‑1 (MQ=255)
                                                                                                                                         cgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGt                                                                                                                                            >  1:681556/1‑289 (MQ=255)
                                                                                                                                          gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGtt                                                                                                                                           >  1:526899/1‑289 (MQ=255)
                                                                                                                                                      tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGt                                                                                                                               >  2:158505/1‑289 (MQ=255)
                                                                                                                                                        cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa                                                                                                                                                                                                                                                                                     >  1:454415/1‑137 (MQ=255)
                                                                                                                                                        cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa                                                                                                                                                                                                                                                                                     <  2:454415/137‑1 (MQ=255)
                                                                                                                                                        cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGtt                                                                                                                              >  2:461360/1‑288 (MQ=255)
                                                                                                                                                         ccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAc                                                                                                                            <  2:404837/289‑1 (MQ=255)
                                                                                                                                                          cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAcc                                                                                                                           <  1:339532/289‑1 (MQ=255)
                                                                                                                                                             tCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGTTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACATTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTCTTCAACAGTGATGAGTCATCGTGGCATCGCGTTGACGATCTTCGGGGTCAAACGGTTTATGACCTAGCGCCATGGGACGCGTTATTTCATTCGTAACTCGTTTCTTTCACGGAAAATGTTAcattg                                                                                                                        >  2:417676/1‑285 (MQ=255)
                                                                                                                                                                         catGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACCGTGATGAGTCATCGTGGCATCTCGTTGACGATCATCGGGGTAATACGGTTAATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGTTCCGTGAGCGGAAAATGTTACATGGTTATCGACGGaa                                                                                                            >  2:558689/1‑289 (MQ=255)
                                                                                                                                                                              aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGa                                                                                                       >  1:717568/1‑289 (MQ=255)
                                                                                                                                                                              aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGa                                                                                                       >  1:696273/1‑289 (MQ=255)
                                                                                                                                                                                  gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGttt                                                                                                   >  2:242317/1‑289 (MQ=255)
                                                                                                                                                                                         atcGACTACAGCCGTGAGGGGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                                                                                             <  2:596685/286‑1 (MQ=255)
                                                                                                                                                                                         aGCGACTACAGCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCGAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                                                                                             <  2:703773/288‑1 (MQ=255)
                                                                                                                                                                                           cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTgg                                                                                          >  1:548255/1‑289 (MQ=255)
                                                                                                                                                                                            gACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGa                                                                                         <  2:230218/289‑1 (MQ=255)
                                                                                                                                                                                             actactTCCGTGAGGGGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTTTAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGaa                                                                                        <  1:490851/289‑1 (MQ=255)
                                                                                                                                                                                              cTACATCCGTGAGGTGAATGTGGTGAAGTCTTGCCGTTTCTGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGaa                                                                                        <  2:219444/288‑1 (MQ=255)
                                                                                                                                                                                              cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGaa                                                                                        <  1:455938/288‑1 (MQ=255)
                                                                                                                                                                                                  aTCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa                                                                                   <  1:461360/289‑1 (MQ=255)
                                                                                                                                                                                                       tGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa                                                                                   >  2:488177/1‑284 (MQ=255)
                                                                                                                                                                                                       tGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa                                                                                   <  1:488177/284‑1 (MQ=255)
                                                                                                                                                                                                       tGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGc                                                                               >  2:704519/1‑288 (MQ=255)
                                                                                                                                                                                                        gAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGcc                                                                              <  1:158505/288‑1 (MQ=255)
                                                                                                                                                                                                        agggTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCTGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCt                                                                             <  1:235177/287‑1 (MQ=255)
                                                                                                                                                                                                         aGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTg                                                                            >  2:436174/1‑289 (MQ=255)
                                                                                                                                                                                                           gTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTggg                                                                          >  2:44634/1‑289 (MQ=255)
                                                                                                                                                                                                           gTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTgg                                                                           <  2:109986/288‑1 (MQ=255)
                                                                                                                                                                                                                            tGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATgcagaagca                                                         >  1:123820/1‑289 (MQ=255)
                                                                                                                                                                                                                            tGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCGGGGTGAAAGATgcagaagca                                                         >  2:560523/1‑289 (MQ=255)
                                                                                                                                                                                                                              cccGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGACAGATgcagacgcaga                                                       >  2:205804/1‑289 (MQ=255)
                                                                                                                                                                                                                                   gtCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCTTAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGaaaa                                                    >  1:19934/1‑287 (MQ=255)
                                                                                                                                                                                                                                   gtCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAAc                                                  >  2:405074/1‑289 (MQ=255)
                                                                                                                                                                                                                                                             tgtTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGCTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACTGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGAGgaagac                        >  1:686233/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                      cTACTTTATAGAGCTTAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACaaaaaa                <  2:19934/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACaaaaaa                <  2:681556/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACa              >  1:445094/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         cTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGACGACACAAAAAACAg             >  2:114153/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                         cTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGc            <  1:242317/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ataGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGa        <  1:427388/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ataGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGa        <  2:123820/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gagCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGATGc     >  1:144557/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 agCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGCCGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGCTCCGGATTCGGGAGGCGGCAGAAACAAAAACAAGCCTGAGGCa    >  2:322269/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   cATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAACACAGCCTGATGCAgg  >  1:574947/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   cATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGATGCAgg  >  2:399527/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   cATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGATGCAgg  >  1:526945/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   cATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAAAACAAAACAGACTGATGCAgg  >  1:382938/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                               
CGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCCTGGGTGAAAGATGCAGAAGCAGAAAAACTGTTCCGGATTCGGGAGGCGGAAGAAACAAAAAACAGCCTGATGCAGG  >  NC_000913/581364‑581927

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