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 208,932 A→C H322Q (CAT→CAG h16_A0195 ← Putative GTPase (G3E family)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479208,9320AC100.0% 148.6 / NA 38H322Q (CAT→CAGh16_A0195Putative GTPase (G3E family)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/20);  total (18/20)

GGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATA  >  AM260479/208805‑209005
                                                                                                                               |                                                                         
ggcAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGaa                                                                      >  1:817580/1‑133 (MQ=255)
ggcAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGaa                                                                      >  1:797344/1‑133 (MQ=255)
ggcAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGaa                                                                      >  2:363615/1‑133 (MQ=255)
    gTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                 <  2:548119/134‑1 (MQ=255)
    gTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                 >  1:548119/1‑134 (MQ=255)
     ttCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcaccac                                                                <  2:116944/134‑1 (MQ=255)
          cAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGc                                                           >  1:245750/1‑134 (MQ=255)
               ccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  2:797344/134‑2 (MQ=255)
               ccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  2:817580/134‑2 (MQ=255)
               ccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  1:363615/134‑2 (MQ=255)
               ccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  1:30305/134‑2 (MQ=255)
                     ggATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCg                                                >  1:887257/1‑134 (MQ=255)
                      gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                 >  2:749090/1‑132 (MQ=255)
                      gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                 <  1:749090/132‑1 (MQ=255)
                               gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                       <  1:523914/133‑1 (MQ=255)
                               gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                       <  2:245750/133‑1 (MQ=255)
                                  gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  2:126406/115‑2 (MQ=255)
                                  gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      >  1:126406/1‑114 (MQ=255)
                                       ggTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGccc                              >  1:819251/1‑134 (MQ=255)
                                                      acacCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGt               >  2:915049/1‑134 (MQ=255)
                                                              cgcgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTc                                                       >  1:663331/1‑86 (MQ=255)
                                                              cgcgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTc                                                       <  2:663331/86‑1 (MQ=255)
                                                              cgcgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                 >  2:476007/1‑92 (MQ=255)
                                                              cgcgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                 <  1:476007/92‑1 (MQ=255)
                                                              cgcgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTccc        <  1:915049/133‑1 (MQ=255)
                                                                    cccGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGtata  <  1:467968/133‑1 (MQ=255)
                                                                        gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  1:852987/77‑2 (MQ=255)
                                                                        gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      >  2:852987/1‑76 (MQ=255)
                                                                                             cgcCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                 <  1:887662/45‑1 (MQ=255)
                                                                                             cgcCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                 >  2:887662/1‑45 (MQ=255)
                                                                                               ccGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      >  1:67229/1‑53 (MQ=255)
                                                                                               ccGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  2:67229/54‑2 (MQ=255)
                                                                                                 gATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTc                                              <  1:215041/60‑1 (MQ=255)
                                                                                                 gATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTc                                              >  2:215041/1‑60 (MQ=255)
                                                                                                              aTCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      <  1:972834/39‑2 (MQ=255)
                                                                                                              aTCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                      >  2:972834/1‑38 (MQ=255)
                                                                                                                 aGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGtata  >  2:678404/1‑88 (MQ=255)
                                                                                                                 aGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGtata  <  1:678404/88‑1 (MQ=255)
                                                                                                                               |                                                                         
GGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATA  >  AM260479/208805‑209005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
GGTGCACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTCGCTGCGGAAGAC  >  AM260479/208790‑209072
                                                                                                                                              |                                                                                                                                            
GGTGCACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAC                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2204:7863:9902/1‑149 (MQ=60)
GGTGCACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAC                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2203:9425:69407/1‑149 (MQ=60)
    CACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCAC                                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2102:12197:8142/1‑150 (MQ=60)
         AGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGC                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1116:5631:12534/150‑1 (MQ=60)
          GGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCG                                                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1116:13644:38812/1‑150 (MQ=60)
              AGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACctgtctctt                                                                                                                          >  GWNJ‑0478:712:GW2002102894th:2:2111:4025:71617/1‑140 (MQ=60)
                   ttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAC                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2207:17308:56681/75‑1 (MQ=60)
                     TCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGT                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2207:2795:54079/1‑150 (MQ=60)
                             CCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATAT                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1105:19537:13269/150‑1 (MQ=60)
                             CCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATAT                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1208:11530:24668/150‑1 (MQ=60)
                              ttctttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATAT                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2213:12284:61068/70‑1 (MQ=60)
                              gtagctccgtctcgtgggctcggagatgtgtataagagacagCCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATAT                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2209:3660:39633/107‑1 (MQ=60)
                                  GAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATctgtc                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1109:10196:44478/1‑145 (MQ=60)
                                  GAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTAC                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2114:4977:60465/150‑1 (MQ=60)
                                   gggctcggagatgtgtataagagacaGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTAC                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1215:9380:65027/123‑1 (MQ=60)
                                     gaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCC                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1115:17809:63225/90‑1 (MQ=60)
                                       TGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCT                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2204:13846:21186/1‑150 (MQ=60)
                                             CGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAG                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2101:12242:38774/149‑1 (MQ=60)
                                                              CCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATActgtctcttatacacatctccgagcccacgagac                                                                         >  GWNJ‑0478:712:GW2002102894th:2:2111:18326:65747/1‑116 (MQ=60)
                                                                            GCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGA                                                            <  GWNJ‑0478:712:GW2002102894th:2:2211:5842:27371/149‑1 (MQ=60)
                                                                                  GCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGG                                                      <  GWNJ‑0478:712:GW2002102894th:2:1215:4368:23974/149‑1 (MQ=60)
                                                                                   cgtgggctcggagatgtgtataagagacaGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGG                                                      <  GWNJ‑0478:712:GW2002102894th:2:2110:18934:55177/119‑1 (MQ=60)
                                                                                               TCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaatcaaaaaaac                                         >  GWNJ‑0478:712:GW2002102894th:2:1203:14762:66053/1‑59 (MQ=60)
                                                                                               TCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttg                                        >  GWNJ‑0478:712:GW2002102894th:2:1104:18550:65073/1‑84 (MQ=60)
                                                                                                                         GCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaccaaagagacgat              >  GWNJ‑0478:712:GW2002102894th:2:2115:14586:83241/1‑61 (MQ=60)
                                                                                                                                 GCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTCGCTGCGG       >  GWNJ‑0478:712:GW2002102894th:2:1208:16880:51596/1‑149 (MQ=60)
                                                                                                                                     GTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCtgtctcttatacacatctccgagcccacgagacggagct  >  GWNJ‑0478:712:GW2002102894th:2:2113:9849:53177/1‑111 (MQ=60)
                                                                                                                                              |                                                                                                                                            
GGTGCACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTCGCTGCGGAAGAC  >  AM260479/208790‑209072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: