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 149 26.7 2337422 87.5% 2045244 105.0

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% 164.5 / NA 42H322Q (CAT→CAGh16_A0195Putative GTPase (G3E family)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/21);  total (21/21)

GTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTA  >  AM260479/208809‑209045
                                                                                                                           |                                                                                                                 
gTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                                                         <  2:1058644/134‑1 (MQ=255)
 ttCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcaccac                                                                                                        >  2:648969/1‑134 (MQ=255)
     ccAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTg                                                                                                    >  1:215348/1‑134 (MQ=255)
      cAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGc                                                                                                   >  2:1061184/1‑134 (MQ=255)
          cccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTc                                                                                               <  1:648969/134‑1 (MQ=255)
           ccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTc                                                                                               >  1:121591/1‑133 (MQ=255)
               gAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                                                            <  2:215348/132‑4 (MQ=255)
               gAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                                                              >  2:108244/1‑129 (MQ=255)
               gAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                                                              <  1:108244/130‑2 (MQ=255)
                aGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                                                         <  1:98518/134‑1 (MQ=255)
                           gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCacac                                                                                          <  1:200820/122‑1 (MQ=255)
                           gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCacac                                                                                          >  2:200820/1‑122 (MQ=255)
                              gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                                                              <  2:483548/115‑2 (MQ=255)
                              gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                                                              >  1:483548/1‑114 (MQ=255)
                                     tCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCtt                                                                    <  1:1061184/134‑1 (MQ=255)
                                             gATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                                                         <  2:1000975/105‑1 (MQ=255)
                                             gATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                                                         >  1:1000975/1‑105 (MQ=255)
                                                       aTGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAac                                                                                                             >  2:1394/1‑75 (MQ=255)
                                                       aTGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAac                                                                                                             <  1:1394/75‑1 (MQ=255)
                                                            cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                                                              >  1:461890/1‑84 (MQ=255)
                                                            cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                                                              <  2:461890/85‑2 (MQ=255)
                                                                    gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                                                            <  2:310240/79‑4 (MQ=255)
                                                                    gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                                                            >  1:310240/1‑76 (MQ=255)
                                                                        ctctTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAac                                                                                                             <  2:1128861/58‑1 (MQ=255)
                                                                        ctctTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAac                                                                                                             >  1:1128861/1‑58 (MQ=255)
                                                                        ctctTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGa                                 >  2:82709/1‑134 (MQ=255)
                                                                              ctccCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGcagcag                                                        <  1:819824/105‑1 (MQ=255)
                                                                              ctccCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGcagcag                                                        >  2:819824/1‑105 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcaccac                                                                                                        >  1:593004/1‑52 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcaccac                                                                                                        <  2:593004/52‑1 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                                                         <  1:63407/67‑1 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                                                         >  2:63407/1‑67 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                                                               >  2:248116/1‑77 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                                                               <  1:248116/77‑1 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGttt                                                     >  2:1016062/1‑103 (MQ=255)
                                                                                   aCTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGttt                                                     <  1:1016062/103‑1 (MQ=255)
                                                                                           ccGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCac                                                                                            >  2:385744/1‑56 (MQ=255)
                                                                                           ccGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCac                                                                                            <  1:385744/56‑1 (MQ=255)
                                                                                                               cTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTa  >  2:1034512/1‑126 (MQ=255)
                                                                                                               cTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTa  <  1:1034512/126‑1 (MQ=255)
                                                                                                                     gACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatat                                                                                 >  1:671221/1‑41 (MQ=255)
                                                                                                                     gACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatat                                                                                 <  2:671221/41‑1 (MQ=255)
                                                                                                                           |                                                                                                                 
GTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTA  >  AM260479/208809‑209045

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGTGCACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTCGCTGCGGAAGACGAACGAC  >  AM260479/208789‑209079
                                                                                                                                               |                                                                                                                                                   
cggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAC                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2104:17448:3314/95‑1 (MQ=60)
cggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAC                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1209:5329:64197/95‑1 (MQ=60)
     gaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCAC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1102:11180:79680/90‑1 (MQ=60)
     GACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCAC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1202:8721:27811/150‑1 (MQ=60)
      agtgatttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCAC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2102:8377:3104/68‑1 (MQ=60)
          agaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGC                                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2214:4115:29308/89‑1 (MQ=60)
          AGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGC                                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1112:2268:58658/1‑150 (MQ=60)
               gattagcgagtgtctcgtgggctcggagatgtgtataagagacagACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCC                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1101:16657:10713/105‑1 (MQ=60)
                GGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCC                                                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1108:14599:51685/1‑149 (MQ=60)
                 GCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACctgtctcttat                                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1101:17148:35390/1‑138 (MQ=60)
                             GCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATA                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2104:19018:61689/150‑1 (MQ=60)
                              CCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATAT                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1106:6239:74014/150‑1 (MQ=60)
                                agcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAG                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1206:7835:33969/86‑1 (MQ=60)
                                   ggcatacgagattagcgagtgtctcgtgggctcgggatgtgtataagagacaGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTAC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1104:5837:6678/98‑1 (MQ=60)
                                   GAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATctgtc                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:1213:5588:66224/1‑145 (MQ=60)
                                   GAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTAC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1106:5761:16018/150‑1 (MQ=60)
                                   GAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTA                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1113:1779:81421/149‑1 (MQ=60)
                                              tagcgagtgtctcgtgggctcggagatgtgtataagagacaGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAG                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1114:16914:21426/108‑1 (MQ=60)
                                                  GGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACctgtctcttatacac                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2207:17793:83188/1‑135 (MQ=60)
                                                        GTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTT                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2212:11227:5309/150‑1 (MQ=60)
                                                                 GATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATctgtctcttatacacatctccgagcccacgagac                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1210:7356:56119/1‑115 (MQ=60)
                                                                    gtgtataagagacagGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATAC                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2115:2441:68155/135‑1 (MQ=60)
                                                                       atgtgtataagagacaGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGC                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2208:8141:67037/134‑1 (MQ=60)
                                                                         CCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGctgtctcttatacacatctccgagcccacgagacactcgct                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1203:10979:66979/1‑109 (MQ=60)
                                                                             GCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCcctgtctcttatacacatctccgagcccacgaaacactcgctaatctcgtatgccgtcttctg                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2215:14818:27567/1‑87 (MQ=60)
                                                                                        GGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttga                                                        >  GWNJ‑0478:712:GW2002102894th:2:2103:3265:18857/1‑82 (MQ=60)
                                                                                                                          GCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaa                      >  GWNJ‑0478:712:GW2002102894th:2:2107:13819:89289/1‑71 (MQ=60)
                                                                                                                          GCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTC                     <  GWNJ‑0478:712:GW2002102894th:2:2211:9682:66525/150‑1 (MQ=60)
                                                                                                                                  GCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTCGCTGCGG              >  GWNJ‑0478:712:GW2002102894th:2:1211:8085:70634/1‑149 (MQ=60)
                                                                                                                                             CCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaacaagagcacggagaaaaaccgcacaggcaccgcg  >  GWNJ‑0478:712:GW2002102894th:2:1105:1618:95155/1‑39 (MQ=60)
                                                                                                                                             CCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaggaaagatagaaacagcggtgatagagtagc    >  GWNJ‑0478:712:GW2002102894th:2:1210:7203:17688/1‑39 (MQ=60)
                                                                                                                                             CCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaaagaaagaagagcacgcagaggc  >  GWNJ‑0478:712:GW2002102894th:2:2109:3858:14259/1‑44 (MQ=60)
                                                                                                                                               |                                                                                                                                                   
TGGTGCACTCAGGCCAGGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAGGATGCCGGACAGGAATTCCTCCAGCTTGCCGTAGTGGAAGGGCTTGTCGCTGCGGAAGACGAACGAC  >  AM260479/208789‑209079

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