Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,507 A→C N720H (AAC→CAC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5070AC97.6% 265.5 / NA 83N720H (AAC→CAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (38/42);  total (40/43)
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.71e-01

TGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGA  >  NC_000913/4187209‑4187802
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       
ggCGGACCAGAGCATGGACAACGGCTTCGCCTAGTCAGCGCGTTCTGGTGCAGCTGTTGGTATCGAGGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGGGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCAc                                                                                                                                                                                                                                                                                                       <  1:119494/300‑1 (MQ=255)
        aGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGGGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCt                                                                                                                                                                                                                                                                                                >  1:674057/1‑300 (MQ=255)
          cctgcacgtcgcgcttccgtactgtcAGCGCGTTCTGGTGCCTCTTTTGGTATCGAAGAAATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCCAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCCGGCGAACGCTACAACTAAGTTATCGATCTCTGGGCTGCGGCGAACGCTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAGCCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACa                                                                                                                                                                                                                                                                                             <  1:671418/276‑1 (MQ=255)
              tGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTAGCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATGTACatgat                                                                                                                                                                                                                                                                                         >  1:550839/1‑301 (MQ=255)
                   cccGGCTTCACCTATGCCGCGCGATTGGGGGCATCTGTGGGTAGGGATGACATGGTCATCCCGGAGAAGAGACACGAAATCATCTCCGAGGGAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCGGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCg                                                                                                                                                                                                                                                                                    <  1:167056/300‑1 (MQ=255)
                     cGGTTTCGCCTATGCAGCGCGTTCTGGGGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGAc                                                                                                                                                                                                                                                                                  <  1:323406/301‑1 (MQ=255)
                                 tGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATATCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCTGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATAAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTCCa                                                                                                                                                                                                                                                                                             >  2:671603/1‑278 (MQ=255)
                                   cAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTgg                                                                                                                                                                                                                                                                    <  2:691100/301‑1 (MQ=255)
                                    aGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGCTAGGAAGAGAAGCAGGTTACCTTCGACAGTATCAACGTGATGGACGACTCCGGTGCGCGTAGt                                                                                                                                                                                                                                                                   >  1:257617/1‑301 (MQ=255)
                                              ggTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCCGcac                                                                                                                                                                                                                                                         >  2:301526/1‑301 (MQ=255)
                                                ggCAACTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAAAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAg                                                                                                                                                                                                                                                       <  1:301436/300‑1 (MQ=255)
                                                      tGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGt                                                                                                                                                                                                                                                 >  2:466273/1‑301 (MQ=255)
                                                      tGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACCGAAACCGTGATTAACCGTTACGGTCAGGAAGAGAAGCAGGTTTCCTTCCAACGCCTGTACATGATGGGCGACCGTGTAGAGAGTGGTTCTGGGGCAAATATTCgg                                                                                                                                                                                                                                                 >  1:151602/1‑300 (MQ=255)
                                                      tGTTGGTAGCGATGATATGGTGAACCAGGAGAAGAAACACGAAAACATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCg                                                                                                                                                                                                                                                  <  2:151649/300‑1 (MQ=255)
                                                               cGATGACATGGTCAACCCGGAGAAGAAACACGAAAGCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                                                                                                                                                                                                                                         <  1:636532/300‑1 (MQ=255)
                                                                gATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGACGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTATGAGGCACAGCTTCTTCAGCTTGCTg                                                                                                                                                                                                                                       >  2:465460/1‑301 (MQ=255)
                                                                 aTGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg                                                                                                                                                                                                                                      <  2:473993/301‑1 (MQ=255)
                                                                           cATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACCGATTCGTCAGCTTGCTGGGAGTCGTGGt                                                                                                                                                                                                                            >  2:673086/1‑301 (MQ=255)
                                                                            aTCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                                                                                                                                                                                                                                            >  1:247522/1‑284 (MQ=255)
                                                                            aTCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                                                                                                                                                                                                                                            <  2:247601/284‑1 (MQ=255)
                                                                            aTCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCGACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGCTTCGTCAGCTTGCTGGTATGCGTGGtt                                                                                                                                                                                                                           >  2:72855/1‑300 (MQ=255)
                                                                             tCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCt                                                                                                                                                                                                                          >  2:53062/1‑301 (MQ=255)
                                                                               ccGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAACCAGGTTTCCTTCAACAGCATCTACATGCTGGCCGACTCCGGAGCGCGTTGTTCTGGGGCACAGATTCGTCAGCTTGCCGGTATGCGTTGTCttg                                                                                                                                                                                                                        >  1:703864/1‑299 (MQ=255)
                                                                                cGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTCCGCGTGGTTCTGCGGCACAGGTCCGTCAGCTTGGTGGATTGCGTGGTCTGAt                                                                                                                                                                                                                       >  1:92201/1‑301 (MQ=255)
                                                                                           cacGAAATCCACTCCGAGGCAGACGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCgg                                                                                                                                                                                                            <  2:550995/301‑1 (MQ=255)
                                                                                             cGAAATCCTCCCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGAt                                                                                                                                                                                                          <  2:92224/301‑1 (MQ=255)
                                                                                                       tcCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATCCGTAAGCTTGCTGGTATGCGTGGTCGGATGGCGAAGCCGGATGGCTCcatc                                                                                                                                                                                                 >  1:577196/1‑300 (MQ=255)
                                                                                                        ccGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCCGCTTGCTGGTATGCGTGGTCTGGTGGCGAAGCCGGATGGCTc                                                                                                                                                                                                     >  1:644043/1‑295 (MQ=255)
                                                                                                        ccGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAGGCCGGAGGGCTCcatc                                                                                                                                                                                                 >  2:258416/1‑299 (MQ=255)
                                                                                                        ccGAGGAAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCGGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTc                                                                                                                                                                                                     <  2:644219/295‑1 (MQ=255)
                                                                                                          ggggCAGAAGCAGAAGTTGCTGAGATTCAGGGGCAGTTCCAGTCTGGTTTGGTAACTGCGGGCGCAAGCTACAGCAAAGTTATCGTTATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGGGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCg                                                                                                                                                                                             <  1:53052/299‑1 (MQ=255)
                                                                                                          gAGGCAGAATCAGAAAGTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACGGCGGTCGAACGCTACAAAAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCg                                                                                                                                                                                             <  2:257698/301‑1 (MQ=255)
                                                                                                                       aagTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGCTGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCAc                                                                                                                                                                                >  2:360753/1‑301 (MQ=255)
                                                                                                                        ccttGCTGACATTCAGGGGCGGTTCACTTCGGTTCAGGTAACGGGGGACGAACGCTACAACAAAGTTAGCGTGATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGGGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCAcc                                                                                                                                                                               <  1:353155/299‑1 (MQ=255)
                                                                                                                           tGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTAAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGCCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGAGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGTAACGACAATctcgagg                                                                                                                                                                            >  2:353260/1‑295 (MQ=255)
                                                                                                                                      aGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGAGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAc                                                                                                                                                                                         <  2:344087/277‑1 (MQ=255)
                                                                                                                                      aGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGAGGTCTGATGGCGAAGCCGGATGGCCCCATCATCGAAAc                                                                                                                                                                                         >  1:343982/1‑277 (MQ=255)
                                                                                                                                           cAGTTCCCGTCTGGTCTGGTCACTGCGGGCGAACGCTAAAACAAAGTTATCGATATCGGGGCTGCGGCGAACGATCGTGTATCAAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTc                                                                                                                                                            <  1:466136/301‑1 (MQ=255)
                                                                                                                                            aGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCgg                                                                                                                                                                                                                                                                             <  2:277989/187‑1 (MQ=255)
                                                                                                                                            aGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCgg                                                                                                                                                                                                                                                                             >  1:277906/1‑187 (MQ=255)
                                                                                                                                               tCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGGTGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGTTCTgc                                                                                                                                                         >  2:568452/1‑299 (MQ=255)
                                                                                                                                                ccAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTACCCGTGACGGTCAGGAAGAGAAGCCGGTTTCCTTCCACAGCATCTACATGATGGCAGACTCCCGTGCGCGTGGTTCTGCGGCACACATTCGTCAGCTTGCTGGTATGCGTGGTATCTTGGCGTAGCAGCGTGGCCCCCTCATCCAAACGCCAATCAcctcatccttccgcgcatgcatgac                                                                                                                                                        >  1:89806/1‑277 (MQ=255)
                                                                                                                                                          gttggaaacgCGGGCGAACGCTACAACAAGGTTATCGATATCTGGGTGGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGt                                                                                                                                             <  1:672900/292‑1 (MQ=255)
                                                                                                                                                           ggggaaCTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCTAACGATCGTGTATCCAAGGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGGTGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                                                                                                                            <  1:568290/297‑1 (MQ=255)
                                                                                                                                                               acatgagtgCGAAAGCTACAACAAAGTGATCGAGATCTGGGCTGCAGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACtt                                                                                                                                         <  1:584749/292‑1 (MQ=255)
                                                                                                                                                                     gggCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                                                                                                  >  2:112919/1‑301 (MQ=255)
                                                                                                                                                                            cGCTACAACAAAGTTATCTATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacccacc                                                                                                                           >  2:730870/1‑300 (MQ=255)
                                                                                                                                                                             gCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAACGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACgg                                                                                                                          <  1:112880/301‑1 (MQ=255)
                                                                                                                                                                              cTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATCCGTCAGCTTGCTGGTATGCGTGGTCTGCTGGCGAAGCCGGATGGCTCCCTCATCGAAACGCCCATCACCGCGCACTTCCGTCAAGGTCTTCCCGTCCCCCAGTACCTCCTCTCCACCCCCggc                                                                                                                         >  1:317361/1‑300 (MQ=255)
                                                                                                                                                                                acaacaAAGTTATCGATATCTGGGCGGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                                                                                                                                       <  2:165395/269‑1 (MQ=255)
                                                                                                                                                                                 caacaaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGGGcc                                                                                                                      >  2:132062/1‑300 (MQ=255)
                                                                                                                                                                                         ttATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGGACGTCCTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAgg                                                                                                              >  2:331461/1‑301 (MQ=255)
                                                                                                                                                                                          aaTCGGTATCGGGGCGGCGGCGAACGATAGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt                                                                                                             <  2:450582/300‑1 (MQ=255)
                                                                                                                                                                                            tCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCACCTCCACCCCCGGTCACCGTAAAGGTCt                                                                                                           >  1:299266/1‑301 (MQ=255)
                                                                                                                                                                                                tatCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCg                                                                                                       >  2:559824/1‑301 (MQ=255)
                                                                                                                                                                                                    ggggCTGCGGCGACCGATCGTGTATCCAAAGCGATGATAGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATa                                                                                                   <  2:704062/300‑1 (MQ=255)
                                                                                                                                                                                                     gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCCCGGTGCTCGTAAAGGTcagaccggaa                                                                                                   >  1:252072/1‑291 (MQ=255)
                                                                                                                                                                                                     gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTc                                                                                                            <  2:252152/291‑1 (MQ=255)
                                                                                                                                                                                                         tGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGGAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcaca                                                                                                                                                                                                                                                        <  1:154312/147‑1 (MQ=255)
                                                                                                                                                                                                         tGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcaca                                                                                                                                                                                                                                                        >  2:154359/1‑147 (MQ=255)
                                                                                                                                                                                                                         gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCGCACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCTTCTGCTTGCCGGTATGCGTGGTCTGATGGCGCAGCCGGATGGCTCCAGCATGGAAACGCCAATCACGGCGAACTTCCGTGACGGTCTTAACGGCATCCAGTCCTTCACAGCCACCCCCGGTGTCAGTAACGGACTGGCGGACACCGcccatacaaccgggaact                                                                              >  1:397097/1‑284 (MQ=255)
                                                                                                                                                                                                                         gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAAc                                                                               <  2:577359/300‑1 (MQ=255)
                                                                                                                                                                                                                           aTCCAAAGCGATGATGGATAACCGGCAAACGGAACCCGTGATTAACCGGGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTcc                                                                            <  2:542770/301‑1 (MQ=255)
                                                                                                                                                                                                                                aaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTa                                                                       <  2:317455/301‑1 (MQ=255)
                                                                                                                                                                                                                                          ggATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTTGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAACCACCGCGAACTTCCGTGAAGGTCTGAACGTACCCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGaaaccctctaaccccgcttccttgagcgcc                                                             >  1:650723/1‑274 (MQ=255)
                                                                                                                                                                                                                                            aTAACCGTCAAACTGAAACCGGGATTGACCGTGACGGTAAAGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTACATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtcgtcg                                                           <  1:559668/301‑1 (MQ=255)
                                                                                                                                                                                                                                               aCCTGCAAAGCAGAACCGGGATTAACCGTGACGGGCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCt                                                        <  1:360637/301‑1 (MQ=255)
                                                                                                                                                                                                                                                 cTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTACGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCCCGTAAAGGTCTGGCGGATCCCGCACTTAAAACGACCAACTCCCGTCACCGTCCTCGTGGTCTgg                                                      >  1:719363/1‑301 (MQ=255)
                                                                                                                                                                                                                                                           aaaCCGTGATTTACCCTGAGGGGCAGGAAGGAGAGCAGGTTTCCCTCCACGGCATCTACTGGGTGGCCGGCTCCGGTGCGGGTGGTTCTGCGGCACAGATTCGGCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgc                                            <  1:258335/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGCCGTCTGGTTGACGTGgcgc                                          >  2:348175/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                     ggagttcagacggtgctcttccgttctggTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGAcc                                    <  2:421321/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccggtcaggTCCGGAAGGGAGGCAGGTTTCCTTCCACAGCATCTACAGGATGGCGGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGCGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCtggtg                               <  2:89829/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ggTGGCGGTCAGGAAGCGAAGCAGGTTTCTTTCCACAGCATCTACAGGATGGCCGGCTCCGGTGCGCGTGGTTCTGGGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGAACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCtggtgg                              <  1:72842/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gcGGACAAGAAGGAAAGAAGGTTTCCTTCCTCGGCATCGCCCGGATGGCCGCCTCCGGTGCGGGTGGGTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTAc                          <  1:132018/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cGGTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGACACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCCGGTTGCCGTGGCGCGGGCCCTGGTGGTTAcc                         >  1:64629/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                               ggTCAGGAAGAGAAGCAGGTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGaaa                                                                                        <  2:612830/237‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   acgaaggcaaGCCGGTTTCCTGCCCCGGCCTCGACATGGTGGCCGCCTCCGGTGCGCTTGGTTCGGCGGCACAGGTTCGTCAGCTTGCTGGTATACGTGGTCTGATTGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCTCCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGa                    <  1:348070/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         agaagaAGGTTTCATTCCACAGCATCTACATGATGGCCGACGCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATtg              <  2:454141/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           aaGCAGGGTTCCTTCCAAAGCATCTACATGATTGCCGACTCCGGTGCGCGGGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCCCCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATtgtg            <  2:397221/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                ggTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGCTACCTGACTCGTCGTCTGGTTGCACTGGAGCTGGACCTGGTGGTTAGCGAAGACGACTGTgggccc       >  1:661435/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                                                ggTTTCCTTCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCCCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAACCCGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGCCGTGGCGCAGGAccagatcggaaaggcgtcgggtaggggaagg       >  1:421195/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                                                 gTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTGCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTAccc      <  2:674244/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       ttCCACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCAGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTTCCCATGa  >  2:256805/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       
TGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGGTCTGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGA  >  NC_000913/4187209‑4187802

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