Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 588,640 T→C 5.2% Y772C (TAT→TGT)  nfrA ← bacteriophage N4 receptor, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913588,6400TC 5.2% 263.6 / 2.8 78Y772C (TAT→TGT) nfrAbacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (41/32);  new base C (2/2);  total (43/35)
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.13e-01

TCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGA  >  NC_000913/588352‑588927
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                
tcgCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGAt                                                                                                                                                                                                                                                                                                  >  1:602407/1‑289 (MQ=255)
     gTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCGGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCtt                                                                                                                                                                                                                                                                                             <  1:133554/289‑1 (MQ=255)
          tGCCGCCATTAAGGGATGAGGCGCTGGCGGCCAGCATGGTTTCGTGTGGGCCGTTGTGGCCGTTCGGCGGCAACTGCTGTTCGACAGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCGCCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTccc                                                                                                                                                                                                                                                                                        <  1:53104/289‑1 (MQ=255)
              gccATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGgc                                                                                                                                                                                                                                                                                    >  2:549067/1‑289 (MQ=255)
                          tGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCATTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTTTCCGTAGCTACGATAGCTTTTCCCTTGCGCTGCGccgccg                                                                                                                                                                                                                                                                        >  1:294590/1‑289 (MQ=255)
                              gcgcTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACa                                                                                                                                                                                                                                                                     <  2:204018/288‑1 (MQ=255)
                                              tGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGAAAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACt                                                                                                                                                                                                                                                    <  2:186115/289‑1 (MQ=255)
                                                      gATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTccc                                                                                                                                                                                                                                                                                        >  2:565547/1‑245 (MQ=255)
                                                      gATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTccc                                                                                                                                                                                                                                                                                        <  1:565547/245‑1 (MQ=255)
                                                      gATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGACAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTgcg                                                                                                                                                                                                                                             >  2:62106/1‑288 (MQ=255)
                                                      gATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCA‑CCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGcc                                                                                                                                                                                                                                           >  2:402003/1‑289 (MQ=255)
                                                           gcCATTTTTGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTCCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCATCAGCATATCGCGTCCGATGCGGTACGCGGCTTTCCGTTGCCCGTAGCTACGATCGTTTCTCTCTGGCGCTTCGCCGCCGCCATTCTTGTTCGCTGGACTCCTTACGCCGGaa                                                                                                                                                                                                                                       >  1:407810/1‑289 (MQ=255)
                                                           gcCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCTCCCGACATTATTCTTAGCGGTACTACTTGCCCCggc                                                                                                                                                                                                                                        >  1:372456/1‑287 (MQ=255)
                                                           gcCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGa                                                                                                                                                                                                                                        >  1:66913/1‑288 (MQ=255)
                                                                tttGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCaa                                                                                                                                                                                                                                  <  2:602407/289‑1 (MQ=255)
                                                                         tCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGa                                                                                                                                                                                                                                                                       >  1:558330/1‑243 (MQ=255)
                                                                         tCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGa                                                                                                                                                                                                                                                                       <  2:558330/243‑1 (MQ=255)
                                                                              ggCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAAt                                                                                                                                                                                                                    <  1:247781/289‑1 (MQ=255)
                                                                                  aCTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGaa                                                                                                                                                                                                                >  1:288783/1‑289 (MQ=255)
                                                                                                   aTGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCCCCGGAACGCAAGCTGATGGAAGAATCGAACCTGCACGTCCAgcggc                                                                                                                                                                                               >  1:95685/1‑289 (MQ=255)
                                                                                                        aaaGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcgaccg                                                                                                                                                                                          <  2:672064/289‑1 (MQ=255)
                                                                                                         aaGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTCTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcgaccga                                                                                                                                                                                         <  1:330848/289‑1 (MQ=255)
                                                                                                                gATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGAcctccac                                                                                                                                                                                  >  2:104555/1‑287 (MQ=255)
                                                                                                                           ggCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAAcggc                                                                                                                                                                       <  1:549067/289‑1 (MQ=255)
                                                                                                                                     gcAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCA‑CCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGa                                                                                                                                                              <  1:402003/287‑1 (MQ=255)
                                                                                                                                                   gACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGAtttt                                                                                                                                               >  2:406095/1‑289 (MQ=255)
                                                                                                                                                     cATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGt                                                                                                                                             >  2:524526/1‑289 (MQ=255)
                                                                                                                                                        cGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACatt                                                                                                                                                                                                                                                                   <  1:46491/168‑1 (MQ=255)
                                                                                                                                                        cGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACatt                                                                                                                                                                                                                                                                   >  2:46491/1‑168 (MQ=255)
                                                                                                                                                        cGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGaa                                                                                                                                                                                                                <  1:592992/219‑1 (MQ=255)
                                                                                                                                                        cGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGaa                                                                                                                                                                                                                >  2:592992/1‑219 (MQ=255)
                                                                                                                                                        cGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGt                                                                                                                                             >  2:190610/1‑286 (MQ=255)
                                                                                                                                                                  aCCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTg                                                                                                                                <  2:407810/289‑1 (MQ=255)
                                                                                                                                                                            tcaCCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCa                                                                                                                       >  1:29888/1‑288 (MQ=255)
                                                                                                                                                                              aCCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAgcg                                                                                                                    <  2:95685/289‑1 (MQ=255)
                                                                                                                                                                                cccGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAgcg                                                                                                                    >  2:228611/1‑287 (MQ=255)
                                                                                                                                                                                 ccGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGTTAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCt                                                                                                                 >  2:224997/1‑289 (MQ=255)
                                                                                                                                                                                 ccGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCt                                                                                                                 <  2:66913/289‑1 (MQ=255)
                                                                                                                                                                                    ttttCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCCGCGGCCACCGACCTCCTCATGCAAACGGCGGCAATTGAAGCGTCGTGGATTTTGTTCTGGGGTCCGTGGGGTTATCAGCGCATGAt                                                                                                              >  1:461139/1‑289 (MQ=255)
                                                                                                                                                                                    ttttCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcgaccga                                                                                                                                                                                         >  1:280344/1‑214 (MQ=255)
                                                                                                                                                                                    ttttCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcgaccga                                                                                                                                                                                         <  2:280344/214‑1 (MQ=255)
                                                                                                                                                                                             gTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGt                                                                                                     >  1:507380/1‑289 (MQ=255)
                                                                                                                                                                                              tATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGt                                                                                                     >  1:287135/1‑288 (MQ=255)
                                                                                                                                                                                              tATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGAATTTGTTCTGGGGTCCGTGGGGTTATCCGCTCCTGATTATCAATGtc                                                                                                    >  1:287131/1‑289 (MQ=255)
                                                                                                                                                                                                   gCAAAGACGAGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCAGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATGATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAACCGAAACTGCACGTACAGCGGCGACCGCCCTACTCATGCCAACGGGGCAAGTCGAAGAGTCGTCGGTTTGGTTCCGGGGTCAATGGGGTTATCACCGTCTGATTATCAATGccaataa                                                                                               >  1:331577/1‑285 (MQ=255)
                                                                                                                                                                                                       aGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCAc                                                                                           <  2:507380/289‑1 (MQ=255)
                                                                                                                                                                                                       aGACGCGGCTATAAACTGAGAGCACGTCGCCTTCCAGCAGCATATTGCGTCCGACGCTGTACTCGGCTTCCAGTTGTCCGTAGCTCCGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACCCATTCCGCCGGCACGCAAGCCGATGCAATAATCGAAACTGACCGCCCACCGGCGACCGACCCCCTCATCCACACGGCGAAGATTGAACCGTTGTGGGTATTGCTCTGCGGCCAGCGGGATTATCAGCGCCTgtgtcgtcatgtccttaatct                                                                                            >  1:116855/1‑268 (MQ=255)
                                                                                                                                                                                                        gACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCAc                                                                                           >  1:218455/1‑288 (MQ=255)
                                                                                                                                                                                                            cgGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGACAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTAt                                                                                                                         >  1:686595/1‑254 (MQ=255)
                                                                                                                                                                                                            cgGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGACAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTAt                                                                                                                         <  2:686595/254‑1 (MQ=255)
                                                                                                                                                                                                             ggCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGcc                                                                                      <  2:218455/288‑1 (MQ=255)
                                                                                                                                                                                                             ggCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGcc                                                                                      <  2:29888/288‑1 (MQ=255)
                                                                                                                                                                                                                 ataaacttgaggCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGAGGTACTCGGCTTCCAGTTGTCCGTAGCGACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCg                                                                                 <  1:685772/278‑1 (MQ=255)
                                                                                                                                                                                                                 ataAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCg                                                                                 <  1:543490/289‑1 (MQ=255)
                                                                                                                                                                                                                 ataAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCg                                                                                 <  2:288783/289‑1 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGa                                                                                                                                                                                                           >  2:47957/1‑159 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGa                                                                                                                                                                                                           <  1:47957/159‑1 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTCTCAATGTCCTCAATCACCAGCAGGGCGTAGtgct                                                                          >  2:308368/1‑288 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCa                                                                                                         >  1:103943/1‑257 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCa                                                                                                         <  2:103943/257‑1 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGtgct                                                                          >  2:687536/1‑288 (MQ=255)
                                                                                                                                                                                                                         agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGtgct                                                                          >  2:198291/1‑288 (MQ=255)
                                                                                                                                                                                                                          gagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGGGCTGGTGATTTTGTTCTGGGGTCAGTGGGGTTCTCAGCGCATGATTATCAATGTCATCAATCACACTCCGGGCGTAGtgctg                                                                         >  1:427821/1‑288 (MQ=255)
                                                                                                                                                                                                                                     cTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCAt                                                             <  1:228611/289‑1 (MQ=255)
                                                                                                                                                                                                                                                 tatTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACg                                                <  2:337928/289‑1 (MQ=255)
                                                                                                                                                                                                                                                      cGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCt                                                                                                                 >  1:247183/1‑220 (MQ=255)
                                                                                                                                                                                                                                                      cGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCt                                                                                                                 <  2:247183/220‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTc                                         >  2:78245/1‑288 (MQ=255)
                                                                                                                                                                                                                                                            aTGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCCTGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCGGGGGTCAGTGGGGTTATCAGCGCCTGATTATCCATGTCATCAATCACCAGCCTGGCGTTGTGCTTGGTCGCGGGCATGTC‑ATCCAGACTCTGGTTCACGt                                     >  1:677086/1‑289 (MQ=255)
                                                                                                                                                                                                                                                             tGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGACAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTa                                    <  1:62106/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ggCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCCCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCCGTGGGGTTAACAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATTTCTACCCAG‑TGCCGGTTCACGTCGGCCAGtttt                          >  1:461188/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                         cTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCCATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGCC‑ATCCAGACGCTGGATCACGTAGGCCAGTTGTCg                        >  2:151726/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         cTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCg                        >  2:479877/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                            ccAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGa                      >  1:637528/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                            ccAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGAt                     <  2:287131/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGAt                     <  2:287135/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      gTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGa           <  2:329153/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATc         <  2:425508/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               aTAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGa  >  2:290113/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                                tAGCTTGTCCCTGGCCCGGCGTGGCCGAAAGTAGTGTTGGCGGGACTGATTGCGCCGGAACGCAGACCGTAGAAAGAATCGAAGCAGAACGGCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGGTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCGTGATTATCAAGGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGa  <  1:78519/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                
TCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATTTTGGCCGTTCAGCGGCAACTGCTGTTCGACGGCGATGAAAAAGATCTGATCGCGCAGCGGCTTCCAGCGCAGACCGGTGCCGGACATCGGATTTTTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTC‑ATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGA  >  NC_000913/588352‑588927

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