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
A0 F0 I1 R1 155 25.8 2285122 88.0% 2010907 104.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 2,843,580 +G intergenic (‑1/+68) pip ← / ← h16_B2516 prolyl aminopeptidase/transcriptional regulator, LysR‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604802,843,5791.G97.9% 178.7 / ‑3.8 48coding (1/957 nt)pipprolyl aminopeptidase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (25/22);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGC  >  AM260480/2843463‑2843702
                                                                                                                      |                                                                                                                           
cgGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACg                                                                                                               >  1:862465/1‑132 (MQ=255)
        cgcCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                     <  2:862465/134‑1 (MQ=255)
        cgcCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                     <  2:445169/134‑1 (MQ=255)
           cGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGc                                                                                                   >  1:1123898/1‑133 (MQ=255)
                                agtagATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATggg                                                                                                                     >  2:608112/1‑94 (MQ=255)
                                agtagATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATggg                                                                                                                     <  1:608112/94‑1 (MQ=255)
                                agtagATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATggg                                                                                                                     <  1:485282/94‑1 (MQ=255)
                                agtagATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATggg                                                                                                                     >  2:485282/1‑94 (MQ=255)
                                        aCATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCaaa                                                                       >  2:787707/1‑132 (MQ=255)
                                          aTGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGAt                                                                                                          <  1:804882/95‑1 (MQ=255)
                                          aTGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGAt                                                                                                          >  2:804882/1‑95 (MQ=255)
                                          aTGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                   <  2:236878/134‑1 (MQ=255)
                                                gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATTGATATGGGTCTACGGAGATCGGGAGc                                                                                                   <  2:588959/96‑1 (MQ=255)
                                                gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCg                                                                                                        <  2:267836/91‑1 (MQ=255)
                                                gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCg                                                                                                        >  1:267836/1‑91 (MQ=255)
                                                gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGc                                                                                                   >  1:588959/1‑96 (MQ=255)
                                                gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGc                                                                                       <  1:649153/108‑1 (MQ=255)
                                                gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGc                                                                                       >  2:649153/1‑108 (MQ=255)
                                                       aCATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCg                                                                                     <  1:308331/103‑1 (MQ=255)
                                                       aCATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCg                                                                                     >  2:308331/1‑103 (MQ=255)
                                                                  gCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCATAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCg                                               >  1:862816/1‑130 (MQ=255)
                                                                  gCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCATAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCg                                               <  2:862816/130‑1 (MQ=255)
                                                                     ggTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGgaga                                                                                                           >  1:591915/1‑67 (MQ=255)
                                                                     ggTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGgaga                                                                                                           <  2:591915/67‑1 (MQ=255)
                                                                              gggCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCAcc                                <  2:1123898/133‑1 (MQ=255)
                                                                                    gATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                     <  1:246319/58‑1 (MQ=255)
                                                                                    gATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                     >  2:246319/1‑58 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGgaga                                                                                                           >  1:797366/1‑51 (MQ=39)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGgaga                                                                                                           <  2:797366/51‑1 (MQ=39)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                     <  2:395659/57‑1 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                     >  1:395659/1‑57 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGc                                                                                                   <  1:263462/59‑1 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGc                                                                                                   >  2:263462/1‑59 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGaa                                                                              <  2:91830/80‑1 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGaa                                                                              >  1:91830/1‑80 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                   <  1:1019147/91‑1 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                   >  2:1019147/1‑91 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGa                                              >  1:424655/1‑112 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGa                                              <  2:424655/112‑1 (MQ=255)
                                                                                     aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCt                        >  1:424822/1‑134 (MQ=255)
                                                                                                       cgcAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                   <  2:1089145/73‑1 (MQ=255)
                                                                                                       cgcAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                   >  1:1089145/1‑73 (MQ=255)
                                                                                                         cAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGa                                              <  2:499797/92‑1 (MQ=255)
                                                                                                         cAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGa                                              >  1:499797/1‑92 (MQ=255)
                                                                                                          aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGa                                              >  2:814146/1‑91 (MQ=255)
                                                                                                          aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGa                                              <  1:814146/91‑1 (MQ=255)
                                                                                                          aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACg   >  2:727482/1‑134 (MQ=255)
                                                                                                           aGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGc  >  1:954680/1‑134 (MQ=255)
                                                                                                                      |                                                                                                                           
CGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGC  >  AM260480/2843463‑2843702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ 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
CCACCGGGCCCGCCGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGCTCCGCGGCGCGCTCGCTGTAGACAAG  >  AM260480/2843434‑2843728
                                                                                                                                                   |                                                                                                                                                    
gagatgtgtataagagaCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGgat                                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1210:2698:41342/133‑4 (MQ=60)
CCACCGGGCCCGCCGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGgat                                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2209:18331:78995/150‑4 (MQ=60)
 gtataagagacagGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATA                                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2212:11677:54136/137‑1 (MQ=60)
 agaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATA                                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2209:12412:7133/89‑1 (MQ=60)
 CACCGGGCCCGCCGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATA                                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2215:10320:66437/150‑1 (MQ=60)
     gatttttttgtgttttttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGG                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2205:11612:24288/57‑1 (MQ=60)
     GGGCCCGCCGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGG                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1102:7435:87141/150‑1 (MQ=60)
      ttgttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGG                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1109:3958:90141/70‑1 (MQ=60)
             gctttatctttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggcttcggagatgtgtataagagacaGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGG                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1104:10879:39439/64‑1 (MQ=60)
              GTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGG                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2201:10593:50114/1‑148 (MQ=60)
                    ggctcggagatgtgtataagagacaGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2102:13886:23929/125‑1 (MQ=60)
                    ggctcggagatgtgtataagagacaGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGACCGGG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1213:17377:61169/125‑1 (MQ=60)
                         CCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCC                                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2216:5911:9780/1‑149 (MQ=60)
                               gctcggagatgtgtataagagacaGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGT                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2115:18839:46863/126‑1 (MQ=60)
                                    cagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCC                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1114:16274:8941/88‑1 (MQ=60)
                                      agaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGG                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1114:7428:93698/90‑1 (MQ=60)
                                                   ataagagacagGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAA                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2104:19751:82047/138‑1 (MQ=60)
                                                    atgtgtataagagacaGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAAC                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1203:12852:64450/134‑1 (MQ=60)
                                                            TAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTctgtc                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1204:16789:101155/1‑145 (MQ=60)
                                                              GTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACctgtctcttatacacatctccgagcccac                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1116:3349:16019/1‑121 (MQ=60)
                                                              GTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGctgtctcttatacacatctccg                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2101:21255:99204/1‑126 (MQ=60)
                                                                      gtgggctcggagatgtgtataagagacaGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAAC                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2211:10594:80440/122‑1 (MQ=60)
                                                                                GCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCATAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACAC                                                                    >  GWNJ‑0478:712:GW2002102894th:2:2201:12982:52256/1‑150 (MQ=60)
                                                                                   GACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgc                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2103:3741:19546/1‑97 (MQ=60)
                                                                                   GACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGgactgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaatatagg                                                                >  GWNJ‑0478:712:GW2002102894th:2:1207:3310:2310/1‑64 (MQ=60)
                                                                                           gtgggctcggagatgtgtataagagacaGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCC                                                         <  GWNJ‑0478:712:GW2002102894th:2:2115:4263:91294/122‑1 (MQ=60)
                                                                                                agtttagcgagtgtctcgtgggctcggagatgtgtataagagacaggcatgGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCC                                                    <  GWNJ‑0478:712:GW2002102894th:2:1204:20925:18786/99‑1 (MQ=60)
                                                                                                  acgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCC                                                   <  GWNJ‑0478:712:GW2002102894th:2:1206:9856:15469/100‑1 (MQ=60)
                                                                                                   GTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaa                                                  >  GWNJ‑0478:712:GW2002102894th:2:2115:15978:9926/1‑81 (MQ=60)
                                                                                                   GTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgc                                                 >  GWNJ‑0478:712:GW2002102894th:2:1206:3103:89154/1‑87 (MQ=60)
                                                                                                   GTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCtgtctcttatacacatctccgagcccacgagacactcgc                                                 >  GWNJ‑0478:712:GW2002102894th:2:1107:10181:71624/1‑111 (MQ=60)
                                                                                                   GTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCctgtctct                                                 >  GWNJ‑0478:712:GW2002102894th:2:1209:5248:47790/1‑142 (MQ=60)
                                                                                                   GTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTC                                                 >  GWNJ‑0478:712:GW2002102894th:2:1209:20387:48706/1‑150 (MQ=60)
                                                                                                                CGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGgatctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaacaccgactaaacacatgaacaggccaacga                                    >  GWNJ‑0478:712:GW2002102894th:2:1114:9258:50834/1‑35 (MQ=60)
                                                                                                                     TCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCtgtctcttatacacatctccgagcccacgagacactcgctaatctc                               >  GWNJ‑0478:712:GW2002102894th:2:2214:15562:41301/1‑104 (MQ=60)
                                                                                                                       GGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCctgtctcttatacacatctccgagccc                              >  GWNJ‑0478:712:GW2002102894th:2:1212:12602:34456/1‑122 (MQ=60)
                                                                                                                         GATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGCT                           >  GWNJ‑0478:712:GW2002102894th:2:2206:15161:12110/1‑150 (MQ=60)
                                                                                                                                          GCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGCTCCGCGGCGCGCTCGCTG          >  GWNJ‑0478:712:GW2002102894th:2:1212:9878:72136/1‑150 (MQ=60)
                                                                                                                                            CGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgt        >  GWNJ‑0478:712:GW2002102894th:2:2111:19872:41862/1‑94 (MQ=60)
                                                                                                                                                 GGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaa   >  GWNJ‑0478:712:GW2002102894th:2:1215:2032:89984/1‑75 (MQ=60)
                                                                                                                                                 GGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCctgtctcttatacacatctccgagcccacgagacactcgctaat  >  GWNJ‑0478:712:GW2002102894th:2:2213:3987:34217/1‑106 (MQ=60)
                                                                                                                                                   |                                                                                                                                                    
CCACCGGGCCCGCCGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACGCTCCGCGGCGCGCTCGCTGTAGACAAG  >  AM260480/2843434‑2843728

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