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 3,615,722 +A intergenic (+48/‑141) prfA → / → h16_A3341 protein chain release factor A (RF‑1)/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,615,7201.A100.0% 121.5 / NA 32intergenic (+46/‑143)prfA/h16_A3341protein chain release factor A (RF‑1)/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (13/19);  total (13/19)

CGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACGCGTAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACT  >  AM260479/3615593‑3615836
                                                                                                                                 |                                                                                                                    
cGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTaaagt                                                                                                                 >  1:514619/1‑131 (MQ=255)
  gCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTaaagtcg                                                                                                               <  1:256175/134‑6 (MQ=255)
  gCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTaaagtcg                                                                                                               <  2:646071/134‑6 (MQ=255)
  gCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTaaagtcg                                                                                                               <  1:563433/134‑6 (MQ=255)
      cATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGa                                                                                                            >  1:418222/1‑134 (MQ=255)
       aTGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGa                                                                                                            >  2:825457/1‑133 (MQ=255)
       aTGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGa                                                                                                            >  1:58435/1‑133 (MQ=255)
          gATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATc                                                                                                        <  1:997707/134‑1 (MQ=255)
                      tCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCtt                                                                                            <  1:825457/134‑1 (MQ=255)
                       cGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTg                                                                                           >  1:647466/1‑134 (MQ=255)
                        ggCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTaaagt                                                                                                                 >  2:836680/1‑107 (MQ=255)
                        ggCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTaaagt                                                                                                                 <  1:836680/110‑4 (MQ=255)
                            cTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAAt                                                                                      >  2:606995/1‑134 (MQ=255)
                                cTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCtt                                                                                  <  2:58435/134‑1 (MQ=255)
                                cTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCtt                                                                                  <  1:685404/134‑1 (MQ=255)
                                 ttCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTaaa                                                                                       >  1:585059/1‑128 (MQ=255)
                                 ttCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTaaa                                                                                       <  2:585059/128‑1 (MQ=255)
                                   cGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGa                                                                                >  1:530849/1‑133 (MQ=255)
                                     gAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAAc                                                                             <  1:606995/134‑1 (MQ=255)
                                     gAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAAc                                                                             <  1:639928/134‑1 (MQ=255)
                                          ccAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCAc                                                                        <  2:504888/134‑1 (MQ=255)
                                               cTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTgg                                                                     >  2:614547/1‑132 (MQ=255)
                                                    cAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAgggg                                                               >  1:22935/1‑133 (MQ=255)
                                                       cTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAg                                                           <  2:530849/134‑1 (MQ=255)
                                                       cTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAg                                                           <  2:418222/134‑1 (MQ=255)
                                                                   ggCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCt                                                 <  2:792483/132‑1 (MQ=255)
                                                                                       acaAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCaaa                           <  2:647466/134‑1 (MQ=255)
                                                                                          aaaTGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACgg                        <  2:514619/134‑1 (MQ=255)
                                                                                          aaaTGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACgg                        <  2:22935/134‑1 (MQ=255)
                                                                                                 gTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAg                                                           <  2:456843/92‑1 (MQ=255)
                                                                                                 gTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAg                                                           >  1:456843/1‑92 (MQ=255)
                                                                                                    cAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGcc                                                                                                              <  1:675473/38‑1 (MQ=255)
                                                                                                    cAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGcc                                                                                                              >  2:675473/1‑38 (MQ=255)
                                                                                                      aCCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTg                                                                                                    <  2:380074/46‑1 (MQ=255)
                                                                                                      aCCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTg                                                                                                    >  1:380074/1‑46 (MQ=255)
                                                                                                           ccTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAg       >  1:872995/1‑134 (MQ=255)
                                                                                                                cGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCAc   <  2:102119/133‑1 (MQ=255)
                                                                                                                cGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACt  <  1:614547/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                    
CGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACGCGTAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACT  >  AM260479/3615593‑3615836

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATCATGGACGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACGCGTAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACTCCCGAGAACATACCTAGGACGAAATCATG  >  AM260479/3615584‑3615865
                                                                                                                                            |                                                                                                                                               
GATCATGGACGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGAT                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1210:3747:57573/1‑150 (MQ=60)
 ATCATGGACGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGAT                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1102:13125:100187/1‑149 (MQ=60)
         aagagacagGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGT                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2205:2713:71546/140‑1 (MQ=60)
          GGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTAT                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1201:16989:26463/150‑1 (MQ=60)
          GGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTAT                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2105:11022:48261/150‑1 (MQ=60)
             ttttgttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtatgagagacagCCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTT                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2110:3988:58260/68‑1 (MQ=60)
                 TGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGT                                                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:2109:21131:92934/1‑150 (MQ=60)
                  GGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTA                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2215:10150:33416/150‑1 (MQ=60)
                     ccatgtcgccgtccatgatcatgtcgatcttgtacagcgtcaggttgaTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAAT                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:1106:12056:84297/49‑150 (MQ=60)
                           GCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTG                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:2104:19944:91567/1‑149 (MQ=60)
                             TGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAA                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2103:18920:18689/1‑149 (MQ=60)
                              GTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAAC                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2204:15926:97836/150‑1 (MQ=60)
                                        GCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTA                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2113:10769:17202/150‑1 (MQ=60)
                                             GGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGT                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2106:18272:32350/150‑1 (MQ=60)
                                             GGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGT                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2109:19350:46833/150‑1 (MQ=60)
                                              GAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGT                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1103:7463:59953/1‑149 (MQ=60)
                                                   tttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagcccgaacaccgcgtaAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGC                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1202:10390:41545/61‑1 (MQ=60)
                                                                      GCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttc                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1203:1511:85459/1‑90 (MQ=60)
                                                                                           GCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCtgtctcttatacacatctccgagccc                                              >  GWNJ‑0478:712:GW2002102894th:2:1212:4416:75958/1‑123 (MQ=60)
                                                                                                ACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaa                                         >  GWNJ‑0478:712:GW2002102894th:2:1209:9152:37139/1‑78 (MQ=60)
                                                                                                     ATGTCGTGCAACCGGCCTTCCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGT                                   >  GWNJ‑0478:712:GW2002102894th:2:2206:15529:98549/1‑150 (MQ=60)
                                                                                                                        CCGCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACTCCCGAGAACATACCT                <  GWNJ‑0478:712:GW2002102894th:2:2106:13796:84767/150‑1 (MQ=60)
                                                                                                                          GCCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaa               >  GWNJ‑0478:712:GW2002102894th:2:1213:18832:14075/1‑73 (MQ=60)
                                                                                                                           CCCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACTCCCGAGAACATACCTAGG             >  GWNJ‑0478:712:GW2002102894th:2:2214:13598:46890/1‑150 (MQ=60)
                                                                                                                            CCCGAACACCGCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaagaagtcacgaagagg            >  GWNJ‑0478:712:GW2002102894th:2:1114:9807:53368/1‑56 (MQ=59)
                                                                                                                                      GCGTAAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACTCCCGAGAACATACCTAGGACGAAATCATG  >  GWNJ‑0478:712:GW2002102894th:2:2204:18983:2831/1‑150 (MQ=60)
                                                                                                                                            |                                                                                                                                               
GATCATGGACGGCGACATGGATGAGCTGCTGTCGGCACTGGCTTCGGAACACCAGGCTGATCAGCTTGCTGCACTTGGCGAAGAGTCGTAAGCACCACAAAATGTCGTGCAACCGGCCTTCCGCCCCGAACACGCGTAAGTCGCCGATATCTTTGTATTTCTTGTAAATGCTTGAAACAGCACTGGTAGGGGTAAGGCGGCTGTCTATAATCGATTTCGACCGGCAAACGGAGAATGGTGCCGGGTAGTCACTCCCGAGAACATACCTAGGACGAAATCATG  >  AM260479/3615584‑3615865

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: