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 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% 145.1 / NA 38intergenic (‑53/‑222)h16_A1136/grpEconserved hypothetical protein/molecular chaperone GrpE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (14/24);  total (14/24)

ATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTC  >  AM260479/1239936‑1240171
                                                                                                                      |                                                                                                                      
aTAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttat                                                                                                         <  1:333139/134‑1 (MQ=255)
   aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                       <  2:963375/133‑1 (MQ=255)
   aaGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTttatta                                                                                                       <  1:233651/133‑1 (MQ=255)
              gACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttttgtc                                                                                                                >  2:582962/1‑108 (MQ=255)
              gACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttttgtc                                                                                                                <  1:582962/112‑5 (MQ=255)
                    cAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGAtt                                                                                     <  2:305690/134‑1 (MQ=255)
                      aTGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGc                                                                                   <  1:146653/134‑1 (MQ=255)
                           gCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                              <  1:269700/134‑1 (MQ=255)
                           gCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                              <  2:503445/134‑1 (MQ=255)
                               caAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttt                                                                                                                    <  1:138277/91‑1 (MQ=255)
                               caAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCAtttt                                                                                                                    >  2:138277/1‑91 (MQ=255)
                                  ggCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                              <  2:600604/127‑1 (MQ=255)
                                  ggCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                              >  1:600604/1‑127 (MQ=255)
                                             aCGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                             <  1:472029/133‑1 (MQ=255)
                                                aCTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGa                                                         >  1:930043/1‑134 (MQ=255)
                                                 cTGGCATCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                        <  2:930043/134‑1 (MQ=255)
                                                 cTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACt                                                                                                     >  2:443587/1‑89 (MQ=255)
                                                 cTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACt                                                                                                     <  1:443587/89‑1 (MQ=255)
                                                     cAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                        <  1:905448/130‑1 (MQ=255)
                                                     cAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAt                                                        >  2:905448/1‑130 (MQ=255)
                                                       gCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCata                                                  >  1:226807/1‑134 (MQ=255)
                                                       gCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCata                                                  >  1:697737/1‑134 (MQ=255)
                                                         ttATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                             <  1:562244/121‑1 (MQ=255)
                                                         ttATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                             >  2:562244/1‑121 (MQ=255)
                                                           aTGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATg                                                <  2:602283/132‑1 (MQ=255)
                                                                       ccTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTaa                                  <  1:1000709/134‑1 (MQ=255)
                                                                              gAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCgtgtg                           <  2:697737/134‑1 (MQ=255)
                                                                                 gagGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGttt                        >  2:565682/1‑134 (MQ=255)
                                                                                        aaaGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTc                 <  2:803575/134‑1 (MQ=255)
                                                                                                  cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                              >  1:723756/1‑63 (MQ=255)
                                                                                                  cATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATgc                                                                              <  2:723756/63‑1 (MQ=255)
                                                                                                   aTCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgcg        >  1:796166/1‑132 (MQ=255)
                                                                                                   aTCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAgcg        <  2:796166/132‑1 (MQ=255)
                                                                                                        cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCt                                                                                             <  1:835034/42‑1 (MQ=38)
                                                                                                        cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCt                                                                                             >  2:835034/1‑42 (MQ=38)
                                                                                                        cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                             >  2:797981/1‑74 (MQ=255)
                                                                                                        cacaACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTGCCGCGTCACGCCAGCGCCActc  <  2:759052/133‑1 (MQ=255)
                                                                                                        cacaACCCCCGGCATTTTTGTCACTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCg                                                             <  1:797981/74‑1 (MQ=255)
                                                                                                          caACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCt                                        >  2:118418/1‑93 (MQ=255)
                                                                                                          caACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCt                                        <  1:118418/93‑1 (MQ=255)
                                                                                                           aaCCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTcc                      <  2:387323/110‑1 (MQ=255)
                                                                                                           aaCCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTcc                      >  1:387323/1‑110 (MQ=255)
                                                                                                                      |                                                                                                                      
ATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTC  >  AM260479/1239936‑1240171

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
GCCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGCCTT  >  AM260479/1239923‑1240188
                                                                                                                                       |                                                                                                                                   
gtgggctcggagatgtgtataagagacagCCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTAC                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1111:8935:10836/121‑1 (MQ=60)
    agagacagGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGG                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2104:15644:78631/142‑1 (MQ=60)
            GATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAG                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1207:7167:42549/150‑1 (MQ=60)
             tcgcttccttgtgttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAG                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1215:15198:27114/60‑1 (MQ=60)
               AAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTG                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1114:14562:76432/149‑1 (MQ=60)
                 cgtgggctcggagatgtgtataagagacagGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGAT                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1213:1666:99519/119‑1 (MQ=60)
                         tttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGC                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2205:17670:60578/72‑1 (MQ=60)
                                CCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGC                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:2107:3129:91085/1‑150 (MQ=60)
                                  AATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAA                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1111:16460:61094/150‑1 (MQ=60)
                                       GGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGG                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1203:3771:23252/150‑1 (MQ=60)
                                              AGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGAT                                                                         >  GWNJ‑0478:712:GW2002102894th:2:2210:10160:46704/1‑150 (MQ=60)
                                                     GTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATAT                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1114:9176:37346/1‑150 (MQ=60)
                                                     GTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATAT                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2113:17632:89357/1‑150 (MQ=60)
                                                        gacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagTCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCC                                                               <  GWNJ‑0478:712:GW2002102894th:2:1108:2321:48070/93‑1 (MQ=60)
                                                        agatgtagctccgtctcgtgggctcggagatgtgtataagagacagAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCACTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCC                                                               <  GWNJ‑0478:712:GW2002102894th:2:2202:20345:17296/104‑1 (MQ=60)
                                                         CACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCC                                                               <  GWNJ‑0478:712:GW2002102894th:2:1215:4962:58870/149‑1 (MQ=60)
                                                         CACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCC                                                               <  GWNJ‑0478:712:GW2002102894th:2:2105:14348:40376/149‑1 (MQ=60)
                                                             ACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCC                                                          <  GWNJ‑0478:712:GW2002102894th:2:2208:19234:100245/150‑1 (MQ=60)
                                                                             gggctcggagatgtgtataagagacagGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTT                                          <  GWNJ‑0478:712:GW2002102894th:2:1108:19296:86757/123‑1 (MQ=60)
                                                                                   ACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCG                                    <  GWNJ‑0478:712:GW2002102894th:2:2215:5049:49907/150‑1 (MQ=60)
                                                                                                GGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCtgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtc                       >  GWNJ‑0478:712:GW2002102894th:2:2115:7152:46608/1‑94 (MQ=60)
                                                                                                 GCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCC                       >  GWNJ‑0478:712:GW2002102894th:2:2202:20164:4984/1‑149 (MQ=60)
                                                                                                         ATGAACCATCTCCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCtgtctctt               >  GWNJ‑0478:712:GW2002102894th:2:1209:6319:80255/1‑141 (MQ=60)
                                                                                                                    CCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaa    >  GWNJ‑0478:712:GW2002102894th:2:1207:9242:22831/1‑71 (MQ=60)
                                                                                                                    CCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCtgtctcttatacacatctccgagcccacgagacggagctacatccc   >  GWNJ‑0478:712:GW2002102894th:2:2106:6713:30258/1‑104 (MQ=60)
                                                                                                                    CCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGCCT   >  GWNJ‑0478:712:GW2002102894th:2:1202:12321:82815/1‑150 (MQ=60)
                                                                                                                    CCACAACCCCCGGCATTTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGCCTT  >  GWNJ‑0478:712:GW2002102894th:2:1102:4821:47166/1‑150 (MQ=60)
                                                                                                                                       |                                                                                                                                   
GCCGGAAAGGCAGATAAAGACGGTCATGACCGCCAATGCGGCCACAAGGCAGAGTGTCACGACTGGCAGCTTATGCATGGCCAACCTTCCGGAAGAGGCACAAAGATGAACCATCTCCACAACCCCCGGCATTTTGTCATTTTTATTACTGGGAATCTTAGTGATTGCGATGCAGCAAAGCAATAAGGCGTCGATTGCATATGCCTGGCCTTTGTAACCGTGTGTTTCCGCGTCACGCCAGCGCCACTCTTGAAAACCTGCGCCTT  >  AM260479/1239923‑1240188

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