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 155 25.8 2285122 88.0% 2010907 104.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,841,200 +G intergenic (‑2736/‑2680) h16_A1684 ← / → ycaC conserved hypothetical protein/YcaC related amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,841,1971.G97.2% 134.3 / ‑2.9 36intergenic (‑2733/‑2683)h16_A1684/ycaCconserved hypothetical protein/YcaC related amidohydrolase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (18/17);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCG  >  AM260479/1841067‑1841319
                                                                                                                                   |                                                                                                                           
gATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAggg                                                                                                                         >  2:1139337/1‑134 (MQ=255)
                  gcggcAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCgg                                                                                                         >  1:889067/1‑133 (MQ=255)
                    ggcAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGAt                                                                                                      >  1:703336/1‑134 (MQ=255)
                     gcAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGAtt                                                                                                     >  1:728113/1‑134 (MQ=255)
                       aaTGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTcc                                                                                                   <  2:728113/134‑1 (MQ=255)
                         tGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCcgc                                                                                                 <  2:889067/134‑1 (MQ=255)
                              ccAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACt                                                                                            <  2:439022/134‑1 (MQ=255)
                                                          gTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGAc                                                                <  1:347488/134‑1 (MQ=255)
                                                           tACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGAcc                                                               >  2:142926/1‑134 (MQ=255)
                                                           tACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGAcc                                                               >  1:792080/1‑134 (MQ=255)
                                                                    gcgcATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCg                                                      >  2:588489/1‑134 (MQ=255)
                                                                    gcgcATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCg                                                      >  2:882968/1‑134 (MQ=255)
                                                                         tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGAcc                                                               <  1:1037170/120‑1 (MQ=255)
                                                                         tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGAGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGAcc                                                               >  2:1037170/1‑120 (MQ=255)
                                                                               ccGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCg                                           <  1:1137011/134‑1 (MQ=255)
                                                                                cGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGc                                          <  1:670061/134‑1 (MQ=255)
                                                                                       ggCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCg                                    <  1:142926/133‑1 (MQ=255)
                                                                                             cgccTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGa                             <  2:792080/134‑1 (MQ=255)
                                                                                                cTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAgg                           >  1:1104617/1‑133 (MQ=255)
                                                                                                 tGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGccc                                     <  1:285901/122‑1 (MQ=255)
                                                                                                 tGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGccc                                     >  2:285901/1‑122 (MQ=255)
                                                                                                 tGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAgg                           >  1:1066809/1‑132 (MQ=255)
                                                                                                   cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGccc                                     >  1:350253/1‑120 (MQ=255)
                                                                                                   cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGccc                                     <  2:350253/120‑1 (MQ=255)
                                                                                                             tGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCcgc                                                                                                 >  2:1134108/1‑50 (MQ=255)
                                                                                                             tGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCcgc                                                                                                 <  1:1134108/50‑1 (MQ=255)
                                                                                                                 aGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATggg           >  1:405391/1‑132 (MQ=255)
                                                                                                                 aGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATggg           <  2:405391/132‑1 (MQ=255)
                                                                                                                   tCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGaa                                                                                              <  1:414186/47‑1 (MQ=39)
                                                                                                                   tCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGaa                                                                                              >  2:414186/1‑47 (MQ=39)
                                                                                                                        cAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCcgcg  >  2:1087115/1‑134 (MQ=255)
                                                                                                                         aGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                       >  2:985781/1‑64 (MQ=255)
                                                                                                                         aGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                       <  1:985781/64‑1 (MQ=255)
                                                                                                                          ggAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAg                                                                                     <  2:218499/49‑1 (MQ=38)
                                                                                                                          ggAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAg                                                                                     >  1:218499/1‑49 (MQ=38)
                                                                                                                            aaGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                       >  2:416902/1‑61 (MQ=255)
                                                                                                                            aaGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                       <  1:416902/61‑1 (MQ=255)
                                                                                                                                   |                                                                                                                           
GATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCG  >  AM260479/1841067‑1841319

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGTGGC  >  AM260479/1841070‑1841348
                                                                                                                                   |                                                                                                                                                    
CTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGG                                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2202:20989:18907/1‑149 (MQ=60)
  CGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATT                                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2109:4411:4556/1‑150 (MQ=60)
   GAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTC                                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2111:4665:41257/1‑150 (MQ=60)
                      tttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCAC                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2216:1472:57724/80‑1 (MQ=60)
                         ttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1208:2884:91224/77‑1 (MQ=60)
                                         CGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCC                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2112:6221:82410/1‑150 (MQ=60)
                                               gtgttcttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagaCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATC                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1211:8565:6560/70‑1 (MQ=60)
                                               acggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATC                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2207:14749:94477/96‑1 (MQ=60)
                                               aagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATC                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1209:7971:6322/91‑1 (MQ=60)
                                                      GGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTC                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1204:12965:7866/150‑1 (MQ=60)
                                                       GTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCC                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2210:14936:75734/150‑1 (MQ=60)
                                                        tcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaggttcagGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCC                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1202:13041:6031/76‑1 (MQ=60)
                                                                           GCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAG                                                         <  GWNJ‑0478:712:GW2002102894th:2:2216:8937:72137/150‑1 (MQ=60)
                                                                            CCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGG                                                        <  GWNJ‑0478:712:GW2002102894th:2:2106:3801:6863/150‑1 (MQ=60)
                                                                              GCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaggcagcgct                                                      >  GWNJ‑0478:712:GW2002102894th:2:1206:13372:73656/1‑63 (MQ=60)
                                                                              GCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGT                                                       >  GWNJ‑0478:712:GW2002102894th:2:2215:15124:18375/1‑149 (MQ=60)
                                                                               CTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGT                                                       >  GWNJ‑0478:712:GW2002102894th:2:2213:15830:30247/1‑148 (MQ=60)
                                                                                 GTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCC                                                   >  GWNJ‑0478:712:GW2002102894th:2:1204:14084:24734/1‑150 (MQ=60)
                                                                                 GTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCC                                                   <  GWNJ‑0478:712:GW2002102894th:2:1201:2695:23135/150‑1 (MQ=60)
                                                                                   AGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCC                                                  <  GWNJ‑0478:712:GW2002102894th:2:1109:12662:67234/149‑1 (MQ=60)
                                                                                               GCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCT                                     >  GWNJ‑0478:712:GW2002102894th:2:1208:5148:44341/1‑150 (MQ=60)
                                                                                                        CATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaa                             >  GWNJ‑0478:712:GW2002102894th:2:1112:9393:67313/1‑79 (MQ=60)
                                                                                                                           GTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCG           >  GWNJ‑0478:712:GW2002102894th:2:2116:13061:43583/1‑148 (MQ=60)
                                                                                                                                  GGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATctgtctcttatacacatctccgagcccacgagacactcgctaatc  >  GWNJ‑0478:712:GW2002102894th:2:2213:6824:91905/1‑105 (MQ=60)
                                                                                                                                  GGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGTGGC  <  GWNJ‑0478:712:GW2002102894th:2:2101:4776:96414/150‑1 (MQ=60)
                                                                                                                                   |                                                                                                                                                    
CTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGTGGC  >  AM260479/1841070‑1841348

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