Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 724,768 +C coding (82/669 nt) lysE → lysine efflux permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479724,7661.C100.0% 122.4 / NA 32R27P (CGC→CCC) lysElysine efflux permease
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (15/17);  total (15/17)

GATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATG  >  AM260479/724639‑724894
                                                                                                                                |                                                                                                                                
gATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACgccgccctgg                                                                                                                            >  2:680292/1‑130 (MQ=255)
                gCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCgg                                                                                                             >  1:46205/1‑134 (MQ=255)
                       gTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTcc                                                                                                      >  2:444090/1‑134 (MQ=255)
                                  gCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGc                                                                                           <  2:924949/134‑1 (MQ=255)
                                         tGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGt                                                                                    >  1:636745/1‑134 (MQ=255)
                                             aaGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTccc                                                                                <  2:636745/134‑1 (MQ=255)
                                                         cATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCg                                                                    <  1:582117/134‑1 (MQ=255)
                                                                 cTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCg                                                                                                              <  2:47379/84‑1 (MQ=255)
                                                                 cTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCg                                                                                                              >  1:47379/1‑84 (MQ=255)
                                                                      cccGGCGTGCAATAATAGCCACCCAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAggg                                                      >  1:100465/1‑134 (MQ=255)
                                                                               cAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGCATGCTTTCGTGCTGCGCCAAGGGTTgcggcg                                              >  2:236740/1‑134 (MQ=255)
                                                                               cAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTgcggcg                                              <  1:680292/134‑1 (MQ=255)
                                                                                  taataGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           <  1:475357/134‑1 (MQ=255)
                                                                                  taataGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           <  2:46205/134‑1 (MQ=255)
                                                                                  taataGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           <  1:444090/134‑1 (MQ=255)
                                                                                   aataGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAgg                                                       <  1:55212/121‑1 (MQ=255)
                                                                                   aataGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAgg                                                       >  2:55212/1‑121 (MQ=255)
                                                                                            cccaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTc                                 <  2:100465/131‑1 (MQ=255)
                                                                                              aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           >  1:851283/1‑122 (MQ=255)
                                                                                              aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           <  2:851283/122‑1 (MQ=255)
                                                                                              aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTg                               >  1:574656/1‑134 (MQ=255)
                                                                                                          gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTtct                                                                                                                    <  1:1012843/37‑1 (MQ=37)
                                                                                                          gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTtct                                                                                                                    >  2:1012843/1‑37 (MQ=37)
                                                                                                          gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCttt                                                                                                        <  2:206696/49‑1 (MQ=38)
                                                                                                          gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCttt                                                                                                        >  1:206696/1‑49 (MQ=38)
                                                                                                          gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAg                                          <  1:663216/111‑1 (MQ=255)
                                                                                                          gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAg                                          >  2:663216/1‑111 (MQ=255)
                                                                                                                      aTGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCtt                                                                                                         <  1:504961/36‑1 (MQ=37)
                                                                                                                      aTGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCtt                                                                                                         >  2:504961/1‑36 (MQ=37)
                                                                                                                       tGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTcc      <  2:594266/134‑1 (MQ=255)
                                                                                                                          aCGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTcc      <  2:313348/131‑1 (MQ=255)
                                                                                                                          aCGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTcc      >  1:313348/1‑131 (MQ=255)
                                                                                                                          aCGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGAt   >  1:912625/1‑134 (MQ=255)
                                                                                                                           cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           >  2:305246/1‑93 (MQ=255)
                                                                                                                           cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                                           <  1:305246/93‑1 (MQ=255)
                                                                                                                           cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCt                                <  1:551897/104‑1 (MQ=255)
                                                                                                                           cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCt                                >  2:551897/1‑104 (MQ=255)
                                                                                                                           cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATg  >  2:331178/1‑134 (MQ=255)
                                                                                                                                |                                                                                                                                
GATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATG  >  AM260479/724639‑724894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
ATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATGCGCTG  >  AM260479/724640‑724899
                                                                                                                                 |                                                                                                                                   
ATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGT                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1102:3090:18494/1‑150 (MQ=60)
        ttttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCC                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2215:20873:80579/67‑1 (MQ=60)
                   tctttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCC                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1216:13414:20581/66‑1 (MQ=60)
                         GTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTT                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2109:6265:7043/1‑150 (MQ=60)
                                                 GGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCctgtctcttatacacatctccgagcccacgagacgg                                                                >  GWNJ‑0478:712:GW2002102894th:2:1102:17813:26000/1‑114 (MQ=60)
                                                      GGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGT                                                          >  GWNJ‑0478:712:GW2002102894th:2:1106:17440:79346/1‑150 (MQ=60)
                                                       GCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGT                                                          <  GWNJ‑0478:712:GW2002102894th:2:2104:8565:50868/150‑1 (MQ=60)
                                                              gacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCG                                                   <  GWNJ‑0478:712:GW2002102894th:2:1104:4334:50580/93‑1 (MQ=60)
                                                                    CCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCA                                             <  GWNJ‑0478:712:GW2002102894th:2:1102:14496:81285/150‑1 (MQ=60)
                                                                             GCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTG                                    <  GWNJ‑0478:712:GW2002102894th:2:2110:7778:66583/150‑1 (MQ=60)
                                                                              CAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGG                                   >  GWNJ‑0478:712:GW2002102894th:2:2209:6040:29123/1‑150 (MQ=60)
                                                                              CAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGG                                   >  GWNJ‑0478:712:GW2002102894th:2:2105:18016:98615/1‑150 (MQ=60)
                                                                                gggctcggagatgtgtataagagacagTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGT                                 <  GWNJ‑0478:712:GW2002102894th:2:1205:18886:48476/123‑1 (MQ=60)
                                                                                ATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaa                                  >  GWNJ‑0478:712:GW2002102894th:2:1116:7075:2390/1‑66 (MQ=60)
                                                                                ATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGT                                 <  GWNJ‑0478:712:GW2002102894th:2:1213:8212:91139/150‑1 (MQ=60)
                                                                                ATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGT                                 <  GWNJ‑0478:712:GW2002102894th:2:1215:12495:61745/150‑1 (MQ=60)
                                                                                 aagagacagCCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTC                                <  GWNJ‑0478:712:GW2002102894th:2:2111:16339:64770/141‑1 (MQ=60)
                                                                                          CCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaa                        >  GWNJ‑0478:712:GW2002102894th:2:1115:18449:16882/1‑79 (MQ=60)
                                                                                            gtgtataagagacagTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGT                      <  GWNJ‑0478:712:GW2002102894th:2:2102:8532:37879/134‑1 (MQ=60)
                                                                                                      agatgtgtataagagacagacgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCC           <  GWNJ‑0478:712:GW2002102894th:2:2110:13123:50734/124‑1 (MQ=60)
                                                                                                          TTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATG       >  GWNJ‑0478:712:GW2002102894th:2:1205:7314:100407/1‑150 (MQ=60)
                                                                                                          TTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATG       >  GWNJ‑0478:712:GW2002102894th:2:2211:13783:12577/1‑150 (MQ=60)
                                                                                                               CATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTctgtctcttatacacatctccgagcccacgagacgg    >  GWNJ‑0478:712:GW2002102894th:2:1110:10780:2512/1‑112 (MQ=60)
                                                                                                               CATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATGCGCTG  >  GWNJ‑0478:712:GW2002102894th:2:1108:15574:89434/1‑150 (MQ=60)
                                                                                                                                 |                                                                                                                                   
ATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTGTGCGCGCTTTCCGATGCGCTG  >  AM260479/724640‑724899

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: