Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 374,683 G→C 5.3% P272P (CCG→CCC) ‡ mhpE → 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,6830GC 5.3% 465.6 / 12.2 131P272P (CCG→CCC) ‡mhpE4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I
Reads supporting (aligned to +/- strand):  ref base G (62/62);  new base C (3/4);  total (65/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 132 total aligned reads displayed.

AAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATC  >  NC_000913/374420‑374955
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                  
aaaGCAGGGCTGAAACCAGAAACGGCAACTGGCAGGCACGCTCACCATAACCGGAGTCTTGGCGTGGCGAACGCTATGGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAAGGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGa                                                                                                                                                                                                                                                                         <  4:255094/275‑1 (MQ=255)
                   aaaCGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGGGTCGGTCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                      >  1:122638/1‑275 (MQ=255)
                      cGCAAACTGGCATGCACGCTCACCAGAACCTGAGTCTTGGCGGGGCGAAATCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTggg                                                                                                                                                                                                                                                    <  1:155836/274‑1 (MQ=255)
                                         tCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTTCGAGTCGATCGCGAAACGCTGGCGCTGGGATTCGAAGGTGTTTACTCg                                                                                                                                                                                                                                 >  2:352615/1‑274 (MQ=255)
                                          cACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGGCGATCGCGAAACGCTGGGTCGGGGATCCGCCGGGGTTTACTCGAg                                                                                                                                                                                                                               >  4:338042/1‑275 (MQ=255)
                                                    tGAGTCTTGGCGTGGCGACCTCCATCGGCGCGGTGGAAGGTGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGt                                                                                                                                                                                                                     <  2:208706/275‑1 (MQ=255)
                                                                                  cggTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTTAAACGGCGGCGGctcggttgggc                                                                                                                                                                                       >  4:319211/1‑268 (MQ=255)
                                                                                    gTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCCCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGACCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCtt                                                                                                                                                                                    >  2:508704/1‑275 (MQ=255)
                                                                                            gggCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTCAGTGCgg                                                                                                                                                                            >  2:278909/1‑275 (MQ=255)
                                                                                              gCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTCATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCCAGCTTCCTGCGTCACTGTGAAACGGCGGCGGGGCGTTATGGCTTAAGTGCGTTg                                                                                                                                                                          >  1:161012/1‑275 (MQ=255)
                                                                                               cTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTCATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTTAAACGGCTGCGGCGCGTTATGGCTTAAGTGCGGTg                                                                                                                                                                          >  2:280514/1‑274 (MQ=255)
                                                                                                  cGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGCAACGGCGGCGGCGCGTGATGGCTTCAGTGCGGTGGata                                                                                                                                                                      >  2:440911/1‑275 (MQ=255)
                                                                                                                   gCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGc                                                                                                                                                    >  2:236067/1‑275 (MQ=255)
                                                                                                                             aaTGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTc                                                                                                                                                                                                                                                                           <  6:494358/148‑1 (MQ=255)
                                                                                                                               tGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGa                                                                                                                                                                                                                                                                         <  5:844385/148‑1 (MQ=255)
                                                                                                                               tGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGAACGGNACGAGTCGATCGCGCAACGATGGCGCTGGGGAGCGCGGGTGTTTAGTCGAGCGTACTGCGTCACTGTGAACCGGNGGCGGGGCGGTATGGACTAGCTGCGGTGGGNNATCTGGTGGAGTCGGGGCaagggaggggg                                                                                                                                         >  3:78937/1‑266 (MQ=255)
                                                                                                                                       ggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                       >  8:579550/1‑142 (MQ=255)
                                                                                                                                       ggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                       <  7:579550/142‑1 (MQ=255)
                                                                                                                                            aGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                            <  2:437860/274‑1 (MQ=255)
                                                                                                                                            aGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                            <  2:338043/274‑1 (MQ=255)
                                                                                                                                                  cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaa                                                                                                                                                                                                                                                                  >  7:237212/1‑136 (MQ=255)
                                                                                                                                                  cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaa                                                                                                                                                                                                                                                                  <  8:237212/136‑1 (MQ=255)
                                                                                                                                                  cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                      >  8:592437/1‑148 (MQ=255)
                                                                                                                                                         cTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACg                                                                                                                                                                                                                                               >  6:818994/1‑148 (MQ=255)
                                                                                                                                                             aaGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGt                                                                                                                                                                                                                                                                                  <  7:1107520/109‑1 (MQ=255)
                                                                                                                                                             aaGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGt                                                                                                                                                                                                                                                                                  >  8:1107520/1‑109 (MQ=255)
                                                                                                                                                               ggtttgaTTGCCGCCGCGGATAAACTGGGCGGGCAGCATGGTACCGAGCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACg                                                                                                         <  1:3029/268‑1 (MQ=255)
                                                                                                                                                               gggtTTATGGCCGCCGCGGATAAACGGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGt                                                                                                        <  2:171450/273‑1 (MQ=255)
                                                                                                                                                                tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGAt                                                                                                                                                                                                                                                                        <  7:689358/116‑1 (MQ=255)
                                                                                                                                                                tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGAt                                                                                                                                                                                                                                                                        >  8:689358/1‑116 (MQ=255)
                                                                                                                                                                tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                           <  6:586863/145‑1 (MQ=255)
                                                                                                                                                                tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                           >  5:586863/1‑145 (MQ=255)
                                                                                                                                                                 gtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGt                                                                                                                                                                                                                                                                                  <  8:590892/105‑1 (MQ=255)
                                                                                                                                                                 gtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGt                                                                                                                                                                                                                                                                                  >  7:590892/1‑105 (MQ=255)
                                                                                                                                                                   ttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGgcg                                                                                                    <  2:518269/275‑1 (MQ=255)
                                                                                                                                                                       tGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                 <  8:421167/148‑1 (MQ=255)
                                                                                                                                                                       tGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                 >  5:1089860/1‑148 (MQ=255)
                                                                                                                                                                         ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                 >  6:566487/1‑146 (MQ=255)
                                                                                                                                                                         ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                 <  5:566487/146‑1 (MQ=255)
                                                                                                                                                                            ccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGTCCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGCTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTCTGGCTTTATGTGGGTGTATGTTCGCGTTTAGCTGGGCAAACGCCGTATGGTTGGCGGCCTGGAGGATCTGATCGTTGACGGTGCGCGGGGTCTg                                                                                           >  4:223560/1‑275 (MQ=255)
                                                                                                                                                                            ccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCtt                                                                                                                                                                                                                            <  5:818994/148‑1 (MQ=255)
                                                                                                                                                                               cgGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGgtgt                                                                                                                                                                                                                                        <  7:288560/133‑1 (MQ=255)
                                                                                                                                                                                   tAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTg                                                                                                                                                                                                                        >  7:821914/1‑145 (MQ=255)
                                                                                                                                                                                          ggCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtgtg                                                                                                                                                                                                              <  8:1125410/148‑1 (MQ=255)
                                                                                                                                                                                          ggCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGAATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGGTCTGCGCAACAACAaata                                                                             >  2:100793/1‑275 (MQ=255)
                                                                                                                                                                                               gcagcaTGGGACCGATCTCTATGCGGTAAGGGAGGGTGCCGACGACCTGGGGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGGGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatga                                                                         <  4:287361/274‑1 (MQ=255)
                                                                                                                                                                                                  gcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAg                                                                                                                                                                                                                                                                             >  8:885494/1‑77 (MQ=255)
                                                                                                                                                                                                  gcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAg                                                                                                                                                                                                                                                                             <  7:885494/77‑1 (MQ=255)
                                                                                                                                                                                                  gcaGGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGGCCTGGTGCGGCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGAc                                                                     <  1:78937/275‑1 (MQ=255)
                                                                                                                                                                                                      gggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                      <  7:232966/96‑1 (MQ=255)
                                                                                                                                                                                                      gggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctg                                                                                                                                                                                                                                                      >  8:232966/1‑96 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACg                                                                                                                                                                                                                                                                               <  7:786645/70‑1 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACg                                                                                                                                                                                                                                                                               >  8:786645/1‑70 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcg                                                                                                                                                                                                                                                                      <  6:1100924/79‑1 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcg                                                                                                                                                                                                                                                                      >  5:1100924/1‑79 (MQ=255)
                                                                                                                                                                                                       ggACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCt                                                                                                                                                                                                                                             <  7:1137631/104‑1 (MQ=255)
                                                                                                                                                                                                         aCCTATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTg                                                                                                       >  2:269535/1‑234 (MQ=255)
                                                                                                                                                                                                          ccGATCTCTATGCGTTAATGGAGGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCgag                                                              <  2:433861/274‑1 (MQ=255)
                                                                                                                                                                                                            gATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACg                                                                                                                                                                                                                                               <  7:493607/97‑1 (MQ=255)
                                                                                                                                                                                                              tctcTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                       >  5:978093/1‑71 (MQ=37)
                                                                                                                                                                                                              tctcTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                       <  6:978093/71‑1 (MQ=38)
                                                                                                                                                                                                               ctctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                                                                                                                                                             >  6:181461/1‑80 (MQ=255)
                                                                                                                                                                                                               ctctATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                                                                                                                                                             <  5:181461/80‑1 (MQ=255)
                                                                                                                                                                                                                 ctATGCGTTAATGGATGCCGCCGACGCCCGGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCa                                                      <  4:100792/275‑1 (MQ=255)
                                                                                                                                                                                                                  tATGCGTTAATGGAGGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCa                                                      <  2:319212/274‑1 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctgg                                                                                                                                                                                                                                                           >  8:887028/1‑78 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctgg                                                                                                                                                                                                                                                           <  7:887028/78‑1 (MQ=255)
                                                                                                                                                                                                                   aTGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAg                                                                                                                                                                                                                               <  7:651764/106‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtg                                                                                                                                                                                                                <  5:928766/119‑1 (MQ=255)
                                                                                                                                                                                                                         tAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGCTGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCTGAGACCGTGCATTTTGTa                                               >  2:527644/1‑274 (MQ=255)
                                                                                                                                                                                                                          aaTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGt                                                                                                                                                                                                                                                                            <  8:164148/54‑1 (MQ=255)
                                                                                                                                                                                                                           aTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                 <  6:1098468/96‑1 (MQ=255)
                                                                                                                                                                                                                                 gccgccGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGata                                                                                                                                                                      >  5:1024389/1‑148 (MQ=255)
                                                                                                                                                                                                                                  ccgccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                                 <  8:283312/89‑1 (MQ=255)
                                                                                                                                                                                                                                       gacgagCTGGTGCGTACGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtg                                                                                                                                                                                                                <  8:1034686/101‑1 (MQ=255)
                                                                                                                                                                                                                                       gacgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTTACtg                                                                                                                                                                                                                >  7:1034686/1‑101 (MQ=255)
                                                                                                                                                                                                                                       gacgacCGGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCa                                 <  2:78820/274‑1 (MQ=255)
                                                                                                                                                                                                                                           acCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGa                                                                                                                                                            <  6:1024389/148‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGaaa                                                                                                                                                                                                           >  7:328480/1‑99 (MQ=255)
                                                                                                                                                                                                                                              tGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGaaa                                                                                                                                                                                                           <  8:328480/99‑1 (MQ=255)
                                                                                                                                                                                                                                               ggTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTTCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATGCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGCTATGATCGTTTACGTGCCGCTGGCTCTGCGCAACAACAAATCATGCTTACTCCCCAGAACGTGCATTTTGTTTGCCCATCCCAgacccaacccc                         >  4:341489/1‑264 (MQ=255)
                                                                                                                                                                                                                                                 tGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTcc                                                                                                                                                                                                                          >  7:225231/1‑81 (MQ=255)
                                                                                                                                                                                                                                                 tGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTggg                                                                                                                                                     >  7:882983/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCaaa                                                                                                                                                 <  8:882983/148‑1 (MQ=255)
                                                                                                                                                                                                                                                      ccGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCgg                                                                                                                                                                            >  7:147513/1‑121 (MQ=255)
                                                                                                                                                                                                                                                      ccGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCgg                                                                                                                                                                            <  8:147513/121‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTgg                                                                                                                                                      >  5:229661/1‑138 (MQ=255)
                                                                                                                                                                                                                                                               gATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                           >  5:592096/1‑50 (MQ=255)
                                                                                                                                                                                                                                                               gATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTgg                                                                                                                                                                                                                                           <  6:592096/50‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACt                                                                                                                                                                                                                                   >  8:817670/1‑56 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                                                                                                                                                          >  5:775160/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                 tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTAt                                                                                                                                                                                          <  6:775160/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCCGGGCAAACGCCGGATGGTTgg                                                                                                                                   <  6:884385/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aCCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTAt   >  2:121397/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     cgggTGCGAGTGGGTCGGGAAAGGGGGGGGCTGGGGTAGGCTGGGGGTTACGCGGGGTTCTTGAGGCAATGGGAAGCGGAGGCGGGGCGTTATGGCGTAAGGGCGGTGGGTATTCTCGTTGGGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGGTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAAGAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTAt   <  4:286551/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTGCGGGTCGTGCGCGAAACGCTGGGGCTGGGATACGCTGGGGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTGAGGGTTTAAGGGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGGGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATc  <  4:440910/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGc                                                                                                                                                                                                                              <  8:311311/57‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGc                                                                                                                                                                                                                              >  7:311311/1‑57 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGtt                                                                                                                                     <  8:146588/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGtt                                                                                                                                     >  7:146588/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                  
AAAGCAGAGCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATC  >  NC_000913/374420‑374955

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