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 152 24.9 2083062 91.6% 1908084 109.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,354,645 +C intergenic (+98/‑15) h16_A1248 → / → h16_A1249 predicted signal transduction protein containingCHASE, PAS, GGDEF and EAL domains/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,354,6431.C100.0% 149.4 / NA 39intergenic (+96/‑17)h16_A1248/h16_A1249predicted signal transduction protein containingCHASE, PAS, GGDEF and EAL domains/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (23/16);  total (23/16)

CTCTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCG  >  AM260479/1354513‑1354767
                                                                                                                                   |                                                                                                                            
ctctCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGccc                                                                                                                           >  1:267335/1‑133 (MQ=255)
                 cAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGccc                                                                                                                           >  2:15337/1‑116 (MQ=255)
                 cAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGccc                                                                                                                           <  1:15337/117‑2 (MQ=255)
                     cGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAgg                                                                                                       <  1:1030430/134‑1 (MQ=255)
                                  aCTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCtgctgc                                                                                           <  1:9058/133‑1 (MQ=255)
                                      cccGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACttt                                                                                                                    >  2:553389/1‑104 (MQ=255)
                                      cccGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACttt                                                                                                                    <  1:553389/104‑1 (MQ=255)
                                                           tAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCg                                                                 >  1:1856/1‑134 (MQ=255)
                                                            aaGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  <  2:367613/100‑1 (MQ=255)
                                                            aaGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  >  1:367613/1‑100 (MQ=255)
                                                                     aaTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTg                                                        >  2:534085/1‑133 (MQ=255)
                                                                      aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  >  1:12328/1‑90 (MQ=255)
                                                                      aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  <  1:194297/90‑1 (MQ=255)
                                                                      aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  >  1:110874/1‑90 (MQ=255)
                                                                      aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  >  2:194297/1‑90 (MQ=255)
                                                                      aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  <  2:12328/90‑1 (MQ=255)
                                                                      aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  <  2:110874/90‑1 (MQ=255)
                                                                           aGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGACCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAgg                                                 >  1:325649/1‑134 (MQ=255)
                                                                                gCTAGATCATGCAGGACCGTACTATTGAAATTATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAg                                                                                      <  2:944836/92‑1 (MQ=255)
                                                                                gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  >  1:939869/1‑80 (MQ=255)
                                                                                gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                  <  2:939869/80‑1 (MQ=255)
                                                                                gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAg                                                                                      >  1:944836/1‑92 (MQ=255)
                                                                                gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCg                                               >  2:723201/1‑131 (MQ=255)
                                                                                    gATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGc                                         <  1:113777/133‑1 (MQ=255)
                                                                                            aGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGATTCGCTGGAGCAGGCGGAGAGCGGCTACCgg                                >  2:845631/1‑134 (MQ=255)
                                                                                             ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCgg                                >  1:538020/1‑133 (MQ=255)
                                                                                               aCCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCaa                             >  1:988996/1‑134 (MQ=255)
                                                                                                  gTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGc                                                                   >  1:843580/1‑93 (MQ=255)
                                                                                                  gTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGc                                                                   <  2:843580/93‑1 (MQ=255)
                                                                                                   tACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTcc                         >  1:496804/1‑134 (MQ=255)
                                                                                                   tACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTcc                         >  2:885618/1‑134 (MQ=255)
                                                                                                       aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                       >  2:41227/1‑132 (MQ=255)
                                                                                                            aaTCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTgg                >  2:210117/1‑134 (MQ=255)
                                                                                                              tCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcg              <  1:210117/134‑1 (MQ=255)
                                                                                                                 tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTAcc                                  >  1:1014937/1‑111 (MQ=255)
                                                                                                                 tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTAcc                                  <  2:1014937/111‑1 (MQ=255)
                                                                                                                      cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTAc                                   <  2:380550/105‑1 (MQ=255)
                                                                                                                      cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTAc                                   >  1:380550/1‑105 (MQ=255)
                                                                                                                      cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                       <  2:833813/117‑1 (MQ=255)
                                                                                                                      cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                       >  1:833813/1‑117 (MQ=255)
                                                                                                                        ggCCATCGACGCCCTGGACTTTTCCAGGATGAAGACCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATc    <  2:325649/134‑1 (MQ=255)
                                                                                                                          ccATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCg  >  1:332933/1‑134 (MQ=255)
                                                                                                                                   |                                                                                                                            
CTCTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCG  >  AM260479/1354513‑1354767

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGAGAACGGTCCCGCTCTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCCCGTCGTCCC  >  AM260479/1354498‑1354784
                                                                                                                                                    |                                                                                                                                           
TGGAGAACGGTCCCGCTCTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCct                                                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1201:19943:66848/1‑148 (MQ=60)
        cgtgggctcggagatgtgtataagagacaGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACT                                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1208:12288:68465/118‑1 (MQ=60)
        cgtgggctcggagatgtgtataagagacaGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACT                                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1202:11596:75357/118‑1 (MQ=60)
              cagCTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGA                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1103:18403:12939/147‑1 (MQ=60)
                       gtgtataagagacaGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCC                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2102:11918:17121/135‑1 (MQ=60)
                                   CCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATC                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2216:8601:41459/150‑1 (MQ=60)
                                        cgtgggctcggagatgtgtataagagacaGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAG                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1112:2546:41545/121‑1 (MQ=60)
                                                GACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGC                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1101:15737:67252/149‑1 (MQ=60)
                                                     CCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGTTCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCctggacctgtctcttatacacatctccgagcccacgagacggagctacatctcgt                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:2116:12004:70833/1‑95 (MQ=60)
                                                           AGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGG                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1101:11164:14835/1‑150 (MQ=60)
                                                            GTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCtgtctcttatacacatctc                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1208:8224:33961/1‑131 (MQ=60)
                                                                      TGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCtgtctcttatacacatctccgagcccacg                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1101:10955:90630/1‑121 (MQ=60)
                                                                      TGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCtgtctcttatacacatctccgagcccacg                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1108:10879:48785/1‑121 (MQ=60)
                                                                           AAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGACCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGC                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1204:4971:67741/1‑150 (MQ=60)
                                                                                GTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCtgtctcttatacacatctccgagcccacgagacggagct                                                            >  GWNJ‑0478:712:GW2002102894th:2:2210:15214:45361/1‑111 (MQ=60)
                                                                                GTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGctgtctcttatacacatctccgagccc                                                             >  GWNJ‑0478:712:GW2002102894th:2:2210:13008:73841/1‑122 (MQ=60)
                                                                                             GAGCTAGGTCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGC                                                >  GWNJ‑0478:712:GW2002102894th:2:2103:6007:7286/1‑149 (MQ=60)
                                                                                               GCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAAT                                             >  GWNJ‑0478:712:GW2002102894th:2:2214:7107:21714/1‑150 (MQ=60)
                                                                                                  AGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCctgtctcttatacacatctccgagccc                                          >  GWNJ‑0478:712:GW2002102894th:2:2205:6964:53715/1‑123 (MQ=60)
                                                                                                  AGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTC                                           <  GWNJ‑0478:712:GW2002102894th:2:1108:3783:66808/149‑1 (MQ=60)
                                                                                                   GATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCT                                         >  GWNJ‑0478:712:GW2002102894th:2:1214:2456:52289/1‑150 (MQ=60)
                                                                                                                 GTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCctgtctctt                           >  GWNJ‑0478:712:GW2002102894th:2:2215:4213:57981/1‑141 (MQ=60)
                                                                                                                      ATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGtctcttatacac                      >  GWNJ‑0478:712:GW2002102894th:2:1209:18593:11856/1‑138 (MQ=60)
                                                                                                                      ATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCtg                      >  GWNJ‑0478:712:GW2002102894th:2:2204:10682:87928/1‑148 (MQ=60)
                                                                                                                          AAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGG                  >  GWNJ‑0478:712:GW2002102894th:2:1207:20755:18023/1‑150 (MQ=60)
                                                                                                                            ATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGAC                <  GWNJ‑0478:712:GW2002102894th:2:1115:11958:28757/150‑1 (MQ=60)
                                                                                                                                          CATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCCCGTCGTCCC  >  GWNJ‑0478:712:GW2002102894th:2:1207:6858:28000/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                           
TGGAGAACGGTCCCGCTCTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCCCGTCGTCCC  >  AM260479/1354498‑1354784

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