Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A13 F1 I0 R1 37 132.8 2417326 95.5% 2308546 265.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,171,648 A→C 100% noncoding (12/120 nt) rrfB → 5S ribosomal RNA of rrnB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,6480AC99.0% 227.7 / ‑5.0 96noncoding (12/120 nt)rrfB5S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (53/41);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTA  >  NC_000913/4171480‑4171931
                                                                                                                                                                         |                                                                                                                                                                                                                                                                                             
ccGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGa                                                                                                                                                                     >  1:376404/1‑289 (MQ=14)
   ggTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGACGGTAGTGGT‑GGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                                                                  >  1:428478/1‑289 (MQ=12)
    gtgtgtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCGGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGTCAGGCATCAAATACAACGAAAg                                                                                                                                                                  >  1:380180/1‑288 (MQ=14)
       tgtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTc                                                                                                                                                              >  1:165517/1‑289 (MQ=17)
        gtAAGCGCTGCTATGCGTGGAGCTAACCGGTACCAGTGAACCGTGGGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATT‑TGAGGAGATTTTCAGCCTGATACAGAGTAAATCAGAACGTAGAAGCGGTCTGATAAAACAGAATTGGCCAGCCGAGCGTAGCGCGGTTGTCCCACCTGACCCCATGC‑CTAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTCGTGCG‑GGGTCTGCCCATGCGAGAGTAGGGAATAGTCAAGCATCAAATAAAACGTAAGGCTc                                                                                                                                                              <  2:75915/288‑1 (MQ=11)
        gtAAGCGCAGCGATGCGTTGAGCTAACGTGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGGTTTTCAGCCTGATACAGATGAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTc                                                                                                                                                              <  1:839898/288‑1 (MQ=12)
        gtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTc                                                                                                                                                              <  2:723835/288‑1 (MQ=17)
        gtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTc                                                                                                                                                              <  2:74582/288‑1 (MQ=17)
         tAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAg                                                                                                                                                            >  1:155576/1‑289 (MQ=18)
             cgcAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                                                        >  1:434995/1‑289 (MQ=21)
              gcAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaa                                                                                                                                                       >  1:376473/1‑289 (MQ=21)
              gcAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGAGACAGATTAAATCAAAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCg                                                                                                                                                         <  2:148919/287‑1 (MQ=18)
               cAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaaa                                                                                                                                                      >  2:1066979/1‑289 (MQ=25)
                 gCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAg                                                                                                                                                     >  1:849406/1‑288 (MQ=25)
                 gCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTAGAAAg                                                                                                                                                     >  1:12096/1‑288 (MQ=18)
                   gATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAc                                                                                                                                                   <  2:434995/288‑1 (MQ=25)
                            aGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCCGGCATCAAATAAAACGAAAGGCCACGCCGAAAGACTGGGCCTTTc                                                                                                                                         >  2:633561/1‑289 (MQ=255)
                               tAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGt                                                                                                                                       <  1:774545/288‑1 (MQ=25)
                                 aCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAc                                                                                                                                                   <  1:883878/274‑1 (MQ=25)
                                 aCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAc                                                                                                                                                   >  2:883878/1‑274 (MQ=25)
                                      tACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCACACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTa                                                                                                                                   >  2:191945/1‑285 (MQ=25)
                                           aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCCGAAGCGGTCTGATAAGACAGAATTTGCCTGGCGGCCGTAGCGCGGAGGTCCCACCTGACCCCATGC‑CGCACTCAGAAGTGAAACGCCGTAGCACCTCTGGTAGTGTG‑GGGTCTCACCATGCCAGAGTCGGGAGCTGTCCGGTACCAAATAAAACGAAAGGCACAGTCCAACGACCGGGCCTTGCGTATTATCttatggt                                                                                                                          >  2:1081747/1‑283 (MQ=11)
                                           aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCACATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTgttgtt                                                                                                                          >  2:876200/1‑289 (MQ=25)
                                             gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGttgtttg                                                                                                                        >  1:582231/1‑289 (MQ=25)
                                               aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGttgtttg                                                                                                                        <  2:15803/287‑1 (MQ=25)
                                                ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCg                                                                                                                     >  2:382698/1‑289 (MQ=25)
                                                ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCg                                                                                                                     >  1:829994/1‑289 (MQ=25)
                                                ccGTGAGGCTTAACCTAACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGAAACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCACCACATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGttttttgtgg                                                                                                                     >  2:933189/1‑287 (MQ=255)
                                                  gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGTTCGTCCCACCTGACCCCATGC‑CGAACTCAGAAATGCAACGCCGTAGCGCCTATGGTAGTGTG‑GGGTCTCCCCATGCTAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAATACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                   <  1:1171336/289‑1 (MQ=255)
                                                  gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                   <  1:1066979/289‑1 (MQ=25)
                                                   tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGGAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGCTGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                                                        >  2:801565/1‑251 (MQ=21)
                                                   tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGGAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                                                        <  1:801565/251‑1 (MQ=21)
                                                   tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGTGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                  >  1:325477/1‑289 (MQ=25)
                                                   tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGACAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                  >  1:370750/1‑289 (MQ=25)
                                                   tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                  >  1:543397/1‑289 (MQ=25)
                                                   tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                  >  2:153831/1‑289 (MQ=25)
                                                     aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa                                                                                                                 >  1:805379/1‑288 (MQ=25)
                                                     aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa                                                                                                                 >  1:1084011/1‑288 (MQ=25)
                                                       gCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc                                                                                                               <  2:817471/288‑1 (MQ=25)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACATATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAATTGCCAGGCATTAGATAAAACGAAAGGCCCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTTTCGGTGAAc                                                                                                               >  2:519353/1‑287 (MQ=255)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  2:582231/288‑1 (MQ=25)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  1:454583/288‑1 (MQ=25)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  1:1077256/288‑1 (MQ=25)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  1:122926/288‑1 (MQ=25)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  1:153831/288‑1 (MQ=25)
                                                        cTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCGGATACAGATGAAATCAGAACGCAGAAGCGGGCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  2:1035140/288‑1 (MQ=25)
                                                           aaCCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTGTGATAAAACAGAATTTGCCTGGCGGACGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTATTGTG‑GGGTCTCCCCATGCGAGAGTAGGGGACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTAGAAAGACTGGGCCTATAGTTTTATCTGTTGTTTGTAGGGGAACGatc                                                                                                           >  1:192978/1‑288 (MQ=11)
                                                            aCCTGACAAAGCCGAAGCTGTTTTGGCGGATA‑AGAGAGGATTTGCAGCCTGATACAGATGAAATCAGAACGCAGAAGCTGTCTGATAAAACTGAAGGTGCCTTGCGGCCGTAGCTCGGTGGGCCCACCTGTCCCCATGC‑CGAACTCAGAAGTGAGACGCCGTAGCGCCGCTGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAGCGAAAGGCTCAGTCGAAAGACTGTGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCtctc                                                                                                         <  2:840896/289‑1 (MQ=11)
                                                            aCCGTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGACGATTTTCAGCCTGATACAGATTAACTCAGAACACAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGTTAACTGCCAGGCATCAAATAAAACGAACGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCTGTGAACGCtctc                                                                                                         <  2:949102/289‑1 (MQ=255)
                                                                 aCAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCTGTCTGATAAAACAGAATTTGCCGGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAAATCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                    <  2:330563/289‑1 (MQ=25)
                                                                 aCAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                    <  2:849406/289‑1 (MQ=25)
                                                                 aCAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAAAACGCAGAAGCGGTCTGATAAAACAGACTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                    <  1:492325/289‑1 (MQ=25)
                                                                       ccGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACCCCGAAGTGAACCGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTCGACAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTTTCGTGAGTAggcc                                                                                              >  2:1073160/1‑287 (MQ=255)
                                                                       cagAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAACTCAGAACGCAGAAGCGGTCTGATAAAACAGCAGTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGTAATCAGAAGGGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGCCCTCCCCATGCGACAGTAGGGAACTGCCAGGCATCAAATAAAACGAACGGCTCAGTCGCAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGTTGACCCACCTCCTGAGCAGGa                                                                                               >  2:762952/3‑288 (MQ=11)
                                                                         gAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaa                                                                                            >  1:358986/1‑289 (MQ=25)
                                                                                ttttGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAGTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTCGGACAAATccgcc                                                                                     >  1:97292/1‑289 (MQ=18)
                                                                                ttttGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                     >  1:146913/1‑289 (MQ=18)
                                                                                 ttggtcGGATGAAGAGATGATTTTCGGCCTGAAAGTGATTAAATCAGAGCGCAGCATACGTCTGAGAAAACAGAGTGTGCCTGGCGGCAGTAGCTCGGTGGTCCCACCTGACCCCATGC‑CTAGATCATAAGTGAAGCGCCGTAACGCCGATGGTAGTGTG‑GGGTTTCCTCATGCGAGAGTAGGGAACTGCCAGGCAACAATTAAAACGAAAGGCTCAGTTGAAAGACTGGGCCTTTCTTTTGAGCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                     <  2:161983/284‑1 (MQ=1)
                                                                                      cGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAg                                                                                <  1:455420/288‑1 (MQ=18)
                                                                                          tG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAt                                                                           >  1:603715/1‑289 (MQ=17)
                                                                                                  agaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACg                                                                    <  2:358986/289‑1 (MQ=17)
                                                                                                            acTGATGCAGATTAAATCAGAACGCAGAAGCGGTCGGATAAAACGGCAGTTGCCTGGCAGCCGTAGAGCCGTGGTGCCATCTGACCCCATGC‑CGAACTCAGAGGTGAAACGCCGCAGCGCCGACGGTAGTGTG‑GAGTCTCCCCATGCGGGAGTAGGGAACTGCCAGGTATTAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGGCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                          <  2:1084011/288‑1 (MQ=11)
                                                                                                                aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTTTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCACATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTCGGTTGTCGGTGAACGCTCTCCTGGTTAGGACAAATCCGGCGGGAGCGGGTTTGAACGTTGAGAAGCAActg                                                      >  2:556940/1‑287 (MQ=11)
                                                                                                                    atatgaaaTCAGAACGCAGAAGCGGTCTGAGAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTT‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCTAAATACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCg                                                  <  2:819348/284‑1 (MQ=12)
                                                                                                                      aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGa                                                >  1:646728/1‑289 (MQ=17)
                                                                                                                       ttAAATCAGAACGCGAAAGCGGTCTGAGAAACCAGAATTGGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGGGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGTGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCg                                                  <  2:165517/286‑1 (MQ=255)
                                                                                                                          aaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTTGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAg                                                                                >  1:900140/1‑253 (MQ=17)
                                                                                                                          aaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTTGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAg                                                                                <  2:900140/253‑1 (MQ=17)
                                                                                                                          aaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                            >  1:1003444/1‑289 (MQ=17)
                                                                                                                          aaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                            >  2:983825/1‑289 (MQ=17)
                                                                                                                          aaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTCGAAAGACTGGGCCTTTTGTTTTATCTGTTGTTTGTCGGTGAACGCTCTTCTGAGTAGGACAAATCCGGCGGGAGCGGTTTTGAACGTTGCGAAGCAACGGCCCggggggc                                            >  1:916039/1‑288 (MQ=12)
                                                                                                                           aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGTACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                           >  2:310505/1‑289 (MQ=17)
                                                                                                                           aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGTACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                           <  1:310505/289‑1 (MQ=17)
                                                                                                                           aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTAAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                           >  1:913872/1‑289 (MQ=17)
                                                                                                                             cAGAACGCAGAAGCTGTCTGATAAAACAGAATTTGCCGTGCGGCCGGAGCGCGGTGGGCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCTATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc                                         <  2:646728/289‑1 (MQ=12)
                                                                                                                              aGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                        <  2:543397/289‑1 (MQ=17)
                                                                                                                                aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                                                        >  2:971605/1‑175 (MQ=17)
                                                                                                                                aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                                                        <  1:971605/175‑1 (MQ=17)
                                                                                                                                aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              <  2:314015/217‑1 (MQ=17)
                                                                                                                                aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                              >  1:314015/1‑217 (MQ=17)
                                                                                                                                aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCgg                                       >  1:1137470/1‑288 (MQ=17)
                                                                                                                                aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGACGTAGCGCCGAAGGTAGTGTG‑GGGTCTCCCCATACGAGCGTAGGGAACTGGCAGGCATCAAATAAAACGAAAGGCTCAGTCGCAAGACTGGGACTTCCGTTTTATCCGTGGTTTGTCGGTGGACGCTCTCCTTAGAAGGAGCAATCCGCCGGGAGCGGATTTGAACGTTGCGGAGCAACTGCCGGGAGGGTggggcg                                      >  2:605910/1‑285 (MQ=11)
                                                                                                                                     aGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                                 >  2:1054926/1‑289 (MQ=17)
                                                                                                                                      gAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGt                                                                                                   >  2:552744/1‑222 (MQ=17)
                                                                                                                                      gAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGt                                                                                                   <  1:552744/222‑1 (MQ=17)
                                                                                                                                            gTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccg                           >  1:914553/1‑288 (MQ=17)
                                                                                                                                             tCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACTCCATGC‑AGACCTCAGAAGTGAGACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTCGAAAGACTGGGCCTTTCGTTTTATCTTTTGTTTGTCGGTGCACTCTCTCCTGAGGAGGACACCTCCGCCGGGCGCGGATTTGAACTTTGCGAAGCAATGGCCCGGAGGGCGTCGGGCCGTAcgcccgc                          >  2:1065658/1‑288 (MQ=11)
                                                                                                                                               tGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                       >  1:315385/1‑289 (MQ=17)
                                                                                                                                                  tAAAACAGAATTGGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCGGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGTGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                    <  1:583596/289‑1 (MQ=17)
                                                                                                                                                   aaaccatgaTTGGCCGGGCGGCGCTAGCGCGGTGGCCGCACCTGA‑CCCAGGCACGAACTCAGAATGGAGACGCCGTAGCGCCGCTGTGAG‑GTGCGCGCATCCACATGCGAGAGTAGGGAATGGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCGTTCGTTTTATCTGTTGTTTGTCGGTGAACTCTCTCCTGATGAGGACAAATCCGCCGGGAGCGGATTGGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAGGCCCGCCATAAAc                   <  2:810854/281‑1 (MQ=11)
                                                                                                                                                           aTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc               <  1:382698/285‑1 (MQ=17)
                                                                                                                                                            tttGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa          <  1:983825/289‑1 (MQ=17)
                                                                                                                                                            tttGCCGGGCAGCCGGAGCGCGGTGGGTCAACCTGACCCCTGGC‑CGAACTCAGAAGTGAAAGGCCGTAACGCGGATGGGAGGTGG‑GGGTAGCCGCTGGGGAGAGTGGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCGTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCGGGGAGCGGATTTGAACGTTGCGAGGCAACGGCCCGGAGGGTGGCGGGGAGGACGCCCGCCATAAACTGCCAGGCa          <  2:570051/289‑1 (MQ=11)
                                                                                                                                                               gCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCCCTCCTGAGTAGGACAAATCCGCCCGGAGCGGCGTTGAACGTTGCGTAGCAACGGCCCGGGGGCTGGCGGGGAGGACGCCCGCAATAAAGTGACAGGccgca       >  2:1162844/1‑285 (MQ=11)
                                                                                                                                                                    gcggcCGTAGCGCGGTGGTCCCACCTGACCACATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTAGAAAGACTGGGCCTTTCGTTTGATCTGTTGTTTGTCGGTGAACGCTCTCATGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAATCAACGGCCCCGAGGGTGGCGGGCAGTACCACCGACACATACTGGCAGGTAACAAGTTa  >  2:833768/1‑289 (MQ=11)
                                                                                                                                                                    gcggcCGGAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTGTCGTTTTATCGGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt   <  1:593284/288‑1 (MQ=12)
                                                                                                                                                                         |                                                                                                                                                                                                                                                                                             
CCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATG‑AGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGC‑CGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTG‑GGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTA  >  NC_000913/4171480‑4171931

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATC  >  NC_000913/4171469‑4171945
                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         
ACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAA                                                                                                                                                                                  >  M01186:134:000000000‑AAJAW:1:1112:10751:11872/1‑301 (MQ=0)
   TTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGACGGTAGTGGTGGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG                                                                                                                                                                               >  M01186:134:000000000‑AAJAW:1:1114:16355:2725/1‑301 (MQ=4)
    TGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCGGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGTCAGGCATCAAATACAACGAAAGG                                                                                                                                                                               >  M01186:134:000000000‑AAJAW:1:1112:10381:14040/1‑300 (MQ=0)
     GATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAAAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGC                                                                                                                                                                              <  M01186:134:000000000‑AAJAW:1:1110:17126:21884/300‑1 (MQ=6)
       TAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCA                                                                                                                                                                           >  M01186:134:000000000‑AAJAW:1:1106:11676:3920/1‑301 (MQ=16)
         GGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGT                                                                                                                                                                         >  M01186:134:000000000‑AAJAW:1:1105:6694:18125/1‑301 (MQ=20)
             GGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAGATAAAACGAAAGGCTCAGTCGAA                                                                                                                                                                     >  M01186:134:000000000‑AAJAW:1:1114:23845:7716/1‑301 (MQ=23)
              GGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAA                                                                                                                                                                    >  M01186:134:000000000‑AAJAW:1:1112:21410:11910/1‑301 (MQ=28)
                 GTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGA                                                                                                                                                                  >  M01186:134:000000000‑AAJAW:1:2108:22825:13492/1‑300 (MQ=31)
                 GTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTAGAAAGA                                                                                                                                                                  >  M01186:134:000000000‑AAJAW:1:1101:12864:10068/1‑300 (MQ=14)
                 GTGTAAGCGCAGCGATGCGTTGAGATAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCAAGCCGAACTCAGAAGTGAAAAGCCTTAGCGCCGCTGGGAGTGTGGGGTCTCCACATGCGATAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGATGCGTCGAAAGA                                                                                                                                                                  >  M01186:134:000000000‑AAJAW:1:1116:18161:25130/1‑300 (MQ=4)
                  TGTAAGCGCAGCGATGCGTTGAGCTAACGTGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGGTTTTCAGCCTGATACAGATGAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC                                                                                                                                                                 <  M01186:134:000000000‑AAJAW:1:2108:5342:7488/300‑1 (MQ=31)
                                 GCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCG                                                                                                                                                      <  M01186:134:000000000‑AAJAW:1:2109:12633:15732/296‑1 (MQ=36)
                                         CTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTG                                                                                                                                          <  M01186:134:000000000‑AAJAW:1:2106:10615:4135/300‑1 (MQ=36)
                                             CCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGT                                                                                                                                     >  M01186:134:000000000‑AAJAW:1:1118:17463:19334/1‑301 (MQ=36)
                                                GTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGG                                                                                                                                  >  M01186:134:000000000‑AAJAW:1:2107:17336:21221/1‑301 (MQ=36)
                                                   CTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGGAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCC                                                                                                                                                           <  M01186:134:000000000‑AAJAW:1:2106:22398:23733/273‑1 (MQ=37)
                                                   CTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGTGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGA                                                                                                                               >  M01186:134:000000000‑AAJAW:1:1110:7318:21849/1‑301 (MQ=36)
                                                   CTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGACAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGA                                                                                                                               >  M01186:134:000000000‑AAJAW:1:1112:13641:8546/1‑301 (MQ=36)
                                                   CTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGA                                                                                                                               >  M01186:134:000000000‑AAJAW:1:1117:22630:15064/1‑301 (MQ=36)
                                                     AATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAA                                                                                                                              >  M01186:134:000000000‑AAJAW:1:2107:18998:4613/1‑300 (MQ=36)
                                                     AATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAA                                                                                                                              >  M01186:134:000000000‑AAJAW:1:2116:16831:5166/1‑300 (MQ=36)
                                                           CCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTGTGATAAAACAGAATTTGCCTGGCGGACGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTATTGTGGGGTCTCCCCATGCGAGAGTAGGGGACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTAGAAAGACTGGGCCTATAGTTTTATCTGTTGTTTGTAGGGGAACGATCT                                                                                                                        >  M01186:134:000000000‑AAJAW:1:1106:16285:21422/1‑300 (MQ=32)
                                                            CGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGTTCGTCCCACCTGACCCCATGCCGAACTCAGAAATGCAACGCCGTAGCGCCTATGGTAGTGTGGGGTCTCCCCATGCTAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAATACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCC                                                                                                                      <  M01186:134:000000000‑AAJAW:1:2118:5350:20169/301‑1 (MQ=33)
                                                            CGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCC                                                                                                                      <  M01186:134:000000000‑AAJAW:1:2115:8271:14420/301‑1 (MQ=36)
                                                                  GCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGT                                                                                                                 <  M01186:134:000000000‑AAJAW:1:1114:23210:19822/300‑1 (MQ=36)
                                                                  GCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGT                                                                                                                 <  M01186:134:000000000‑AAJAW:1:1104:20156:18158/300‑1 (MQ=36)
                                                                  GCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGT                                                                                                                 <  M01186:134:000000000‑AAJAW:1:1105:4920:16958/300‑1 (MQ=36)
                                                                  GCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGT                                                                                                                 <  M01186:134:000000000‑AAJAW:1:2115:21420:21519/300‑1 (MQ=36)
                                                                  GCATAACCTTACAACGCAGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAAATTAAATCAGAACGCAGAAGCCGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGT                                                                                                                 <  M01186:134:000000000‑AAJAW:1:1106:25466:8839/300‑1 (MQ=34)
                                                                         CTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTTCCCTGAGTAGGACAAA                                                                                                         >  M01186:134:000000000‑AAJAW:1:1111:17476:22283/1‑301 (MQ=35)
                                                                           TACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAAAACGCAGAAGCGGTCTGATAAAACAGACTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATC                                                                                                       <  M01186:134:000000000‑AAJAW:1:1116:9122:2678/301‑1 (MQ=35)
                                                                                CGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAGTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTCGGACAAATCCGCCG                                                                                                  >  M01186:134:000000000‑AAJAW:1:1104:17043:1301/1‑301 (MQ=35)
                                                                                CGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCG                                                                                                  >  M01186:134:000000000‑AAJAW:1:1105:14910:12891/1‑301 (MQ=36)
                                                                                          GTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATT                                                                                        >  M01186:134:000000000‑AAJAW:1:1119:8232:12669/1‑301 (MQ=28)
                                                                                                GCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC                                                                                   <  M01186:134:000000000‑AAJAW:1:1114:9845:20470/300‑1 (MQ=25)
                                                                                                                     GCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAG                                                             >  M01186:134:000000000‑AAJAW:1:2101:25236:20841/1‑301 (MQ=16)
                                                                                                                         GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTTGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC                                                                                   >  M01186:134:000000000‑AAJAW:1:2110:14903:6597/1‑275 (MQ=17)
                                                                                                                         GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTG                                                         >  M01186:134:000000000‑AAJAW:1:2113:9254:13309/1‑301 (MQ=16)
                                                                                                                         GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTCGAAAGACTGGGCCTTTTGTTTTATCTGTTGTTTGTCGGTGAACGCTCTTCTGAGTAGGACAAATCCGGCGGGAGCGGTTTTGAACGTTGCGAAGCAACGGCCCGGGGGGCG                                                         >  M01186:134:000000000‑AAJAW:1:2110:9492:16805/1‑301 (MQ=15)
                                                                                                                          ATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTAAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGG                                                        >  M01186:134:000000000‑AAJAW:1:2110:18468:15346/1‑301 (MQ=16)
                                                                                                                               GATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGG                                                    >  M01186:134:000000000‑AAJAW:1:2117:21531:18765/1‑300 (MQ=16)
                                                                                                                               GATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGT                                                                                                                 >  M01186:134:000000000‑AAJAW:1:1110:10573:14082/1‑239 (MQ=18)
                                                                                                                               GATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCC                                                                                                                                                           <  M01186:134:000000000‑AAJAW:1:2112:26889:12810/197‑1 (MQ=19)
                                                                                                                                    AATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGTACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAC                                              <  M01186:134:000000000‑AAJAW:1:1110:16253:12037/301‑1 (MQ=16)
                                                                                                                                     ATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCC                                                                                                      <  M01186:134:000000000‑AAJAW:1:1117:23533:21433/244‑1 (MQ=18)
                                                                                                                                       CAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACC                                                                                                                                                                                                                                                                                        <  M01186:134:000000000‑AAJAW:1:1118:8291:8354/64‑1 (MQ=0)
                                                                                                                                           ACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGC                                        >  M01186:134:000000000‑AAJAW:1:2110:16488:15792/1‑300 (MQ=16)
                                                                                                                                              CAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATA                                    >  M01186:134:000000000‑AAJAW:1:1110:7021:14879/1‑301 (MQ=16)
                                                                                                                                                 AAGCGGTTTAATAAAACACAATTTGCCTGGCGGCCGTGGCGCGGTGGTCCCACCTGACCCCAGGCCGAACTCAGAAGTGAACCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGCCTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAC                                 <  M01186:134:000000000‑AAJAW:1:2102:19831:13079/301‑1 (MQ=15)
                                                                                                                                                          GATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAG                                                                                                                                                                                                                                        >  M01186:134:000000000‑AAJAW:1:2110:15765:15442/1‑93 (MQ=0)
                                                                                                                                                           ATAAAACAGAATTGGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCGGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGTGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCAT                       <  M01186:134:000000000‑AAJAW:1:1118:24067:20337/301‑1 (MQ=16)
                                                                                                                                                                    AATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAT                  <  M01186:134:000000000‑AAJAW:1:1112:21159:15585/297‑1 (MQ=16)
                                                                                                                                                                     TTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGC             <  M01186:134:000000000‑AAJAW:1:2112:17844:20995/301‑1 (MQ=16)
                                                                                                                                                                             GGCGGCCGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCAGGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTGTCGTTTTATCGGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGC      <  M01186:134:000000000‑AAJAW:1:1119:23217:6090/300‑1 (MQ=16)
                                                                                                                                                                                 GCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATC  >  M01186:134:000000000‑AAJAW:1:1103:15862:13339/1‑300 (MQ=17)
                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         
ACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATC  >  NC_000913/4171469‑4171945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑