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 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% 98.7 / NA 30intergenic (‑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 (13/17);  total (13/17)

AGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTC  >  AM260479/735389‑735647
                                                                                                                                     |                                                                                                                              
aGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAg                                                                                                                               >  2:854461/1‑134 (MQ=255)
    aCAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAgggg                                                                                                                            >  2:706650/1‑132 (MQ=255)
    aCAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAgggg                                                                                                                            >  1:729207/1‑132 (MQ=255)
     cAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAggggca                                                                                                                          <  2:478999/134‑4 (MQ=255)
     cAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAggggca                                                                                                                          <  1:706650/134‑4 (MQ=255)
           gTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCt                                                                                                                     <  1:785493/134‑1 (MQ=255)
                cAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                <  2:196846/134‑1 (MQ=255)
                cAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                <  2:729207/134‑1 (MQ=255)
                 aaaaaaCGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAggggcag                                                                                                                         >  1:195256/1‑119 (MQ=255)
                 aaaaaaCGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAggggcag                                                                                                                         <  2:195256/123‑5 (MQ=255)
                    aaaCGCGTCGCCTGAGGGGCTGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGc                                                                                                                   <  1:51161/127‑1 (MQ=255)
                    aaaCGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGc                                                                                                                   >  2:51161/1‑127 (MQ=255)
                             gCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGa                                                                                                   <  2:376306/134‑1 (MQ=255)
                                       ccggcAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCGGTGCCTGCGTTAACCCACTATTGAACGATCCCCa                                                                                         <  1:329778/134‑1 (MQ=255)
                                          gcAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                <  1:565989/108‑1 (MQ=255)
                                          gcAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                >  2:565989/1‑108 (MQ=255)
                                                                tggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCt                                                                                                                     <  1:433793/81‑1 (MQ=255)
                                                                tggtTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCt                                                                                                                     >  2:433793/1‑81 (MQ=255)
                                                                        aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCAt                                                        >  1:456747/1‑134 (MQ=255)
                                                                        aCGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCAt                                                        >  1:446403/1‑134 (MQ=255)
                                                                          ggtggtGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATcg                                                      <  2:446403/134‑1 (MQ=255)
                                                                           gtggtgCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATcgc                                                     <  1:854461/134‑1 (MQ=255)
                                                                                   tGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                <  2:258483/67‑1 (MQ=255)
                                                                                   tGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                >  1:258483/1‑67 (MQ=255)
                                                                                          ccGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAtt                                                                                                     >  1:86916/1‑71 (MQ=255)
                                                                                          ccGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAtt                                                                                                     <  2:86916/71‑1 (MQ=255)
                                                                                          ccGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGgc                                      >  2:888198/1‑134 (MQ=255)
                                                                                          ccGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGgc                                      >  1:74532/1‑134 (MQ=255)
                                                                                            ggCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAtt                                                                                                     >  2:513094/1‑69 (MQ=0)
                                                                                            ggCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAtt                                                                                                     <  1:513094/69‑1 (MQ=0)
                                                                                                 ggTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                >  2:412484/1‑53 (MQ=11)
                                                                                                 ggTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGtt                                                                                                                <  1:412484/53‑1 (MQ=11)
                                                                                                              aCCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCctgctg                                                                  >  1:880808/1‑86 (MQ=255)
                                                                                                              aCCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCctgctg                                                                  <  2:880808/86‑1 (MQ=255)
                                                                                                               cccTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGc                 <  1:888198/134‑1 (MQ=255)
                                                                                                                        gCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGt        >  1:325692/1‑134 (MQ=255)
                                                                                                                              tgAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTttcttc  <  2:325692/134‑1 (MQ=255)
                                                                                                                                     |                                                                                                                              
AGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTC  >  AM260479/735389‑735647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ 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
CTTTTTGTCCGAGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTCTGCGCTTGCGCGGGGCG  >  AM260479/735378‑735664
                                                                                                                                                   |                                                                                                                                            
CTTTTTGTCCGAGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGgc                                                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2115:6762:76773/1‑147 (MQ=60)
     TGTCCGAGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGC                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2113:16358:75401/149‑1 (MQ=60)
             CGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTA                                                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1112:3937:52568/1‑149 (MQ=60)
               ACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACC                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2112:19336:44524/150‑1 (MQ=60)
               ACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGgcctgtctcttatacac                                                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1105:3723:89719/1‑132 (MQ=60)
                     gctccgtctcgtgggctcggagatgtgtataagagacagATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAT                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1213:19413:45542/111‑1 (MQ=55)
                     AGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTAT                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2203:14206:33376/150‑1 (MQ=60)
                      GTAACCAAAAACGCGTCGCCTGAGGGGCTGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTG                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1102:4230:50903/150‑1 (MQ=60)
                           caagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAAC                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1210:15467:25199/83‑1 (MQ=31)
                           ataagagacaGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAAC                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2105:4673:65226/139‑1 (MQ=60)
                                     atacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGC                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1216:18758:94521/99‑1 (MQ=41)
                                                 GCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCGGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCC                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1207:17736:21161/150‑1 (MQ=60)
                                                                    GTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATC                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1212:1925:75455/1‑150 (MQ=60)
                                                                    GTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATC                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1212:10215:16018/1‑150 (MQ=60)
                                                                               TCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgc                                                             >  GWNJ‑0478:712:GW2002102894th:2:1201:13496:42496/1‑97 (MQ=41)
                                                                                     GGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGC                                                       <  GWNJ‑0478:712:GW2002102894th:2:2209:21105:81941/150‑1 (MQ=60)
                                                                                      GTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgt                                                      >  GWNJ‑0478:712:GW2002102894th:2:1107:17580:86797/1‑101 (MQ=45)
                                                                                      GTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCT                                                      >  GWNJ‑0478:712:GW2002102894th:2:1107:4195:8036/1‑150 (MQ=60)
                                                                                                          GTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGctgtctcttatacacatctccgagcccacgagac                                  >  GWNJ‑0478:712:GW2002102894th:2:2207:18866:47862/1‑116 (MQ=59)
                                                                                                                    GTTCGACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTT                        >  GWNJ‑0478:712:GW2002102894th:2:1204:12458:92927/1‑150 (MQ=60)
                                                                                                                         ACCCTGCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTC                   <  GWNJ‑0478:712:GW2002102894th:2:2207:10807:93028/150‑1 (MQ=60)
                                                                                                                              GCTACGCCATGTGAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaacggtg               >  GWNJ‑0478:712:GW2002102894th:2:1101:13332:16485/1‑63 (MQ=1)
                                                                                                                                          GAGTCAGGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTCTGCGCTTGCGCGGGGCG  >  GWNJ‑0478:712:GW2002102894th:2:2202:16958:40607/1‑150 (MQ=60)
                                                                                                                                                   |                                                                                                                                            
CTTTTTGTCCGAGCGACAAAAAGTAACCAAAAAACGCGTCGCCTGAGCGGCCGGCAATCAATTCCAGGGTGTTGGTGGTTCGCACGGTGGTGCTTGCCGTTCCGGCGTGGTTGCCCGTTCGACCCTGCTACGCCATGTGAGTCAGGGCAGTGCCTGCGTTAACCCACTATTGAACGATCCCCATGCAGGGCGCAGGCCGCCTGCTGCGAGCGGCATCGCGGGGACGCCTACGGCTGCTGCGCGGGCTCAAAGTGCATCGTCGGTTTCTTCTGCGCTTGCGCGGGGCG  >  AM260479/735378‑735664

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: