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 146 27.9 2327082 89.7% 2087392 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,240,057 +T intergenic (‑57/‑218) h16_A1136 ← / → grpE conserved hypothetical protein/molecular chaperone GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,240,0531.T100.0% 164.1 / NA 43intergenic (‑53/‑222)h16_A1136/grpEconserved hypothetical protein/molecular chaperone GrpE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (22/21);  total (22/21)

AAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGC  >  AM260479/1239939‑1240166
                                                                                                                   |                                                                                                                 
aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                  <  2:247042/133‑1 (MQ=255)
aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                  <  2:277784/133‑1 (MQ=255)
       tCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAAt                                                                                          >  2:827232/1‑134 (MQ=255)
           gACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCt                                                                                        >  1:785950/1‑132 (MQ=255)
           gACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCt                                                                                        >  2:341025/1‑132 (MQ=255)
             ccGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAtt                                                                                                         <  2:205854/113‑1 (MQ=255)
             ccGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCAtt                                                                                                         >  1:205854/1‑113 (MQ=255)
                    tGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGc                                                                              >  2:160950/1‑133 (MQ=255)
                     gCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGa                                                                            >  1:1083360/1‑134 (MQ=255)
                     gCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGa                                                                            >  2:490893/1‑134 (MQ=255)
                        gCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                         <  1:695458/134‑1 (MQ=255)
                                          aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                        <  1:248160/133‑1 (MQ=255)
                                          aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                        <  2:785950/133‑1 (MQ=255)
                                          aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGt                                                       >  1:206489/1‑134 (MQ=255)
                                          aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGt                                                       >  1:421312/1‑134 (MQ=255)
                                           cGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTc                                                      >  2:195515/1‑134 (MQ=255)
                                              cTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   <  1:1062353/134‑1 (MQ=255)
                                               tGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   <  1:178533/133‑1 (MQ=255)
                                               tGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   >  2:178533/1‑133 (MQ=255)
                                                        aTGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATg                                           <  1:938505/132‑1 (MQ=255)
                                                         tGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGcc                                         <  2:206489/133‑1 (MQ=255)
                                                         tGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGcc                                         <  1:827232/133‑1 (MQ=255)
                                                         tGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGcc                                         <  2:421312/133‑1 (MQ=255)
                                                                    ccTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATg                                           >  2:671765/1‑120 (MQ=255)
                                                                    ccTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATg                                           <  1:671765/120‑1 (MQ=255)
                                                                       tCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCg                          >  1:793891/1‑134 (MQ=255)
                                                                             agagGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGtt                    >  2:538253/1‑134 (MQ=255)
                                                                                    caAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGt             >  1:94518/1‑134 (MQ=255)
                                                                                     aaaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                         <  1:989852/73‑1 (MQ=255)
                                                                                     aaaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                         >  2:989852/1‑73 (MQ=255)
                                                                                     aaaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTc            <  1:195515/134‑1 (MQ=255)
                                                                                      aaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttat                                                                                                    <  1:503791/45‑1 (MQ=38)
                                                                                      aaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttat                                                                                                    >  2:503791/1‑45 (MQ=38)
                                                                                       aGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   <  2:446089/93‑1 (MQ=255)
                                                                                       aGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   >  1:446089/1‑93 (MQ=255)
                                                                                             aCCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgc    <  1:490893/134‑1 (MQ=255)
                                                                                               cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   >  1:537996/1‑85 (MQ=255)
                                                                                               cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                   <  2:537996/85‑1 (MQ=255)
                                                                                                aTCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgcgc  >  1:387666/1‑133 (MQ=255)
                                                                                                aTCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgcg   >  1:88636/1‑132 (MQ=255)
                                                                                                aTCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgcg   <  2:88636/132‑1 (MQ=255)
                                                                                                       caACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCata                                             >  1:750899/1‑83 (MQ=255)
                                                                                                       caACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCata                                             <  2:750899/83‑1 (MQ=255)
                                                                                                                   |                                                                                                                 
AAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGC  >  AM260479/1239939‑1240166

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
GAATGCGCCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGCCT  >  AM260479/1239917‑1240187
                                                                                                                                             |                                                                                                                                  
GAATGCGCCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTT                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1111:6640:3179/149‑1 (MQ=60)
                  GATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTT                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2113:13168:78913/1‑148 (MQ=60)
                   ttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAG                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1212:11378:20470/75‑1 (MQ=60)
                    TAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTgtctc                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1113:17069:87425/1‑145 (MQ=60)
                            GTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGAT                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2213:17604:18884/1‑150 (MQ=60)
                                             agcttagcgagtgtctcgtgggctcggaggatgtgtataagagacagTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGG                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2207:21027:20889/103‑1 (MQ=60)
                                             GGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGG                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2106:14498:93658/150‑1 (MQ=60)
                                                 ACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTC                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1208:3650:98729/1‑150 (MQ=60)
                                                 ACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTC                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1113:16806:90710/1‑150 (MQ=60)
                                                              gaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCC                                                              <  GWNJ‑0478:712:GW2002102894th:2:1214:21023:48838/90‑1 (MQ=60)
                                                               CACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCC                                                              <  GWNJ‑0478:712:GW2002102894th:2:1115:16309:100088/149‑1 (MQ=60)
                                                                   ACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCC                                                         <  GWNJ‑0478:712:GW2002102894th:2:1110:3864:45517/150‑1 (MQ=60)
                                                                   ACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCC                                                         <  GWNJ‑0478:712:GW2002102894th:2:2212:16719:12235/150‑1 (MQ=60)
                                                                           catacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAAC                                                 <  GWNJ‑0478:712:GW2002102894th:2:1107:10853:41147/98‑1 (MQ=60)
                                                                             TATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAAC                                                 <  GWNJ‑0478:712:GW2002102894th:2:2205:4231:26828/148‑1 (MQ=60)
                                                                             TATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAAC                                                 <  GWNJ‑0478:712:GW2002102894th:2:2107:11287:52587/148‑1 (MQ=60)
                                                                              tagcgagtgtctcgtgggctcggagatgtgtataagagacagCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCG                                               <  GWNJ‑0478:712:GW2002102894th:2:2101:11364:94157/107‑1 (MQ=60)
                                                                              ATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGT                                              >  GWNJ‑0478:712:GW2002102894th:2:2112:17868:21208/1‑150 (MQ=60)
                                                                              ATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCG                                               <  GWNJ‑0478:712:GW2002102894th:2:2115:5051:88083/149‑1 (MQ=60)
                                                                                           CTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTC                                 >  GWNJ‑0478:712:GW2002102894th:2:1107:18026:79239/1‑150 (MQ=60)
                                                                                              CCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCctgtctcttatacacatctccgagccc                              >  GWNJ‑0478:712:GW2002102894th:2:1211:11346:31526/1‑123 (MQ=60)
                                                                                                      GGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCctgtctcttatacacatctccgagcccacgagac                       >  GWNJ‑0478:712:GW2002102894th:2:1215:6457:84922/1‑115 (MQ=60)
                                                                                                       GCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCC                      >  GWNJ‑0478:712:GW2002102894th:2:1206:4241:9723/1‑149 (MQ=60)
                                                                                                       GCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCC                      >  GWNJ‑0478:712:GW2002102894th:2:1107:15251:46201/1‑149 (MQ=60)
                                                                                                          CAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTC                  <  GWNJ‑0478:712:GW2002102894th:2:1113:19269:35083/150‑1 (MQ=60)
                                                                                                              GATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTActgtctcttatacacatctccgagcccacgagacact              >  GWNJ‑0478:712:GW2002102894th:2:2111:3684:96254/1‑113 (MQ=60)
                                                                                                                  AACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAAC          <  GWNJ‑0478:712:GW2002102894th:2:1213:2303:57861/150‑1 (MQ=60)
                                                                                                                        CTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGC    >  GWNJ‑0478:712:GW2002102894th:2:2206:10106:25779/1‑150 (MQ=60)
                                                                                                                         ctCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaacaactactgataaaaaacaaaaat    >  GWNJ‑0478:712:GW2002102894th:2:2114:13755:61680/3‑47 (MQ=60)
                                                                                                                          CCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttc  >  GWNJ‑0478:712:GW2002102894th:2:2206:9962:89522/1‑90 (MQ=60)
                                                                                                                                             |                                                                                                                                  
GAATGCGCCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGCCT  >  AM260479/1239917‑1240187

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