Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 986,411 G→T 5.3% S191Y (TCT→TAT) ‡ ompF ← outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,4110GT 5.3% 732.5 / 21.3 206S191Y (TCT→TAT) ‡ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base G (97/97);  new base T (6/5);  total (104/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 208 total aligned reads displayed.

ACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGT  >  NC_000913/986138‑986681
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
             cTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGACCGGAGTGCAgtg                                                                                                                                                                                                                                                                  >  1:149682/1‑275 (MQ=255)
              ttGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCCGCTGCCCCATAAGCCCCAACGATACCAACGCCTTCGTCTTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAttg                                                                                                                                                                                                                                                                  >  2:515633/1‑274 (MQ=255)
               tGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCTCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATCCCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCCACCACGTCCCCGTTAGAACGGCGTGCAgtgt                                                                                                                                                                                                                                                                 >  4:171417/1‑274 (MQ=255)
                           gCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTCGCTGATAGAACCGCCACCACCGTCGCCGTTAGCACGGCGTTCAGGTTCCGGCTCGttttt                                                                                                                                                                                                                                                    >  4:335691/1‑275 (MQ=255)
                                                aaTTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCg                                                                                                                                                                                                                                 <  2:313195/273‑1 (MQ=255)
                                                aaTTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                                                                                                                                                                                                <  2:33489/274‑1 (MQ=255)
                                                               gtgttAGCGTTACGGGTTGCACCGTAGTTGGCGTCCAGGGAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGATGTTGAGCTTCTTGCAGGTTGGTACGGTAAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                                                                                                                                                                                                <  4:44004/271‑1 (MQ=255)
                                                                                  cACCGTAGTTCGCTGCCAGGTAGATGGTTTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAg                                                                                                                                                                                             <  2:147214/275‑1 (MQ=255)
                                                                                         gTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCTCCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt                                                                                                                                                                                      <  2:437700/275‑1 (MQ=255)
                                                                                                           gttgttCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTACAGTGTCACGCTCGTTTTTACACAGGTACTGAACAGCCAAGTTCAGGCCATCAACAGGACAAAAGAGGTTGGAGTTACGAAAGGTAGCAACGCAGCCaa                                                                                                                                                                    >  2:271143/1‑275 (MQ=255)
                                                                                                                               gACCAGTAGCCCACTGTTCAGATGTTTTACAGTTGACAAGAGGTTGAGCTTCGTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTGCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCa                                                                                                                                                <  1:84501/275‑1 (MQ=255)
                                                                                                                                aCCAGTAGCTCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTTGTACGGCCAGCTGCACCATAAGCAGCAACGATACCAAAGCCTCCGTATTCGTAGCCGATAGAACCGCCTACACCGTCACCGTTACAACGTGCTGCAGTGTTAAGCTCGTTTTTACCCCGGTACTTAACAACGAAGTTGAGGCCATCACCCATACCTAATAAGTTGGAGTTCCGATAGGTAGCCCACCACCCACCACGACCCACGAACAAGTcc                                                                                                                                                >  4:84384/1‑273 (MQ=255)
                                                                                                                                  cAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCACCAACACGACCAACGAAGAAGTCATCg                                                                                                                                             <  2:244394/275‑1 (MQ=255)
                                                                                                                                    gtggCCCACTGTTCAGCTTTTTTACCTTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACCGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCATATTCGAAGCTGATAGAACCGCCAGAACCGTCGCCGTTAGAACTGCGTGCAGTGTCACGCACGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTCGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCt                                                                                                                                           <  4:267804/272‑1 (MQ=255)
                                                                                                                                    gTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCt                                                                                                                                           <  2:92124/275‑1 (MQ=255)
                                                                                                                                       gCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTa                                                                                                                                        <  2:171418/275‑1 (MQ=255)
                                                                                                                                        cccACTGTTAAGTTTTTTTCCCGTGGCCATGATGTGGAGCTTTTTGAAGGTTGGTACGGTCAGCTGAACCATAAGCACCAAAGATACCCAAGCCTTCGTATTCGGAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGGGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTAt                                                                                                                                       <  4:515632/275‑1 (MQ=255)
                                                                                                                                            cTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACgg                                                                                                                                                                                                                                                                           >  7:489897/1‑139 (MQ=255)
                                                                                                                                                    cTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgtgt                                                                                                                                                                                                                                                                 >  6:326332/1‑141 (MQ=255)
                                                                                                                                                      tttttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGt                                                                                                                                                  <  2:472378/250‑1 (MQ=255)
                                                                                                                                                         ttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCaaa                                                                                                                      >  4:79185/1‑275 (MQ=255)
                                                                                                                                                         ttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAAAGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCACCCAGACCAAAGAAGTTGGCGTTAAGATAGGTACCAACTCCGCCACACCGACAAACGAAGAAGTAATGGATGTATGCAGGATcaccacacaa                                                                                                                      >  4:351956/1‑271 (MQ=255)
                                                                                                                                                             cGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                                                                                                                                                                                                 <  8:506436/148‑1 (MQ=255)
                                                                                                                                                             cGTTGCCAAGAGGTTGAGCTTCTGGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCt                                                                                                                                           <  1:103853/250‑1 (MQ=255)
                                                                                                                                                               tggCCAAGAGGTTGGGCTTCTTGCAGGTTGGTACGGTAAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTgg                                                                                                                <  2:66131/273‑1 (MQ=255)
                                                                                                                                                                    aaGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTCTAACgg                                                                                                                                                                                                                                                                           >  5:1119842/1‑115 (MQ=255)
                                                                                                                                                                    aaGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGc                                                                                                                                                                                                                                                                      >  8:581/1‑120 (MQ=255)
                                                                                                                                                                    aaGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAgg                                                                                                                                                                                                                                             <  6:169947/145‑1 (MQ=255)
                                                                                                                                                                       agGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                                                        <  2:377377/275‑1 (MQ=255)
                                                                                                                                                                           tGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCACCGTTCTAACg                                                                                                                                                                                                                                                                            >  6:756844/1‑107 (MQ=255)
                                                                                                                                                                                    tGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGt                                                                                                                                                                                                                                                        <  7:1119539/118‑1 (MQ=255)
                                                                                                                                                                                    tGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTa                                                                                                                                                                                                                                                   <  7:884344/123‑1 (MQ=255)
                                                                                                                                                                                           tGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                                                                                                                                                                                                   <  6:352672/148‑1 (MQ=255)
                                                                                                                                                                                              tACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCg                                                                                                                                                                                                                                                                         <  6:722444/91‑1 (MQ=255)
                                                                                                                                                                                               aCGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCt                                                                                                                                                                                                                                                           <  5:185912/104‑1 (MQ=255)
                                                                                                                                                                                               aCGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                                                                                                                                                                                                     >  6:821136/1‑126 (MQ=255)
                                                                                                                                                                                               aCGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAAc                                                                                                                                                                                                                                     <  5:821136/126‑1 (MQ=255)
                                                                                                                                                                                                  gTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                                                                                                                                                                            <  5:951982/148‑1 (MQ=255)
                                                                                                                                                                                                   tCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGa                                                                                                                                                                                                                                                                               >  5:473832/1‑80 (MQ=255)
                                                                                                                                                                                                   tCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGa                                                                                                                                                                                                                                                                               <  6:473832/80‑1 (MQ=255)
                                                                                                                                                                                                   tCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                                                                                                                                                                                           >  8:951333/1‑148 (MQ=255)
                                                                                                                                                                                                       cTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCg                                                                                                                                                                                                                                                         >  5:317586/1‑98 (MQ=255)
                                                                                                                                                                                                             cATAAGCACCAACTGTACCAATGACGTGGAATGCGTGGCGGGTTGAACCGTCAACGCCGTCGCCGGTGGACCGCGGGGCAGTGTGACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCATACCAAAGAAATTGGGGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                                                  <  4:362175/275‑1 (MQ=255)
                                                                                                                                                                                                             cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCTGGTACTGAACAGCGa                                                                                                                                                                                                                                >  5:331173/1‑117 (MQ=255)
                                                                                                                                                                                                             cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                                                                                                                                                                                                >  6:869922/1‑117 (MQ=255)
                                                                                                                                                                                                              aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAg                                                                                                                                                                                                                                   <  6:1048436/113‑1 (MQ=255)
                                                                                                                                                                                                                  gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGc                                                                                                                                                                                                                                                                      <  7:1050331/74‑1 (MQ=255)
                                                                                                                                                                                                                    aCCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                                           >  2:293091/1‑275 (MQ=255)
                                                                                                                                                                                                                     ccAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGa                                                                                                                                                                                                                                       >  8:346158/1‑102 (MQ=255)
                                                                                                                                                                                                                      cAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtt                                                                                                                                                                                                                                                       <  7:328595/85‑1 (MQ=255)
                                                                                                                                                                                                                      cAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                                                                                                                                                                                                     <  5:1031584/135‑1 (MQ=255)
                                                                                                                                                                                                                       aaCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                                                                                                                                                                                                 <  8:297752/90‑1 (MQ=255)
                                                                                                                                                                                                                       aaCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCAATCAACCAGACCAAAGACGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGACGAAGTCATCGATGTATGCAGTATCCCCACCAAATTCTGGCAGCATATCGTTGTAACCCAGTGCATCCACAACCACACCGAAATTACCGCCGTaaa                                                        >  1:99499/1‑274 (MQ=255)
                                                                                                                                                                                                                        aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                                                                                                                                                                                                 >  6:591228/1‑121 (MQ=255)
                                                                                                                                                                                                                        aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                                                                                                                                                                                                 <  5:591228/121‑1 (MQ=255)
                                                                                                                                                                                                                         cGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                                                                                                                                                                                             <  7:754369/124‑1 (MQ=255)
                                                                                                                                                                                                                                aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgtgt                                                                                                                                                                                                                                                                 >  8:1077372/1‑65 (MQ=255)
                                                                                                                                                                                                                                aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                                                                                                                                                                                  >  8:46359/1‑128 (MQ=255)
                                                                                                                                                                                                                                 aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGttt                                                                                                                                                                                                                                                      <  7:669364/75‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCg                                                                                                                                                                                                                                                                         >  7:1016562/1‑55 (MQ=255)
                                                                                                                                                                                                                                  aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgt                                                                                                                                                                                                                                                                   <  6:551765/61‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgt                                                                                                                                                                                                                                                                   >  5:551765/1‑61 (MQ=255)
                                                                                                                                                                                                                                  aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                                                                                                                                                                                               >  6:290763/1‑129 (MQ=255)
                                                                                                                                                                                                                                  aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                                             >  1:122008/1‑275 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaa                                                                                                                                                                                                        >  5:154631/1‑118 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaa                                                                                                                                                                                                        <  6:154631/118‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                                                                                                                                                                                  <  7:168362/124‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                                                                                                                                                                                  >  8:168362/1‑124 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                                                                                                                                                                              <  7:281206/128‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                                                                                                                                                                              >  8:281206/1‑128 (MQ=255)
                                                                                                                                                                                                                                    ccTTCGTATTCGCAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTCCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAACTAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCCACGAAGAATTCATCCCCTTATGCCGTATAACCCCCAAATTCTCGCCGCATACCGGTGTAACCCAGATCATCATAAACCACACCGTAGTTACGGCAGTAATCGAAAGaaacgaa                                           >  4:43098/1‑270 (MQ=255)
                                                                                                                                                                                                                                      ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCCCGCTCGCTTTTTCCCCGGTCCTGAACAGCGAAGGTCAGGCCATCAACAAGACCAAAGAAGTCGGAGATACGATGAGTATCAACGCCGCCGAAACGACCAACGAAGACGCCATCGACGCAGCCAGTGTCCACACCAAATTCTGGCAGCATATCAGCGTAAACCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCACAAGAACAAActt                                         >  4:482486/1‑273 (MQ=255)
                                                                                                                                                                                                                                      ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                                                                                                                                                                                                               >  8:658436/1‑109 (MQ=255)
                                                                                                                                                                                                                                      ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                                                                                                                                                                                  >  7:349086/1‑122 (MQ=255)
                                                                                                                                                                                                                                       tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAc                                                                                                                                                                                                                <  5:873639/107‑1 (MQ=255)
                                                                                                                                                                                                                                       tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccg                                                                                                                                                                        >  6:1105271/1‑147 (MQ=255)
                                                                                                                                                                                                                                           attcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                                                                                                                                                                                               >  7:551936/1‑120 (MQ=255)
                                                                                                                                                                                                                                                taGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGAc                                                                                                                                                              >  5:755856/1‑148 (MQ=255)
                                                                                                                                                                                                                                                  gCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAac                                                                                                                                                                   >  8:841769/1‑141 (MQ=255)
                                                                                                                                                                                                                                                      aTAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAg                                                                                                                                                                                    >  5:1096583/1‑120 (MQ=255)
                                                                                                                                                                                                                                                           aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                                                                                                                                                                                                 <  5:212867/86‑1 (MQ=255)
                                                                                                                                                                                                                                                           aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                                                                                                                                                                                                 >  6:212867/1‑86 (MQ=255)
                                                                                                                                                                                                                                                           aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                                                                                                                                                                                               >  7:834424/1‑104 (MQ=255)
                                                                                                                                                                                                                                                            ccGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACgaaga                                                                                                                                                     <  7:159332/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg                                                                                                                                                                                        <  7:999792/106‑1 (MQ=255)
                                                                                                                                                                                                                                                                cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACga                                                                                                                                                        <  7:592550/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACgaaga                                                                                                                                                     >  6:948618/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                     ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                                                                                                                                                                                                >  8:1023404/1‑61 (MQ=255)
                                                                                                                                                                                                                                                                     ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGa                                                                                                                                                                                                                                <  5:283867/61‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cGTCGCCGTTAGAACGGCGTGAAGTGGCACACTCGTTTTTACCCAGGTACTGAGCGGCGAAGTTCAGGCCATCCACCGGACCAAGGAAGTTGGAGTGACGCTAGGTAGCAACGCCGCCAACACGACCAACGAAAAAGTCATCGCTATATGCAGTATAACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAATAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGAcgc         <  4:227688/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                                                                                                       >  4:200303/1‑180 (MQ=255)
                                                                                                                                                                                                                                                                        tCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                                                                                                                                                                                               >  5:706910/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                        acGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccg                                                                                                                                                                        >  6:704157/2‑114 (MQ=255)
                                                                                                                                                                                                                                                                        acGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCa                                                                                                                                                            >  5:639314/2‑126 (MQ=255)
                                                                                                                                                                                                                                                                         cGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACga                                                                                                                                                        >  8:1119893/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                          gCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCATGCCATCAACCAGACCAAAGAAGtt                                                                                                                                                                                                 >  8:778260/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                          gCCGTTAGAACGAGTTGCGATGTAACGCTNNNTTTTACCCCGGTACTGAACAGCGAAGTTCAGGCCATCAACCAAACCAAAGAAGTTGAGGTTACGTTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATGCTGGCAGCCTATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGttt     <  3:149682/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ccGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                                                   <  2:200304/180‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGt  >  2:491303/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                               tAGAACGGCGTGCAGTGTCACGCTCGTTTTTAGCCAGGTACTGAACATCGTAGTTCAGGCCATCAACCAGACCAaaga                                                                                                                                                                                                     <  8:701095/78‑1 (MQ=38)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
ACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGATCGGCGTAGCGTTACGGGTTTCACCGTAGTTCGCTGCCAGGTAGATGTTGTTCGCGTCGTACTTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGT  >  NC_000913/986138‑986681

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