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 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 782,333 +A intergenic (‑683/‑39) h16_A0723 ← / → h16_A0724 hypothetical membrane protein/hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479782,3331.A100.0% 124.0 / NA 32intergenic (‑683/‑39)h16_A0723/h16_A0724hypothetical membrane protein/hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (17/15);  total (17/15)

GTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCATGGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATC  >  AM260479/782209‑782449
                                                                                                                              |                                                                                                                    
gTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGa                                                                                                               >  2:386293/1‑134 (MQ=255)
gTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGa                                                                                                               >  2:171027/1‑134 (MQ=255)
        tGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGcc                                                                                                       >  1:927578/1‑134 (MQ=255)
         gTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGc                                                                                                        <  1:822021/132‑1 (MQ=255)
               cGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGgcag                                                                                                 >  1:238236/1‑133 (MQ=255)
               cGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGgcag                                                                                                 >  1:264361/1‑133 (MQ=255)
                 gTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              <  2:986560/134‑1 (MQ=255)
                 gTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              <  1:171027/134‑1 (MQ=255)
                 gTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              <  2:264361/134‑1 (MQ=255)
                                gAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGa                                                                               >  2:496795/1‑134 (MQ=255)
                                            tACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCAtt                                                                   >  1:410987/1‑134 (MQ=255)
                                               aTCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGAt                                                                              <  1:386946/120‑1 (MQ=255)
                                               aTCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGAt                                                                              >  2:386946/1‑120 (MQ=255)
                                               aTCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGc                                                                >  2:894495/1‑134 (MQ=255)
                                               aTCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGc                                                                >  2:699358/1‑134 (MQ=255)
                                                                          gaagaaATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCAGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTa                                       >  1:983846/1‑132 (MQ=255)
                                                                             gaaATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCg                                  >  2:518281/1‑134 (MQ=255)
                                                                               aaTCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              <  1:30617/72‑1 (MQ=255)
                                                                               aaTCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              >  2:30617/1‑72 (MQ=255)
                                                                                   aGGCCTGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCg                            >  2:961325/1‑134 (MQ=255)
                                                                                    ggCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGt                           <  1:518281/134‑1 (MQ=255)
                                                                                    ggCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGt                           <  1:699358/134‑1 (MQ=255)
                                                                                    ggCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCAGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGt                           <  2:983846/134‑1 (MQ=255)
                                                                                       cGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCt                        <  1:386293/134‑1 (MQ=255)
                                                                                                   aGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATgg            <  1:961325/134‑1 (MQ=255)
                                                                                                    gAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGt           <  2:927578/134‑1 (MQ=255)
                                                                                                    gAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGt           <  1:496795/134‑1 (MQ=255)
                                                                                                             ggggATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCt       >  2:796935/1‑129 (MQ=255)
                                                                                                             ggggATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCt       <  1:796935/129‑1 (MQ=255)
                                                                                                             ggggATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATc  >  1:966425/1‑134 (MQ=255)
                                                                                                                gATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              <  2:919539/39‑1 (MQ=255)
                                                                                                                gATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCaa                                                                                              >  1:919539/1‑39 (MQ=255)
                                                                                                                              |                                                                                                                    
GTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCATGGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATC  >  AM260479/782209‑782449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
CCCCGCGCAGGCCTTGGTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCATGGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCAAAACG  >  AM260479/782193‑782470
                                                                                                                                             |                                                                                                                                          
CCCCGCGCAGGCCTTGGTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGA                                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2208:14743:46071/150‑1 (MQ=60)
  CCGCGCAGGCCTTGGTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACC                                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2210:3182:88555/150‑1 (MQ=60)
       agGGCCTTGGTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGC                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2215:17163:52039/148‑1 (MQ=60)
        agctccgtctcgtgggctcggagatgtgtataagagacagGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCC                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1207:12250:101100/110‑1 (MQ=60)
         GGCCTTGGTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCC                                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2211:16905:100179/1‑150 (MQ=60)
                GTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGC                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1114:10522:95274/1‑149 (MQ=60)
                GTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGC                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1201:18017:64271/1‑149 (MQ=60)
                        TGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCG                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2204:13928:39033/148‑1 (MQ=60)
                                cgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGA                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1105:11083:15569/102‑1 (MQ=60)
                                GGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGA                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1113:17524:10287/150‑1 (MQ=60)
                                             CCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTG                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1210:10188:4547/1‑150 (MQ=60)
                                                GAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCT                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1209:9656:67019/150‑1 (MQ=60)
                                                                           CCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCAGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTAC                                                           >  GWNJ‑0478:712:GW2002102894th:2:2214:7435:22014/1‑148 (MQ=60)
                                                                                      cgtgggctcggagatgtgtataagagacagGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGC                                              <  GWNJ‑0478:712:GW2002102894th:2:1108:1628:8350/120‑1 (MQ=60)
                                                                                        cgtgggctcggagatgtgtataagagacagATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGC                                              <  GWNJ‑0478:712:GW2002102894th:2:1116:7881:52579/118‑1 (MQ=60)
                                                                                                cggagatgtgtataagagacagATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTA                                    <  GWNJ‑0478:712:GW2002102894th:2:2102:10559:18680/128‑1 (MQ=60)
                                                                                                   AGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGG                                 <  GWNJ‑0478:712:GW2002102894th:2:2113:6373:76301/150‑1 (MQ=60)
                                                                                                   AGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGG                                 <  GWNJ‑0478:712:GW2002102894th:2:1216:12607:96834/150‑1 (MQ=60)
                                                                                                      CCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGC                              <  GWNJ‑0478:712:GW2002102894th:2:1209:6058:63412/150‑1 (MQ=60)
                                                                                                              agagacagATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCG                      <  GWNJ‑0478:712:GW2002102894th:2:2209:3120:74816/142‑1 (MQ=60)
                                                                                                              GCCACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCG                      >  GWNJ‑0478:712:GW2002102894th:2:2213:19150:26642/1‑150 (MQ=60)
                                                                                                                 ACAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaa                    >  GWNJ‑0478:712:GW2002102894th:2:2211:15533:54159/1‑69 (MQ=60)
                                                                                                                  CAGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGC                  <  GWNJ‑0478:712:GW2002102894th:2:2212:5689:96582/150‑1 (MQ=60)
                                                                                                                   AGAATATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCA                 <  GWNJ‑0478:712:GW2002102894th:2:1214:7890:75702/150‑1 (MQ=60)
                                                                                                                       gacagGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGG             <  GWNJ‑0478:712:GW2002102894th:2:2101:21139:32081/145‑1 (MQ=60)
                                                                                                                        ATCAGGGGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGG             <  GWNJ‑0478:712:GW2002102894th:2:2203:14088:66526/149‑1 (MQ=60)
                                                                                                                              GGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCA       >  GWNJ‑0478:712:GW2002102894th:2:2209:17179:13824/1‑149 (MQ=60)
                                                                                                                              GGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCA       >  GWNJ‑0478:712:GW2002102894th:2:2215:16004:95761/1‑149 (MQ=60)
                                                                                                                              GGGATGGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCA       >  GWNJ‑0478:712:GW2002102894th:2:2207:14008:94125/1‑149 (MQ=60)
                                                                                                                                   GGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCAAAACG  >  GWNJ‑0478:712:GW2002102894th:2:2210:13164:91901/1‑149 (MQ=60)
                                                                                                                                   GGGGAGCGATGAGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCAAAACG  >  GWNJ‑0478:712:GW2002102894th:2:2209:6301:31812/1‑149 (MQ=60)
                                                                                                                                             |                                                                                                                                          
CCCCGCGCAGGCCTTGGTTCTGCTTGTACATCGGTGGGCGCCTGACCAGAATGGAGCCAGTACATCGATGCATCGCCGCTGCATTGCCAAGAAGAAATCAGGCCGGCGCTGCCACAGAATATCAGGGGGATGGGGAGCATGGCCTGGACCATAGCCCGGCAGCAACGCCGGGGTCAGTGATGTCGCGGCATTGGCTGGACCTGGAGCGTCTCCGGGCCTACAGCGTTGCCGTGCTGATCATCTATGGTGCCTTGATCGTGGCATGGGCATTCAAAACG  >  AM260479/782193‑782470

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