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 148 28.0 2229078 88.9% 1981650 107.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,998,635 +C intergenic (‑221/+527) h16_A3705 ← / ← h16_A3706 transcriptional regulator, LysR‑family/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,998,6331.C100.0% 145.6 / NA 38intergenic (‑219/+529)h16_A3705/h16_A3706transcriptional regulator, LysR‑family/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (23/15);  total (23/15)

ATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCC  >  AM260479/3998510‑3998761
                                                                                                                            |                                                                                                                                
atatGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGc                                                                                                                         >  2:793861/1‑134 (MQ=255)
        gAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCa                                                                                                                 <  1:859536/134‑1 (MQ=255)
                      cTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCgccgccgcc                                                                                                   >  2:598763/1‑134 (MQ=255)
                            aTCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCt                                                                                                 <  2:80316/130‑1 (MQ=255)
                            aTCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCt                                                                                                 >  1:80316/1‑130 (MQ=255)
                            aTCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCg                                                                                             >  1:646321/1‑134 (MQ=255)
                            aTCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCg                                                                                             >  1:777894/1‑134 (MQ=255)
                                 aTGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCg                                                                                         >  2:340971/1‑133 (MQ=255)
                                 aTGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGACACCCTCGCCGCCGCCCTGTCGGCCg                                                                                         >  2:150943/1‑133 (MQ=255)
                                  tGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCa                                                                                       >  1:984002/1‑134 (MQ=255)
                                       cAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGcc                                                                                                                                >  1:267859/1‑87 (MQ=255)
                                       cAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGcc                                                                                                                                <  2:267859/87‑1 (MQ=255)
                                                   gTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGa                                                                       <  2:984002/133‑1 (MQ=255)
                                                   gTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGACACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGa                                                                       <  1:150943/133‑1 (MQ=255)
                                                        tcgtcgAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGc                                                                                           <  1:953598/108‑1 (MQ=255)
                                                        tcgtcgAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGc                                                                                           >  2:953598/1‑108 (MQ=255)
                                                          gtcgAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAg                                                               >  2:918770/1‑134 (MQ=255)
                                                          gtcgAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAg                                                               >  1:1096410/1‑134 (MQ=255)
                                                             gAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTc                                                            >  2:641754/1‑134 (MQ=255)
                                                             gAAGTGACGGCTCAGGCATCCATGCAACAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTATCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGtt                                                             <  1:422049/133‑1 (MQ=255)
                                                                  gACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACCGCAAGTTCGAtgt                                                       >  2:803253/1‑134 (MQ=255)
                                                                             cATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACg                                            <  1:641754/134‑1 (MQ=255)
                                                                               tCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGgc                                          <  2:646321/134‑1 (MQ=255)
                                                                                            gaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTc                                                                                                                     >  1:995192/1‑46 (MQ=38)
                                                                                            gaAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTc                                                                                                                     <  2:995192/46‑1 (MQ=38)
                                                                                              aaGCCTGAATCATGAAAACCAACCGCTCCGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAg                           >  2:927163/1‑134 (MQ=255)
                                                                                              aaGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAg                           >  1:840249/1‑134 (MQ=255)
                                                                                               aGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGt                          >  2:253812/1‑134 (MQ=255)
                                                                                                gCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGtt                         <  1:918770/134‑1 (MQ=255)
                                                                                                 ccTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGtt                                                             >  1:998451/1‑97 (MQ=255)
                                                                                                 ccTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGtt                                                             <  2:998451/97‑1 (MQ=255)
                                                                                                 ccTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTc                        >  2:869724/1‑134 (MQ=255)
                                                                                                   tGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCg                                                                                    <  2:440520/72‑1 (MQ=255)
                                                                                                   tGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCg                                                                                    >  1:440520/1‑72 (MQ=255)
                                                                                                         aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTc                        <  2:775732/126‑1 (MQ=255)
                                                                                                         aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTc                        >  1:775732/1‑126 (MQ=255)
                                                                                                         aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCg                  >  2:668624/1‑132 (MQ=255)
                                                                                                         aTGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCg                  <  1:668624/132‑1 (MQ=255)
                                                                                                                     cGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGgcgc    <  2:777894/134‑1 (MQ=255)
                                                                                                                       cTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACGCCGACGGCGccc  >  1:832526/1‑134 (MQ=255)
                                                                                                                            |                                                                                                                                
ATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCC  >  AM260479/3998510‑3998761

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
ATTGTGGATCCTGATATTTATTGATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCCGCACCGGCAAGTTCGATCCGTAC  >  AM260479/3998487‑3998784
                                                                                                                                                     |                                                                                                                                                     
atgtgtataagagacagTTATTGATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCC                                                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2215:6608:76405/133‑1 (MQ=60)
ATTGTGGATCCTGATATTTATTGATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCC                                                                                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:2215:14926:33603/1‑150 (MQ=60)
 gggctcggagatgtgtataagagacaGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCC                                                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2215:5866:64651/123‑1 (MQ=60)
        TCCTGATATTTATTGATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCT                                                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1114:13875:76168/1‑150 (MQ=60)
                              GGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCC                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2203:19356:90040/150‑1 (MQ=60)
                               nnnnnngnnnncggagatgtgtataagagacagATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCC                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2212:9050:101304/116‑1 (MQ=60)
                                    GCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGG                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:1107:17216:18139/1‑150 (MQ=60)
                                    GCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGG                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2107:1536:90285/1‑150 (MQ=60)
                                    GCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGG                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2115:6359:29309/1‑150 (MQ=60)
                                          TATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCAC                                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2211:1792:30133/1‑150 (MQ=60)
                                               AGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCctgtctcttatacacatctccgagcccacg                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1116:1808:69727/1‑117 (MQ=60)
                                                     ataagagacagATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGC                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2207:13097:57441/138‑1 (MQ=60)
                                                                  GCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGT                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2216:10241:12075/1‑150 (MQ=60)
                                                                         AGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGACACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATG                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1111:7354:35666/149‑1 (MQ=60)
                                                                                   CGAAGTGACGGCTCAGGCATCCATGCAACAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTATCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGA                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1211:3938:15377/149‑1 (MQ=60)
                                                                                                   GCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCA                                                    <  GWNJ‑0478:712:GW2002102894th:2:2107:13147:63612/150‑1 (MQ=60)
                                                                                                    CATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaa                                                    >  GWNJ‑0478:712:GW2002102894th:2:2211:3182:90689/1‑76 (MQ=60)
                                                                                                      TCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGT                                                 >  GWNJ‑0478:712:GW2002102894th:2:2201:10062:66359/1‑150 (MQ=60)
                                                                                                         ATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACtgtctcttatacacatctccg                                               >  GWNJ‑0478:712:GW2002102894th:2:2210:17048:8954/1‑128 (MQ=60)
                                                                                                           GCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACctgtctcttatacacatctccgagcccacgagacactcgctaatctcg                                            >  GWNJ‑0478:712:GW2002102894th:2:1210:7192:15576/1‑102 (MQ=60)
                                                                                                                 GAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTAC                                      >  GWNJ‑0478:712:GW2002102894th:2:2115:11255:18650/1‑150 (MQ=60)
                                                                                                                      AGCCTGAATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGA                                 <  GWNJ‑0478:712:GW2002102894th:2:2209:16646:53249/150‑1 (MQ=60)
                                                                                                                             ATCATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCC                          <  GWNJ‑0478:712:GW2002102894th:2:2108:10494:28935/150‑1 (MQ=60)
                                                                                                                               CATGAAAACCAACCGCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACGCCGACGGCGCCCG                        >  GWNJ‑0478:712:GW2002102894th:2:2201:5016:15497/1‑150 (MQ=60)
                                                                                                                                             GCTACGCCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCtgtctcttatacacatctc          >  GWNJ‑0478:712:GW2002102894th:2:1207:8918:88441/1‑131 (MQ=60)
                                                                                                                                                     CTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCCGCACCGGCAAGTTCGATCCGTAC  >  GWNJ‑0478:712:GW2002102894th:2:1206:1968:46439/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
ATTGTGGATCCTGATATTTATTGATATGAGGGAAATGCCTAATATCTAGCCATCGGATGCCGCAATGCAACATAGTCTTTCGTCGAAGTGACGGCTCAGGCATCCATGCAATAGAGAAAGCCTGAATCATGAAAACCAACCGCTACGCCTGATCGCTGTCGCCACCCTCGCCGCCGCCCTGTCGGCCGCACCGGCCTTTGCCAAGAACGGCAAGTTCGATGTGTTCACCGACGGCGCCAAGGTAGCCAAGTTCGATCCGTACACCGACGGCGCCCGCACCGGCAAGTTCGATCCGTAC  >  AM260479/3998487‑3998784

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: