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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,515,428 +C intergenic (+367/‑456) tRNA‑Pro (CGG) → / → pilB tRNA‑Pro/Type IV pilus assembly protein PilB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,515,4271.C96.8% 112.5 / ‑4.2 31intergenic (+366/‑457)tRNA‑Pro (CGG)/pilBtRNA‑Pro/Type IV pilus assembly protein PilB
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (11/19);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCAT  >  AM260479/3515296‑3515549
                                                                                                                                    |                                                                                                                          
cgcgGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGcct                                                                                                                         >  1:346213/1‑132 (MQ=255)
 gcggcgGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGcctt                                                                                                                        <  1:183305/134‑4 (MQ=255)
 gcggcgGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGcctt                                                                                                                        <  1:997211/134‑4 (MQ=255)
           aaTCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtct                                                                                                               >  1:429169/1‑134 (MQ=255)
                  cttCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtcttc                                                                                                             <  2:411266/129‑1 (MQ=255)
                          aGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCatgat                                                                                                 <  1:67496/134‑1 (MQ=255)
                            gTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATgg                                                                                               >  2:466521/1‑134 (MQ=255)
                            gTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATgg                                                                                               <  2:1023389/134‑1 (MQ=255)
                               gCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCa                                                                                            <  2:429169/134‑1 (MQ=255)
                                   aaGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATccc                                                                                        <  2:346213/134‑1 (MQ=255)
                                             cGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAACCTGGCGGGCCGCCTGCCCGGGCATCTTCTTCCCGGCCTCAACCGGGCTTCTGGGGTTAAGGCTGGCATCCCTTtctctcgt                                                                              >  1:218857/1‑131 (MQ=255)
                                             cGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtcttc                                                                                                             <  1:884432/103‑1 (MQ=255)
                                             cGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtcttc                                                                                                             >  2:884432/1‑103 (MQ=255)
                                             cGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtct                                                                                                               >  2:794744/1‑101 (MQ=255)
                                             cGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtct                                                                                                               <  1:794744/101‑1 (MQ=255)
                                                   agaAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGt                                                                                              <  1:840412/112‑1 (MQ=255)
                                                   agaAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGt                                                                                              >  2:840412/1‑112 (MQ=255)
                                                       gTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGAt                                                                    <  2:31862/134‑1 (MQ=255)
                                                                      aGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCa                                                                                            <  2:355587/95‑1 (MQ=255)
                                                                      aGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCa                                                                                            >  1:355587/1‑95 (MQ=255)
                                                                        ccTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCt                                                   <  1:442771/134‑1 (MQ=255)
                                                                           gCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTtgt                                                >  1:94329/1‑134 (MQ=255)
                                                                               gCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTgc                                              <  2:94329/132‑1 (MQ=255)
                                                                               gCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTgc                                              <  1:93995/132‑1 (MQ=255)
                                                                                      ggCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTtg                                     <  2:736990/134‑1 (MQ=255)
                                                                                        gctGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTtgat                                   <  2:218857/132‑1 (MQ=255)
                                                                                              aaaCTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTgcgc                                            >  2:839060/1‑119 (MQ=255)
                                                                                              aaaCTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTgcgc                                            <  1:839060/119‑1 (MQ=255)
                                                                                                 cTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtc                                                                                                                <  2:552883/48‑1 (MQ=39)
                                                                                                 cTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGtc                                                                                                                >  1:552883/1‑48 (MQ=39)
                                                                                                                    ggCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAAtt       <  2:795033/134‑1 (MQ=255)
                                                                                                                        tCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCa   >  1:428911/1‑134 (MQ=255)
                                                                                                                         cTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCAt  >  2:662083/1‑134 (MQ=255)
                                                                                                                            gCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCAt  <  2:428911/131‑1 (MQ=255)
                                                                                                                                    |                                                                                                                          
CGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCAT  >  AM260479/3515296‑3515549

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ 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
GTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCATCACGCAGTGTTATGCACGCACTTTTG  >  AM260479/3515281‑3515575
                                                                                                                                                    |                                                                                                                                                   
GTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCctt                                                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1206:7969:24482/1‑147 (MQ=60)
         CGCCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2109:8419:100914/150‑1 (MQ=60)
           CCTCCGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTT                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1213:2359:20948/1‑150 (MQ=60)
               CGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGC                                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2215:3110:30835/150‑1 (MQ=60)
               CGCGGCGGTGAAATCTGCTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGC                                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1113:15596:80484/150‑1 (MQ=60)
                          agatgtgtataagagacagGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATG                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2203:5269:94975/131‑1 (MQ=60)
                            atgtgtataagagacagGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGT                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2207:11847:69769/133‑1 (MQ=60)
                                        GAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1104:10736:61255/150‑1 (MQ=60)
                                           agagacagAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCT                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2205:16929:96290/142‑1 (MQ=60)
                                             GGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCcttcatcctgtttcttatacacatctccgagcccacgagacggagct                                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1111:2592:75586/1‑103 (MQ=60)
                                               acgtatctattattttttcttgttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATC                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1101:5799:18693/51‑1 (MQ=60)
                                                              taacttgtgttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagattgtataagagacaggcCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGG                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1215:9590:59253/63‑1 (MQ=60)
                                                                      GTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCctgtctcttatacacatctccgagc                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1206:6159:83649/1‑125 (MQ=60)
                                                                           CGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTG                                                                         >  GWNJ‑0478:712:GW2002102894th:2:1106:12864:34375/1‑150 (MQ=60)
                                                                                      GCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGA                                                              <  GWNJ‑0478:712:GW2002102894th:2:1213:20739:94941/150‑1 (MQ=60)
                                                                                             GGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTG                                                       <  GWNJ‑0478:712:GW2002102894th:2:2205:7964:87072/150‑1 (MQ=60)
                                                                                             CGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAAC                                                         <  GWNJ‑0478:712:GW2002102894th:2:1106:4010:32276/148‑1 (MQ=60)
                                                                                                 CACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaa                                                    >  GWNJ‑0478:712:GW2002102894th:2:2102:7038:75259/1‑78 (MQ=60)
                                                                                                                        GCTCCTGCGCGGGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCAT                            >  GWNJ‑0478:712:GW2002102894th:2:1213:12363:19258/1‑150 (MQ=60)
                                                                                                                                   GGCATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCATCACGCAGTGTC                 >  GWNJ‑0478:712:GW2002102894th:2:1208:6714:79559/1‑150 (MQ=60)
                                                                                                                                      ATCTTGCTGACGGCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTctgtctcttatacacatct              >  GWNJ‑0478:712:GW2002102894th:2:2213:19386:58417/1‑131 (MQ=60)
                                                                                                                                                 GCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCATCACGCAGTGTTATGCACGCACTTTTG  >  GWNJ‑0478:712:GW2002102894th:2:2107:19744:86737/1‑150 (MQ=60)
                                                                                                                                                 GCCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCATCACGCAGTGTTATGCACGCACTTTTG  >  GWNJ‑0478:712:GW2002102894th:2:2213:15577:15773/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                                   
GTATTCGGCCGCCTCCGCGGCGGTGAAATCTGCGTTCCGCGAGGTGGCCAAAGACTTCCTCGTGCGAGAAGTCCTCGTCCTTGAAAGCCTGCTCGCCCACCGGCGTGAAAAACTTGATGCGCTCCTGCGCGGGCATCTTGCTGACGGCTTCATCAGGTCTTCTGCGGTCATGATGGTCATCCCTTTCTCTCTTCATCCTCGATCGCCAGGTGCAGCGGCTTGTGCGCCTTCTTGATGAACTGCATGAAGGGCTTCCTGCAAATTGGCATCACGCAGTGTTATGCACGCACTTTTG  >  AM260479/3515281‑3515575

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: