Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,926,900 T→G 7.7% D647A (GAC→GCC)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,9000TG 7.7% 368.7 / 22.5 128D647A (GAC→GCC) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base T (49/67);  new base G (5/6);  total (54/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-03

Only 100 of 133 total aligned reads displayed.

GTTGATGACGACCTTCACCTTTGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTAC  >  NC_000913/1926626‑1927156
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                
gtTGATGACGACCTTCACCTTTGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                                                                                                                                                                                                                                  <  2:221218/275‑1 (MQ=255)
     tgtcgacCTTGACCTTTGGAAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                                                                                                                                             <  4:208122/272‑1 (MQ=255)
      gacgacCTTTACCATTGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCAGGGCTATCTGCGCTGGTAAGTTAATAACCTGGACAGCATGAGGCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                                                                                                                                                                                             <  2:444070/274‑1 (MQ=255)
        ctcctttCACCATTGGCCAATACCGGGGCAACGCTGTTTCTGCCGTTGCCGAACGCGATCCGGGCTTTCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATGATCTTGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag                                                                                                                                                                                                                                                          <  4:537812/270‑1 (MQ=255)
                    ttGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                                                                                                                                                                                                                              <  1:39401/275‑1 (MQ=255)
                      ggCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                                                                                                                                            >  2:531254/1‑275 (MQ=255)
                                 tGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCCTATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATTGCCTGATTCCATGTCGGTACAGAGCCATAAAAGATGGTCATCGGttt                                                                                                                                                                                                                                 >  4:70423/1‑275 (MQ=255)
                                       aGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATGTAATAACCTGGACAGCATGAGCCCGGAGCTGCGGTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                                                                                                                                                                                            <  2:384544/274‑1 (MQ=255)
                                                      tGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCACTAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGa                                                                                                                                                                                                            >  2:272583/1‑275 (MQ=255)
                                                             cgcgATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCNAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCANGGCCTGAGTCCATTTCGGTACCGAGCAATAAAAAATGGCAATCGGNTNTCCGCTCGCCCAATTTAGCCAACCCttt                                                                                                                                                                                                     >  3:43083/1‑275 (MQ=255)
                                                                 aTCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                                                                                                                                                                                  <  2:170215/274‑1 (MQ=255)
                                                                               cgcTGTTTATTTAAAAACCTTGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGACCGCAGACGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATTTTCCATCGCTACGCCTTCGTACGAGTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTGAGCGACCCATTTTGCCGGTTCCCAATATTg                                                                                                                                                                                   <  2:520397/275‑1 (MQ=255)
                                                                                gcTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACCCTGAAACATTATCTGGAAAATACTTAGGCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                                                                                                                                                                                   <  2:29532/274‑1 (MQ=255)
                                                                                gcTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                                                                                                                                                                                   <  2:39063/274‑1 (MQ=255)
                                                                                gcTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                                                                                                                                                                                   <  2:356986/274‑1 (MQ=255)
                                                                                gcTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCAGATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATATAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                                                                                                                                                                                   <  4:310077/274‑1 (MQ=255)
                                                                                     ttATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTAGCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCt                                                                                                                                                                              <  1:115409/274‑1 (MQ=255)
                                                                                              aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATctt                                                                                                                                                                     >  2:399159/1‑272 (MQ=255)
                                                                                                   ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAACAGCTGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTTCACCTGGGAATCTTTCAAAc                                                                                                                                                               >  2:16202/1‑275 (MQ=255)
                                                                                                         cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCTAGCGCGACCAGAAAAGCATATTCCATCGCTACGCTTTCGTACGATTTAAAGCGACCAGATTTATCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTATCGACCCATTTTGCCGGTTCCCAATATTCCACCTGAGAACCGTGCAAACCGGTCg                                                                                                                                                         >  2:188150/1‑275 (MQ=255)
                                                                                                         cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                                                                                                                                                                                          >  2:356087/1‑194 (MQ=255)
                                                                                                             aGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATc                                                                                                                                                                                                                                      <  4:356086/194‑1 (MQ=255)
                                                                                                              gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCCCCTGAGAATCGTGCAAACCGGTCGTTAcc                                                                                                                                                    >  2:450235/1‑275 (MQ=255)
                                                                                                                 cGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCt                                                                                                                                                                                                                           >  2:101241/1‑201 (MQ=255)
                                                                                                                   gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAg                                                                                                                                                  >  2:369724/1‑272 (MQ=255)
                                                                                                                     gCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                                                                                                                                                                                       <  4:101240/201‑1 (MQ=255)
                                                                                                                        gCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCTACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTTCAAACCGGTCGTTACCAGTCAATGCg                                                                                                                                          >  1:99391/1‑275 (MQ=255)
                                                                                                                        gCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAAt                                                                                                                                             <  4:402842/272‑1 (MQ=255)
                                                                                                                                   tGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCCTTCAAACCGGTCGTTACCAGTAACTGCGGATAAGCCt                                                                                                                                 >  4:505097/1‑273 (MQ=255)
                                                                                                                                     aaaCATTATCTGGAAAAGACTTAGTCCGCAGGCGTAGAGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCTACCAGAGTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAACGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtgtg                                                                                                                              <  2:81616/274‑1 (MQ=255)
                                                                                                                                         aTTAGCTGTAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGt                                                                                                                          <  2:70424/274‑1 (MQ=255)
                                                                                                                                            accTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACg                                                                                                                      <  2:409715/273‑1 (MQ=255)
                                                                                                                                                ggAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaa                                                                                                                                                                                                                                                 <  5:637696/148‑1 (MQ=255)
                                                                                                                                                    aaaaCTTAGTCCGCAGGCGTAGCGGGTAAGGTTCCCTGGGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATa                                                                                                              <  4:450234/272‑1 (MQ=255)
                                                                                                                                                       aCTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATCTTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGAAAATGCGGATAAGCCTGTTCGATGACGTTGCCATACgg                                                                                                           >  4:508708/1‑275 (MQ=255)
                                                                                                                                                        cTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACggg                                                                                                          <  4:277643/275‑1 (MQ=255)
                                                                                                                                                           aGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTTCGGTGACGTTGTACTCCGGGCTg                                                                                                       >  4:25331/1‑275 (MQ=255)
                                                                                                                                                            gtgcgcgGGCGAAGCGGGTAAGTGGCCCCGCGCCAGCGCAACAAAAAAACAAAATTCAATCGCTAGGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCGCGCGCAATTTAGCGACCCAGTTGGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGt                                                                                                      <  4:16202/272‑1 (MQ=255)
                                                                                                                                                              ccGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                                                                                                                                                                                                 <  7:854727/118‑1 (MQ=255)
                                                                                                                                                              ccGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                                                                                                                                                                                                 >  8:854727/1‑118 (MQ=255)
                                                                                                                                                                 cAGGGGTACGGGGTAATGTTCCCGGCGCCCGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGAGTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCtt                                                                                                 <  4:282439/275‑1 (MQ=255)
                                                                                                                                                                 cAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTTCGGTGACGTTGTCATACGGGTTGTAGCt                                                                                                  >  1:126264/1‑274 (MQ=255)
                                                                                                                                                                  aGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCACACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGTTTTAGCttt                                                                                                >  1:127093/1‑275 (MQ=255)
                                                                                                                                                                  aGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAc                                                                                                                                                     >  2:74367/1‑222 (MQ=255)
                                                                                                                                                                      gTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGt                                                                                                                                                 <  4:74366/222‑1 (MQ=255)
                                                                                                                                                                        aGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGt                                                                                          <  1:17911/275‑1 (MQ=255)
                                                                                                                                                                        aGCGAGTAAGTTGCCCTGCGCCAGCGCGACCAGAAAAACATCTTCCATCGGTACGCCTTCGTACGATTTAAAGCGACNAGATTTACCGCCATGGCCTTAGTCCATGTCGGTACAGAGCAATAAAAGATGGGCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGt                                                                                          <  3:126264/275‑1 (MQ=255)
                                                                                                                                                                             tgtaTGTTCCCTGGGCAAGCGAGACAAGAAAAGCTGTGTCCGTCGCGACGCCATCGGACGATCTCAACCGCCCAAATTTACCGCCATGGCCTGCGTCCATGGCGGTAAAGAGCAATAAAAGATGGTCATCTGTTTTGAGCGAGCGCAATGTAGCGACCCATTTTGCCGGTTCCCAATATTGCACGTGAGAAGCGTGCAAGCCGGTCGGTACCGGTAGATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGATTTTCATGTATTCg                                                                                     <  4:44885/272‑1 (MQ=255)
                                                                                                                                                                                    tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaa                                                                                                                                                                                                                                                 >  5:785311/1‑112 (MQ=255)
                                                                                                                                                                                    tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaa                                                                                                                                                                                                                                                 <  6:785311/112‑1 (MQ=255)
                                                                                                                                                                                               gcgcGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                                                                                                                                                                                                 >  2:206928/1‑85 (MQ=255)
                                                                                                                                                                                                       aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATc                                                                                                                                                                                                                                      <  5:215334/104‑1 (MQ=255)
                                                                                                                                                                                                       aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATc                                                                                                                                                                                                                                      >  6:215334/1‑104 (MQ=255)
                                                                                                                                                                                                         aaaaGCATTCTCCAGCGTGACGCCTTGGTNCAATTTAAAGCCACCAGATTTACCGNCCTGGCCTGAGTCCAGGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTNCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTAcccc                                                         <  3:99391/275‑1 (MQ=255)
                                                                                                                                                                                                                     cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATGTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACt                                             <  4:369723/275‑1 (MQ=255)
                                                                                                                                                                                                                        cGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTTTAGCTTTTCATGTCCTCGTATTCTTCCGGATCCTCCGGGTTACCCCACTCTTCAAACCCa                                           >  4:157598/1‑274 (MQ=255)
                                                                                                                                                                                                                        cGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCAt                                                                                            <  2:298587/225‑1 (MQ=255)
                                                                                                                                                                                                                         gCTACGCCTTCGTACGATTTAAAGCGACCAGATTTCCCGCCATGGCCGGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAcc                                         <  1:43083/275‑1 (MQ=255)
                                                                                                                                                                                                                         gCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCAcc                                         <  2:508709/275‑1 (MQ=255)
                                                                                                                                                                                                                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                                                                                                                                                                                  <  4:170244/120‑1 (MQ=255)
                                                                                                                                                                                                                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg                                                                                                                                                                      >  5:432732/1‑148 (MQ=255)
                                                                                                                                                                                                                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTAACGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg                                                                                                                                                                      <  6:432732/148‑1 (MQ=255)
                                                                                                                                                                                                                            aggCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACt                                             <  2:157599/266‑1 (MQ=255)
                                                                                                                                                                                                                            aCGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAttt                                                                                                                                                                                                                 <  6:493068/104‑1 (MQ=255)
                                                                                                                                                                                                                            aCGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAttt                                                                                                                                                                                                                 >  5:493068/1‑104 (MQ=255)
                                                                                                                                                                                                                             cGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTg                                     >  2:123344/1‑275 (MQ=255)
                                                                                                                                                                                                                              gCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                                                                                                                                                                                             <  7:29500/106‑1 (MQ=255)
                                                                                                                                                                                                                               ccTTCCTACGATTAAAAGCGAACAGAGTTACGCACATGGCCTGAGTCCATGTCGGTAAAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCAACTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGt                                   <  4:297891/275‑1 (MQ=255)
                                                                                                                                                                                                                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                                                                                                                                                                                               <  6:742657/68‑1 (MQ=255)
                                                                                                                                                                                                                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                                                                                                                                                                                  <  5:135646/129‑1 (MQ=255)
                                                                                                                                                                                                                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                                                                                                                                                                                  >  6:135646/1‑129 (MQ=255)
                                                                                                                                                                                                                                    gTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATTTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAagag                              >  2:258628/1‑275 (MQ=255)
                                                                                                                                                                                                                                      aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                                                                                                                                                                               <  7:1052496/80‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg                                                                                                                                                                      >  8:320543/1‑137 (MQ=255)
                                                                                                                                                                                                                                      aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg                                                                                                                                                                      <  7:320543/137‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGACATCGGTTTTc                                                                                                                                                                                                                               >  8:1052496/1‑80 (MQ=255)
                                                                                                                                                                                                                                            tacaGCAACCATATTTACCGCCATGGCCTGAGTCCATGTCGGAACAGAGCAATAAAAGATGGCCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCAAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGAt                      <  2:487193/272‑1 (MQ=255)
                                                                                                                                                                                                                                            tAAAGCGATCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGAt                      <  2:410109/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                                                                                                                                                                                             <  8:563628/85‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                                                                                                                                                                                             >  7:563628/1‑85 (MQ=255)
                                                                                                                                                                                                                                                     ccgATTTCCCGCCATTGCCTGAGTCCATTTCGGTACAGCGCAATAAAAGATTGTCATCCGATTTCAGCTCGCGCAATTTTGCGACCCCTTTTGCCTGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGc             <  4:87359/273‑1 (MQ=255)
                                                                                                                                                                                                                                                     cAGAATTAACGCCATGGCCTGAGTCCATGTCGGTACAGGGCAATAAAAGATGTGCATCGGTTGTCAGCGCGTGCAATTGAGAGACCCATGTTGCCGGTTCCCAAGATTGCAACTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGc             <  4:258627/275‑1 (MQ=255)
                                                                                                                                                                                                                                                       gATTTACCGCCATGGCCGTAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGAAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAat                                                                                 <  1:54385/205‑1 (MQ=255)
                                                                                                                                                                                                                                                         tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCAcc                                                                                                                                                                               >  5:353031/1‑109 (MQ=255)
                                                                                                                                                                                                                                                         tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCg         <  4:123343/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                                                                                                                                                 >  6:356189/1‑87 (MQ=255)
                                                                                                                                                                                                                                                              cgacATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTa                                                                                                                                                      <  7:838848/126‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                                                                                                                                                                                  <  6:983731/101‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCt                                                                                                                                                                              <  8:578794/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccctggccgtgttCCATTTTGGNNCATGGNAATAACAGTTGGGCACCCGTTTTCGGGCGCCGACATTTTGCGACCCCTTTTGCCGGTTCCTAATCTTGGACCTGTGAAGCGTGCAATCCGGTCGTTACCAGGAGATGCGGATAAGCCTGTGCGTTGTCGTTGTCATACGTGCTGTAGCTTTTCATGTACGCGAAATATTGCGGATCCTGCGGGTTTCCCCACTCTTCAAACTCCCCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTac  <  3:5413/263‑1 (MQ=255)
                                                                                                                                                                                                                                                                acATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                                                                                                                                    >  8:1144004/2‑129 (MQ=255)
                                                                                                                                                                                                                                                                acATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATg                                                                                                                                            >  6:1126699/2‑137 (MQ=255)
                                                                                                                                                                                                                                                                      ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                                                                                                                                    >  7:872462/1‑123 (MQ=255)
                                                                                                                                                                                                                                                                      ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                                                                                                                                    <  8:872462/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tcaggccATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGt                                                                                                                                                        >  5:793105/6‑117 (MQ=255)
                                                                                                                                                                                                                                                                        tcaggccATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGt                                                                                                                                                        <  6:793105/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCaa                                                                                                                                                                 <  6:904482/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCaa                                                                                                                                                                 >  5:904482/1‑103 (MQ=255)
                                                                                                                                                                                                                                                                                  tCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAat                                                                                                                                                                                       >  8:901112/1‑76 (MQ=255)
                                                                                                                                                                                                                                                                                  tCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAg                                                                                                                                    <  7:911080/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                
GTTGATGACGACCTTCACCTTTGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTAC  >  NC_000913/1926626‑1927156

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