Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 374,689 A→G 7.4% R274R (CGA→CGG) ‡ 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,6890AG 7.4% 432.3 / 17.2 143R274R (CGA→CGG) ‡mhpE4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I
Reads supporting (aligned to +/- strand):  ref base A (63/68);  new base G (5/6);  total (69/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-02

Only 100 of 143 total aligned reads displayed.

GCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATC  >  NC_000913/374428‑374955
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                           
gCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaaa                                                                                                                                                                                                                                                                >  2:678056/1‑275 (MQ=255)
     acagaGAAGCGCAAACTAGCATGCACGCTCACCATAACCTGAGGCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAGACTGGGCTGGCAGCATGGGACCGTTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgct                                                                                                                                                                                                                                                            <  2:83567/270‑1 (MQ=255)
                                 tCACCATAACCTGAGTCTTGGCGTGGTGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGCTAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCTGCCGACGACCTGGTTCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGCTTACTCg                                                                                                                                                                                                                               >  2:174823/1‑274 (MQ=255)
                                               gTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCCTGGGACCGATCTCTATGCGTTACTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGCTCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGGTCCGCt                                                                                                                                                                                                                                           >  2:690142/1‑248 (MQ=255)
                                               gTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGCTCGACCGGTACGAGTCGAACGCGAAACGCTGGCGCTGGGATACCCGGGTGGTTACTCGGACTTCCGGCGTctc                                                                                                                                                                                                                >  4:147405/1‑273 (MQ=255)
                                               gTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGACGACGACCTGGTGCGTCAGTTGNCGGAACGAACGGTCCGCGTCCAACGCGAAACGCAGGCGCTGGGATACGCGGGGGTTTACTCGAGAttacgtagtccc                                                                                                                                                                                                                >  3:162265/1‑265 (MQ=255)
                                                 cTTGGCGTGGCGAACCCTATCGGGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGGGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACtg                                                                                                                                                                                                              <  2:632541/275‑1 (MQ=255)
                                                        tggcgAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCAACGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATCCCGCCGACGACCTGGTGCGTCCGTTGCAGGATTGACCGGTTAGGGCCGATCGCGTAACGCAGGCGCTGGGAGACGCTGGTGTTTACGCTAGCTTCCTGCGTCACTGTAAAAcg                                                                                                                                                                                                       >  4:510727/1‑275 (MQ=255)
                                                                            gTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCGGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTANGCGGTAATGGATGCCGCCGCCGACCTGGTGCGTCCGTTGGAGGAGCGACCGGGACGGCTCGATCGCGAANCGCTGGGGCGGGGGATAGCTGGGGGTCACTGGGGCATGCGGCGTCNCNGTGCAAGGgcggcggcggggtgtgggct                                                                                                                                                                                   >  3:79529/1‑264 (MQ=255)
                                                                                       cTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCTAGCTTCCTTCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGGGGTgg                                                                                                                                                                        >  2:243977/1‑275 (MQ=255)
                                                                                               gAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtctc                                                                                                                                                                >  2:370221/1‑275 (MQ=255)
                                                                                                                gcgGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACg                                                                                                                                               <  2:510729/275‑1 (MQ=255)
                                                                                                                     aaTGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTc                                                                                                                                                                                                                                                                          <  5:64048/148‑1 (MQ=255)
                                                                                                                      aTGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                               >  4:296576/1‑189 (MQ=255)
                                                                                                                        gggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGAt                                                                                                                                                                                                                                                                       <  5:569763/148‑1 (MQ=255)
                                                                                                                        gggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATc                                                                                                                                                                                                                                                                      <  7:285661/149‑1 (MQ=255)
                                                                                                                         ggcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATcc                                                                                                                                                                                                                                                                     <  7:840780/149‑2 (MQ=255)
                                                                                                                          gcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGATcct                                                                                                                                                                                                                                                                    <  8:140074/149‑3 (MQ=255)
                                                                                                                          gcgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCt                                                                                                                                                                                                                           <  2:296577/189‑1 (MQ=255)
                                                                                                                           cgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcgcg                                                                                                                                                                                                                                                                   <  7:706359/149‑1 (MQ=255)
                                                                                                                           cgcgggcgcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATcgc                                                                                                                                                                                                                                                                    >  6:656242/1‑148 (MQ=255)
                                                                                                                           cgcgggcgaCGGTAAGCCCCCGGTGGAGATGTTTACTGCGGCGGCGGATAAACTGGGCTGGCAGCATGGGACCGAGCTCTATGCGTTAATGGAGTCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGGGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTg                                                                                                                                    <  4:189757/275‑1 (MQ=255)
                                                                                                                                  gcAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgc                                                                                                                                                                                                                                                             >  6:555704/1‑148 (MQ=255)
                                                                                                                                    aGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                            <  2:147406/274‑1 (MQ=255)
                                                                                                                                    aGGTAACGAACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                            <  2:11847/274‑1 (MQ=255)
                                                                                                                                       tAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcg                                                                                                                                                                                                                                                        <  5:778799/148‑1 (MQ=255)
                                                                                                                                          cGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgctgg                                                                                                                                                                                                                                                   <  8:231142/149‑1 (MQ=255)
                                                                                                                                           gCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTggg                                                                                                                                                                                                                                                  <  8:492722/149‑1 (MQ=255)
                                                                                                                                                gCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATCAGCCGGTGTTTACTCTAGCTTACTGCGTCCATGTGAAACGGGGGCGGCGCGTTATGGCTTAAGt                                                                                                                                                                               >  4:37387/1‑211 (MQ=255)
                                                                                                                                                gCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCg                                                                                                                                                                                                                    >  4:345997/1‑174 (MQ=255)
                                                                                                                                                gCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGGGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGAttg                                                                                                               >  2:410985/1‑273 (MQ=255)
                                                                                                                                                     aaGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGgt                                                                                                                                                                                                                                        <  7:226412/149‑1 (MQ=255)
                                                                                                                                                       gtgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCTATCGCTATGCGGTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCTATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCGGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGt                                                                                                        <  1:135587/275‑1 (MQ=255)
                                                                                                                                                       gtgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGt                                                                                                        <  2:652354/275‑1 (MQ=255)
                                                                                                                                                        tgtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGttt                                                                                                                                                                                                                                    >  8:476482/1‑149 (MQ=255)
                                                                                                                                                         gtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGCTATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGCTTTGAACGTTGCCGtt                                                                                                       >  2:200758/1‑273 (MQ=255)
                                                                                                                                                         gtTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTgg                                                                                                      <  2:655053/275‑1 (MQ=255)
                                                                                                                                                           ttATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTAc                                                                                                                                                                                                                                  <  7:328576/149‑1 (MQ=255)
                                                                                                                                                                 ccgccgCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAg                                                                                                                                                                                                                             >  5:651120/1‑148 (MQ=255)
                                                                                                                                                                  cggagcgGTTAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCCGGTGCGGCCGTTGGAGGATCGACCTGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGAt                                                                                              <  4:678054/270‑1 (MQ=255)
                                                                                                                                                                    ccgcGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTTATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTCCGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGCGTTTACTCGAACTTCCCGCGTCACCGTGACACGGCGGCGGCGCGTTCTGGCTTCAGTCCGGGAGGTGTTCACGTTGAGCTGGGCAAACCACGGCGGGGTGGCGGCAGGGAGGAAGTGGTCGTTTGAGGGGCGCTgggtgtg                                                                                           >  4:592115/1‑270 (MQ=255)
                                                                                                                                                                          aTAACCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGTATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGcaac                                                                                     <  1:162266/275‑1 (MQ=255)
                                                                                                                                                                            aaaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgc                                                                                                                                                                                                                                                             <  6:761107/106‑1 (MQ=255)
                                                                                                                                                                            aaaCTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgc                                                                                                                                                                                                                                                             >  5:761107/1‑106 (MQ=255)
                                                                                                                                                                                   gCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGCCCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACGGTGa                                                                                                                                                                                                           >  5:665611/1‑148 (MQ=255)
                                                                                                                                                                                        cagcaTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGc                                                                                                                                                                                                                                            <  6:686754/110‑1 (MQ=255)
                                                                                                                                                                                           caTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCAc                                                                                                                                                                                                                <  6:159249/135‑1 (MQ=255)
                                                                                                                                                                                           caTGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCAc                                                                                                                                                                                                                >  5:159249/1‑135 (MQ=255)
                                                                                                                                                                                                gACCGATCTCTATGCGTTAATGGATGTCGCCGACGACCTTGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCg                                                                <  2:625805/274‑1 (MQ=255)
                                                                                                                                                                                                 aCCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaaa                                                                                                                                                                                                                                                                >  6:381538/1‑82 (MQ=255)
                                                                                                                                                                                                 aCCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGaaa                                                                                                                                                                                                                                                                <  5:381538/82‑1 (MQ=255)
                                                                                                                                                                                                 aCCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgct                                                                                                                                                                                                                                                     >  4:91973/1‑92 (MQ=255)
                                                                                                                                                                                                  ccGATCTCGATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGTTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCgag                                                              <  1:74387/274‑1 (MQ=255)
                                                                                                                                                                                                    gATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGCGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTa                                                                                                                                                                                         <  7:683915/149‑1 (MQ=255)
                                                                                                                                                                                                    gATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGtt                                                                                                                                                                                          >  5:793499/1‑148 (MQ=255)
                                                                                                                                                                                                     aTCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGa                                                                                                                                                                                                                                                 <  2:91974/92‑1 (MQ=255)
                                                                                                                                                                                                          tATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCtt                                                                                                                                                                                   >  7:198284/1‑149 (MQ=255)
                                                                                                                                                                                                              cGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACg                                                                                                                                                                                                                                             <  2:742575/87‑1 (MQ=255)
                                                                                                                                                                                                              cGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAg                                                                                                                                                                                >  5:862739/1‑148 (MQ=255)
                                                                                                                                                                                                               gTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACTCGTACCGGTCGa                                                                                                                                                                                                                                                                        >  7:473370/1‑60 (MQ=37)
                                                                                                                                                                                                                ttAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGgcgc                                                                                                   <  1:228179/223‑1 (MQ=255)
                                                                                                                                                                                                                  aaTGGATGTCGCCGACGACCTGGTTCGTCCGTTGCAGGATCGACCGGTACGAGCCGATCGCGAAACGCTTGCGCTGGGATACGCCGGTGTTTACTAGAGCTCACTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGTCTTCAGTGCGGTGGATATTCTCGCTGAGCTGGGCAAACTCCGGATGCTTGTCGTCCAGGCGGATTTGATCGTCTACGCGGCGTTGGACCTGCCCACATCCAAAAAATGAATACTTCCTGGAACGTGCATTTTTTAt                                              >  4:557707/1‑274 (MQ=255)
                                                                                                                                                                                                                  aaTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCg                                                                                                                                                                                                                               >  2:587242/1‑97 (MQ=255)
                                                                                                                                                                                                                  aaTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATAtt                                                                                                                                                                   >  1:96529/1‑157 (MQ=255)
                                                                                                                                                                                                                  aaTGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTggcgg                                                                                                                                >  4:660780/1‑192 (MQ=255)
                                                                                                                                                                                                                    tGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGa                                                                                                                                                                                                                                                 >  4:352035/1‑77 (MQ=255)
                                                                                                                                                                                                                    tGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAg                                                                                                                                                                                >  4:799067/1‑142 (MQ=255)
                                                                                                                                                                                                                      gaggccgccGACGACCTGNTGCGTCCGTTGCAGGTGCGACCGGTACGAGTCGTTCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCg                                                                                                                                                               <  3:96528/154‑1 (MQ=255)
                                                                                                                                                                                                                      gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCt                                                                                                                                                                                                                           <  4:587240/97‑1 (MQ=255)
                                                                                                                                                                                                                      gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                            <  2:660782/192‑1 (MQ=255)
                                                                                                                                                                                                                      gATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGcc                                           <  2:483890/273‑1 (MQ=255)
                                                                                                                                                                                                                        tGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACg                                                                                                                                                                                                                                             <  2:352037/77‑1 (MQ=255)
                                                                                                                                                                                                                        tGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAACTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTTACGTGGCGCGGGAACTGCGCCACAACCAAAAATGATTACTGCCCAGAACGGGCCATTTGTATGcccca                                        >  4:445657/1‑272 (MQ=255)
                                                                                                                                                                                                                           cgccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATAtt                                                                                                                                                                   <  5:66308/148‑1 (MQ=255)
                                                                                                                                                                                                                            gccgACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATtc                                                                                                                                                                  >  6:435025/1‑148 (MQ=255)
                                                                                                                                                                                                                                acgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGtt                                                                                                                                                             <  8:734345/149‑1 (MQ=255)
                                                                                                                                                                                                                                acgacCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGt                                                                                                                                                              >  5:22212/1‑148 (MQ=255)
                                                                                                                                                                                                                                    ccTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGTGGTGGATATTCTCGtt                                                                                                                                                             <  8:715259/145‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGTGGTGGATATTCTCGtt                                                                                                                                                             >  7:715259/1‑145 (MQ=255)
                                                                                                                                                                                                                                     cTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTCAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCt                                                                                                                                                        >  7:371/1‑149 (MQ=255)
                                                                                                                                                                                                                                             tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctg                                                                                                                                                                                                                                                           >  4:772855/1‑43 (MQ=255)
                                                                                                                                                                                                                                             tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAAc                                                                                                                                                >  7:674612/1‑149 (MQ=255)
                                                                                                                                                                                                                                             tCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAACTGCGGTGGATATTCTCGTTGAGCTGGGCAAAc                                                                                                                                                <  8:622691/149‑1 (MQ=255)
                                                                                                                                                                                                                                              ccGTTGCAGTATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCGGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACg                                                                                                                                               >  8:288461/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 ttGCAGGATCGACCGGTACGAGTCGATCGCGAAAcgctggcgc                                                                                                                                                                                                                                                       <  2:772857/43‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCag                                                                                                                            <  4:732099/165‑1 (MQ=255)
                                                                                                                                                                                                                                                     gcgATCGACTCGTACCGGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGtt                                                                                                                                                             <  8:579131/126‑1 (MQ=255)
                                                                                                                                                                                                                                                        aTCGACTCGTACCGGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGt                                                                                                                                      >  5:464340/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        aTCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGtt                                                                                                                                     >  7:672647/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCagg                                                                                                                           >  1:76633/1‑158 (MQ=255)
                                                                                                                                                                                                                                                          cGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTg                                                                                                                                    <  5:435025/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATAtt                                                                                                                                                                   <  6:380452/115‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATAtt                                                                                                                                                                   >  5:380452/1‑115 (MQ=255)
                                                                                                                                                                                                                                                             ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTggcg                                                                                                                                 >  6:286641/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             ccGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATc  <  2:445659/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                gTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGcc                                                                                                                              <  5:714191/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCagga                                                                                                                          <  7:233778/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCagg                                                                                                                           <  5:446906/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCagg                                                                                                                           <  5:513982/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gAGTCGATGGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCaggag                                                                                                                         >  7:834589/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                           
GCTGAAACCAGAAACGCAAACTGGCATGCACGCTCACCATAACCTGAGTCTTGGCGTGGCGAACTCTATCGCGGCGGTGGAAGAGGGCTGCGACCGAATCGACGCCAGCCTCGCGGGAATGGGCGCGGGCGCAGGTAACGCACCGCTGGAAGTGTTTATTGCCGCCGCGGATAAACTGGGCTGGCAGCATGGGACCGATCTCTATGCGTTAATGGATGCCGCCGACGACCTGGTGCGTCCGTTGCAGGATCGACCGGTACGAGTCGATCGCGAAACGCTGGCGCTGGGATACGCTGGTGTTTACTCGAGCTTCCTGCGTCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATC  >  NC_000913/374428‑374955

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: