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
A4 F4 I1 R1 17 37.4 1246952 98.2% 1224506 142.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,176,112 G→A P241L (CCC→CTC)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,1120GA100.0% 137.6 / NA 42P241L (CCC→CTC) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

CTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGG  >  NC_000913/3175972‑3176254
                                                                                                                                            |                                                                                                                                              
cTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATg                                                                                                                                        <  5:74403/149‑1 (MQ=255)
 tGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATgg                                                                                                                                       >  1:29223/1‑149 (MQ=255)
              taaataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTaaa                                                                                                                          <  3:3650/149‑1 (MQ=255)
                aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGt                                                                                                                        >  1:176155/1‑149 (MQ=255)
                aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGt                                                                                                                        >  2:128257/1‑149 (MQ=255)
                 ataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGtt                                                                                                                       >  2:178448/1‑149 (MQ=255)
                   aaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGAt                                                                                                                                         <  3:123861/129‑1 (MQ=255)
                   aaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGAt                                                                                                                                         >  4:123861/1‑129 (MQ=255)
                     agaagCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGaa                                                                                                                   >  6:18393/1‑149 (MQ=255)
                           gTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGc                                                                                                             >  8:32900/1‑149 (MQ=255)
                                    aGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTaaa                                                                                                     <  7:101719/148‑1 (MQ=255)
                                    aGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTaaa                                                                                                     <  1:128257/148‑1 (MQ=255)
                                     gTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACt                                                                                                   >  7:41469/1‑149 (MQ=255)
                                        gCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGca                                                                                                <  2:29223/149‑1 (MQ=255)
                                           ttCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGcaca                                                                                              >  7:60278/1‑148 (MQ=255)
                                              tGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAgg                                                                                          >  8:110811/1‑149 (MQ=255)
                                                caGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTc                                                                                        >  7:33016/2‑149 (MQ=255)
                                                                 aaaCGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCCGCCGAGCGCGATGGTATCCAGCGTAAAGTGGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCcagcag                                                                       <  3:19599/149‑1 (MQ=255)
                                                                               cAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTc                                                          <  1:95439/148‑1 (MQ=255)
                                                                               cAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTc                                                          >  2:95439/1‑148 (MQ=255)
                                                                                     aTGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCga                                                   <  7:32900/149‑1 (MQ=255)
                                                                                     aTGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCga                                                   <  8:41469/149‑1 (MQ=255)
                                                                                        cACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTc                                                     >  4:24550/1‑144 (MQ=255)
                                                                                          ccTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATcc                                              <  3:24550/149‑1 (MQ=255)
                                                                                              ggtggtCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                                                                                    <  7:132353/91‑1 (MQ=255)
                                                                                              ggtggtCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                                                                                    >  8:132353/1‑91 (MQ=255)
                                                                                                tggtCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTcc                                        <  2:176155/149‑1 (MQ=255)
                                                                                                     aGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCtgat                                    <  1:178448/148‑1 (MQ=255)
                                                                                                             aTCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGAc                           >  7:8690/1‑149 (MQ=255)
                                                                                                             aTCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGAc                           <  5:18393/149‑1 (MQ=255)
                                                                                                             aTCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGAc                           <  8:60278/149‑1 (MQ=255)
                                                                                                                    tGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCtgatga                                  <  2:61346/135‑1 (MQ=255)
                                                                                                                    tGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCtgatga                                  >  1:61346/1‑135 (MQ=255)
                                                                                                                        aCTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCaa                 <  8:8690/148‑1 (MQ=255)
                                                                                                                          tCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTa              >  7:85940/1‑149 (MQ=255)
                                                                                                                          tCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTa              >  3:48366/1‑149 (MQ=255)
                                                                                                                            gagaGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAtt            >  2:64927/1‑149 (MQ=255)
                                                                                                                            gagaGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAt             <  8:85940/148‑1 (MQ=255)
                                                                                                                             agagTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAttt           <  4:48366/149‑1 (MQ=255)
                                                                                                                              gagTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAttt           >  7:2678/1‑148 (MQ=255)
                                                                                                                                    gCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGc    <  8:28164/149‑1 (MQ=255)
                                                                                                                                       aGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCgg  >  3:1973/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGG  >  NC_000913/3175972‑3176254

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGA  >  NC_000913/3175963‑3176255
                                                                                                                                                     |                                                                                                                                               
CGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAG                                                                                                                                                <  NB501016:228:HFK2FAFXY:3:21405:5245:5729/151‑1 (MQ=60)
        GCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGG                                                                                                                                        <  NB501016:228:HFK2FAFXY:3:11609:20328:10701/151‑1 (MQ=60)
         CTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGT                                                                                                                                       >  NB501016:228:HFK2FAFXY:1:11112:8696:7822/1‑151 (MQ=60)
                      GTAAATAAAGACGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAG                                                                                                                          <  NB501016:228:HFK2FAFXY:2:11102:5593:17355/151‑1 (MQ=60)
                        AAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTT                                                                                                                        >  NB501016:228:HFK2FAFXY:1:21311:15275:10493/1‑151 (MQ=60)
                           TAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATG                                                                                                                                         <  NB501016:228:HFK2FAFXY:2:21303:8756:20070/131‑1 (MQ=60)
                                            CAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAC                                                                                                     <  NB501016:228:HFK2FAFXY:1:21204:9451:2618/150‑1 (MQ=60)
                                            CAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAC                                                                                                     <  NB501016:228:HFK2FAFXY:4:21504:13235:3989/150‑1 (MQ=60)
                                             AGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTG                                                                                                   >  NB501016:228:HFK2FAFXY:4:11509:6917:6418/1‑151 (MQ=60)
                                                   CTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACAC                                                                                              >  NB501016:228:HFK2FAFXY:4:11607:10547:5953/1‑150 (MQ=60)
                                                        GCAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCG                                                                                        >  NB501016:228:HFK2FAFXY:4:11504:5629:14791/1‑151 (MQ=60)
                                                                         GAAACGGGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCCGCCGAGCGCGATGGTATCCAGCGTAAAGTGGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGG                                                                       <  NB501016:228:HFK2FAFXY:2:11111:25019:7727/151‑1 (MQ=60)
                                                                                      CCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCC                                                          <  NB501016:228:HFK2FAFXY:1:21102:13033:17962/151‑1 (MQ=60)
                                                                                             GATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAG                                                   <  NB501016:228:HFK2FAFXY:4:11504:11225:13548/151‑1 (MQ=60)
                                                                                             GATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAG                                                   <  NB501016:228:HFK2FAFXY:2:11208:15183:3283/151‑1 (MQ=60)
                                                                                                  ACCTCGGGGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCA                                              <  NB501016:228:HFK2FAFXY:2:11201:5707:19174/151‑1 (MQ=60)
                                                                                                      CGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACT                                                                                                    <  NB501016:228:HFK2FAFXY:4:21608:20958:19335/93‑1 (MQ=60)
                                                                                                             CAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATG                                    <  NB501016:228:HFK2FAFXY:1:21312:15134:8711/150‑1 (MQ=60)
                                                                                                                     CATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACC                           <  NB501016:228:HFK2FAFXY:3:11409:6954:9118/151‑1 (MQ=60)
                                                                                                                     CATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACC                           >  NB501016:228:HFK2FAFXY:4:11405:13070:7505/1‑151 (MQ=60)
                                                                                                                            ATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAA                                  >  NB501016:228:HFK2FAFXY:1:11301:23478:5119/1‑137 (MQ=60)
                                                                                                                                  CTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAT              >  NB501016:228:HFK2FAFXY:2:11302:22457:3288/1‑151 (MQ=60)
                                                                                                                                  CTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAT              >  NB501016:228:HFK2FAFXY:4:21408:11930:11052/1‑151 (MQ=60)
                                                                                                                                      AGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTG           >  NB501016:228:HFK2FAFXY:4:11402:2505:7705/1‑150 (MQ=60)
                                                                                                                                               CAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGA  >  NB501016:228:HFK2FAFXY:2:11102:21707:1063/1‑150 (MQ=60)
                                                                                                                                                   CGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCG                          <  NB501016:228:HFK2FAFXY:3:11410:21962:17866/122‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                               
CGCGGCGAGCTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGA  >  NC_000913/3175963‑3176255

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