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 735,524 +G intergenic (‑409/+175) h16_A0685 ← / ← h16_A0686 hypothetical membrane spanning protein/Predictedpermease/CBS‑domain‑containing membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479735,5211.G100.0% 117.5 / NA 31intergenic (‑406/+178)h16_A0685/h16_A0686hypothetical membrane spanning protein/Predictedpermease/CBS‑domain‑containing membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/15);  total (16/15)

ACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTC  >  AM260479/735393‑735638
                                                                                                                                 |                                                                                                                     
aCAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAgggg                                                                                                                   >  2:100865/1‑132 (MQ=255)
  aaaaaGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAggggcag                                                                                                                <  1:1119201/134‑5 (MQ=255)
       gTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCt                                                                                                            <  1:581159/134‑1 (MQ=255)
            cAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                       <  2:127270/134‑1 (MQ=255)
                      gTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAt                                                                                             <  2:1091307/134‑1 (MQ=255)
                         gCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGa                                                                                          <  2:265601/134‑1 (MQ=255)
                         gCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGa                                                                                          <  1:252303/134‑1 (MQ=255)
                              aGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACg                                                                                       <  1:703509/132‑1 (MQ=255)
                              aGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACg                                                                                       <  1:596714/132‑1 (MQ=255)
                              aGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGAt                                                                                     >  2:244784/1‑134 (MQ=255)
                                  gccggcAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGACTGCGTTAACCCACTATTGAACGATccc                                                                                  >  2:560381/1‑133 (MQ=255)
                                                          ggtggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCa                                                                                                                        <  1:67981/71‑1 (MQ=35)
                                                          ggtggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCa                                                                                                                        >  2:67981/1‑71 (MQ=35)
                                                            tggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGa                                                                                          >  2:273902/1‑99 (MQ=255)
                                                            tggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGa                                                                                          <  1:273902/99‑1 (MQ=255)
                                                            tggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGgccgcc                                                              <  1:478252/127‑1 (MQ=255)
                                                            tggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGgccgcc                                                              >  2:478252/1‑127 (MQ=255)
                                                                  gCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCgg                                                  >  2:388449/1‑133 (MQ=255)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCa                                                                                                                        <  1:770273/61‑1 (MQ=37)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCa                                                                                                                        >  2:770273/1‑61 (MQ=37)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGc                                                                                                          <  2:1099869/75‑1 (MQ=255)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGc                                                                                                          >  1:1099869/1‑75 (MQ=255)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCAt                                               >  2:1073840/1‑134 (MQ=255)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCAt                                               >  2:192880/1‑134 (MQ=255)
                                                                    aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCAt                                               >  1:560775/1‑134 (MQ=255)
                                                                            gCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGa                                       >  1:27822/1‑134 (MQ=255)
                                                                                         gCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGgctgc                          <  1:388449/134‑1 (MQ=255)
                                                                                         gCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGgctgc                          <  1:1073840/134‑1 (MQ=255)
                                                                                           gTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGctgctg                        <  2:167752/134‑1 (MQ=255)
                                                                                                tGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCgg                   >  1:246766/1‑134 (MQ=255)
                                                                                                            ccTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCa       >  1:623071/1‑134 (MQ=255)
                                                                                                                cTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCAtcgt   >  2:962074/1‑134 (MQ=255)
                                                                                                                 tACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCAtcgtc  >  1:312392/1‑134 (MQ=255)
                                                                                                                    gCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGgctgct                         >  1:754948/1‑108 (MQ=255)
                                                                                                                    gCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGgctgct                         <  2:754948/108‑1 (MQ=255)
                                                                                                                            aGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGa                                       <  1:627567/86‑1 (MQ=255)
                                                                                                                            aGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGa                                       >  2:627567/1‑86 (MQ=255)
                                                                                                                                 |                                                                                                                     
ACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTC  >  AM260479/735393‑735638

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
CGAGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTCTGCGCTTGCGCGGGGCGCAGCC  >  AM260479/735387‑735669
                                                                                                                                          |                                                                                                                                                 
acgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCG                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1104:8889:34461/101‑1 (MQ=43)
aagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCG                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2113:3159:27032/92‑1 (MQ=22)
       CAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCC                                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2214:14702:84402/150‑1 (MQ=60)
        agcgagtgtctcgtgggctcggagatgtgtataagagacagCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCC                                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2204:9683:17997/108‑1 (MQ=58)
            AGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAT                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2101:5718:43203/150‑1 (MQ=60)
                              ggagatgtgtataagagacagATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCGCATGCAG                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1116:16451:47805/129‑1 (MQ=60)
                              CGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAG                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1114:7588:42678/150‑1 (MQ=60)
                                   GAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGC                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2108:1513:76841/148‑1 (MQ=60)
                                   GAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGC                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2103:20464:41051/148‑1 (MQ=60)
                                                          GGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCAT                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1213:5658:40264/150‑1 (MQ=60)
                                                           GTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGctgtctcttatacacatctccgagcccacgagacactcgctaatc                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2213:10888:93809/1‑105 (MQ=40)
                                                           GTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATC                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1216:5250:34868/1‑150 (MQ=60)
                                                                   GGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGAC                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1103:1519:79622/1‑150 (MQ=60)
                                                                        agcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagaCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTAC                                                               <  GWNJ‑0478:712:GW2002102894th:2:1201:10543:77737/89‑1 (MQ=60)
                                                                                 gtctcgtgggctcggagatgtgtataagagacagCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTG                                                       <  GWNJ‑0478:712:GW2002102894th:2:2105:7912:35070/116‑1 (MQ=60)
                                                                                       CCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGG                                                 >  GWNJ‑0478:712:GW2002102894th:2:1114:10205:14727/1‑150 (MQ=60)
                                                                                              GGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAG                                          <  GWNJ‑0478:712:GW2002102894th:2:1206:5652:50373/150‑1 (MQ=60)
                                                                                              GGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAG                                          <  GWNJ‑0478:712:GW2002102894th:2:2212:20275:68303/150‑1 (MQ=60)
                                                                                               cgtgggctcggagatgtgtataagagacagATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAG                                          <  GWNJ‑0478:712:GW2002102894th:2:2113:18701:63984/119‑1 (MQ=60)
                                                                                                   GGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCAT                                     >  GWNJ‑0478:712:GW2002102894th:2:2105:9575:9265/1‑150 (MQ=60)
                                                                                                    gggctcggagatgtgtataagagacagGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCAT                                     <  GWNJ‑0478:712:GW2002102894th:2:1215:11087:14694/122‑1 (MQ=60)
                                                                                                        CCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCG                                >  GWNJ‑0478:712:GW2002102894th:2:1203:10923:80050/1‑150 (MQ=60)
                                                                                                           GTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGgctgtctcttat                              >  GWNJ‑0478:712:GW2002102894th:2:2110:16727:48240/1‑137 (MQ=60)
                                                                                                                                      GGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTCTGCGCTTGCGCGGGGCGCAGCC  >  GWNJ‑0478:712:GW2002102894th:2:2113:2928:80399/1‑149 (MQ=60)
                                                                                                                                          |                                                                                                                                                 
CGAGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTCTGCGCTTGCGCGGGGCGCAGCC  >  AM260479/735387‑735669

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