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
A1 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 547,879 2 bp→GA coding (595‑596/630 nt) h16_A0522 ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479547,8790AG100.0% 135.4 / NA 35L199S (CTC→TCC) h16_A0522conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/18);  total (17/18)
*AM260479547,8800GA100.0% 133.2 / NA 35L199S (CTC→TCC) h16_A0522conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

GCAGCTCTTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCAACAGGTTGTAGCGCGCGCCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGC  >  AM260479/547754‑547988
                                                                                                                              ||                                                                                                            
gCAGCTCTTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGgacagcag                                                                                                        <  2:866757/133‑1 (MQ=255)
    ctctTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCca                                                                                                     >  2:512540/1‑132 (MQ=255)
     tctTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAg                                                                                                  >  2:104922/1‑134 (MQ=255)
        tCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGcc                                                                                               >  2:998860/1‑134 (MQ=255)
                  gTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGc                                                                                     >  2:457668/1‑134 (MQ=255)
                       aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGaa                                                                                >  2:185542/1‑134 (MQ=255)
                       aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGaa                                                                                >  1:715792/1‑134 (MQ=255)
                          agcagcGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAAtt                                                                              <  1:512540/133‑1 (MQ=255)
                           gcagcGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            <  2:715792/134‑1 (MQ=255)
                            cagcGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCa                                                                           <  2:207912/134‑1 (MQ=255)
                                 aGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCgtg                                                                      <  1:457668/134‑1 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGgacagcag                                                                                                        >  2:664623/1‑80 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGgacagcag                                                                                                        <  1:664623/80‑1 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTg                                                                                      >  1:907449/1‑98 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTg                                                                                      <  2:907449/98‑1 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            >  1:333011/1‑108 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            <  2:333011/108‑1 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTa                                                                 >  1:458405/1‑119 (MQ=255)
                                                     aTGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTa                                                                 <  2:458405/119‑1 (MQ=255)
                                                          tCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTa                                                                 >  2:433637/1‑114 (MQ=255)
                                                          tCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTa                                                                 <  1:433637/114‑1 (MQ=255)
                                                                      ggTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCcg                                 >  1:437525/1‑134 (MQ=255)
                                                                       gTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTgcg                                                      >  1:451639/1‑112 (MQ=255)
                                                                       gTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTgcg                                                      <  2:451639/112‑1 (MQ=255)
                                                                       gTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCcgc                                <  1:998860/134‑1 (MQ=255)
                                                                        ttGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCTCAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCcgcg                               <  2:437525/134‑1 (MQ=255)
                                                                             gcgcgcgcCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGt                           <  1:641206/133‑1 (MQ=255)
                                                                             gcgcgcgcCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGt                           <  1:185542/133‑1 (MQ=255)
                                                                                       gCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCgcgca                             >  2:580992/1‑121 (MQ=255)
                                                                                       gCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCgcgca                             <  1:580992/122‑2 (MQ=255)
                                                                                       gCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTgccgcc                 >  2:152421/1‑134 (MQ=255)
                                                                                       gCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTgccgcc                 >  2:58824/1‑134 (MQ=255)
                                                                                              agcagcGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGctg          <  1:58824/134‑1 (MQ=255)
                                                                                                      ggcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGc  <  2:1011942/134‑1 (MQ=255)
                                                                                                         gAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            <  2:208494/57‑1 (MQ=255)
                                                                                                         gAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            <  2:672098/57‑1 (MQ=255)
                                                                                                         gAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            >  1:672098/1‑57 (MQ=255)
                                                                                                         gAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                            >  1:208494/1‑57 (MQ=255)
                                                                                                                              ||                                                                                                            
GCAGCTCTTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCAACAGGTTGTAGCGCGCGCCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGC  >  AM260479/547754‑547988

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGCAGCTCTTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCAACAGGTTGTAGCGCGCGCCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAGCGTGCGTCGATGCGGCGCAGGCGCTTG  >  AM260479/547752‑548022
                                                                                                                                |                                                                                                                                               
agctccgtctcgtgggctcggagatgtgtataagagacaGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTG                                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2111:7578:73080/111‑1 (MQ=60)
          TCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAAT                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2112:10458:75005/1‑150 (MQ=60)
                           CAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGG                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1216:4387:62412/149‑1 (MQ=60)
                                  GAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGC                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1212:20176:69178/150‑1 (MQ=60)
                                       agacagGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGC                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1213:15928:34522/144‑1 (MQ=60)
                                        GTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTctgtctcttatacacatctccg                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2208:11583:80894/1‑128 (MQ=60)
                                        GTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCtgtctcttat                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1207:17500:33652/1‑139 (MQ=60)
                                        GTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCC                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:1212:19584:73170/1‑150 (MQ=60)
                                                         GAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGC                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1213:15154:55848/1‑150 (MQ=60)
                                                          AATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCctgtctct                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1212:11627:36000/1‑142 (MQ=60)
                                                              GTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCtgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaa                                                              >  GWNJ‑0478:712:GW2002102894th:2:1211:8273:68818/1‑82 (MQ=60)
                                                                   CAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCtgtctcttatacacatctccgagcccacgagacggagctacatctcgtatg                                                        >  GWNJ‑0478:712:GW2002102894th:2:1110:12966:43276/1‑99 (MQ=60)
                                                                        GGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCaatgcgcctacac                                                   <  GWNJ‑0478:712:GW2002102894th:2:2215:2020:5382/150‑14 (MQ=60)
                                                                            GTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCaatgcgcctacaccac                                                <  GWNJ‑0478:712:GW2002102894th:2:2104:19774:42970/149‑17 (MQ=60)
                                                                              AGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGC                                              <  GWNJ‑0478:712:GW2002102894th:2:2109:3286:34536/149‑1 (MQ=60)
                                                                              AGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGC                                              <  GWNJ‑0478:712:GW2002102894th:2:1113:6712:94343/149‑1 (MQ=60)
                                                                                       CCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaataaaagttatcagatgttt                                     >  GWNJ‑0478:712:GW2002102894th:2:2111:10298:46279/1‑53 (MQ=60)
                                                                                            TGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCtgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgc                                >  GWNJ‑0478:712:GW2002102894th:2:2110:9914:17302/1‑88 (MQ=60)
                                                                                            TGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCtgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgc                                >  GWNJ‑0478:712:GW2002102894th:2:1115:8000:27064/1‑88 (MQ=60)
                                                                                               cgtctcgtgggctcggagatgtgtataagagacAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAG                             <  GWNJ‑0478:712:GW2002102894th:2:1202:13537:83937/117‑1 (MQ=60)
                                                                                               CAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAG                             <  GWNJ‑0478:712:GW2002102894th:2:1104:7488:7253/150‑1 (MQ=60)
                                                                                                                          ACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAGCtgtctcttatacacatctccgagccc  >  GWNJ‑0478:712:GW2002102894th:2:2212:8530:36623/1‑124 (MQ=60)
                                                                                                                                |                                                                                                                                               
TGGCAGCTCTTCGGTGAACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCAACAGGTTGTAGCGCGCGCCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAGCGTGCGTCGATGCGGCGCAGGCGCTTG  >  AM260479/547752‑548022

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: