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
A2 F15 I1 R1 9 36.2 1317874 96.5% 1271748 132.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 758,061 G→T G125* (GGA→TGA)  sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,0610GT100.0% 168.4 / NA 48G125* (GGA→TGA) sdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (25/23);  total (25/23)

CCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGG  >  NC_000913/757918‑758200
                                                                                                                                               |                                                                                                                                           
ccccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACa                                                                                                                                         >  5:2840/1‑148 (MQ=255)
 cccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACaa                                                                                                                                        >  1:58187/1‑148 (MQ=255)
 cccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACaa                                                                                                                                        >  2:157044/1‑148 (MQ=255)
  ccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGa                                                                                                                                           <  1:153073/144‑1 (MQ=255)
  ccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGa                                                                                                                                           >  2:153073/1‑144 (MQ=255)
  ccGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACaaaa                                                                                                                                      <  2:58187/149‑1 (MQ=255)
          ggCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCg                                                                                                                                  >  1:5165/1‑145 (MQ=255)
          ggCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCg                                                                                                                                  <  2:5165/145‑1 (MQ=255)
                          ggCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACa                                                                                                              <  1:110277/149‑1 (MQ=255)
                               gaagaTTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGc                                                                                                         <  6:2840/149‑1 (MQ=255)
                                    ttGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCagag                                                                                                    <  1:157044/149‑1 (MQ=255)
                                          ttCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAgcgc                                                                                              <  2:151817/149‑1 (MQ=255)
                                                     ccAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTAc                                                                                                               <  6:56021/121‑1 (MQ=255)
                                                     ccAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTAc                                                                                                               >  5:56021/1‑121 (MQ=255)
                                                         gTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACggg                                                                               >  1:169631/1‑149 (MQ=255)
                                                                ggTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCa                                                                                                       <  7:34276/118‑1 (MQ=255)
                                                                ggTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCa                                                                                                       >  8:34276/1‑118 (MQ=255)
                                                                      ccGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTa                                                                  >  6:168685/1‑149 (MQ=255)
                                                                       cgcgATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAg                                                                                                             >  3:92400/1‑105 (MQ=255)
                                                                       cgcgATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAg                                                                                                             <  4:92400/105‑1 (MQ=255)
                                                                        gcgATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTAtt                                                                <  2:169631/149‑1 (MQ=255)
                                                                          gATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCt                                                                             >  3:118722/1‑134 (MQ=255)
                                                                          gATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCt                                                                             <  4:118722/134‑1 (MQ=255)
                                                                                   gtAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGAt                                                                                                           <  2:63687/95‑1 (MQ=255)
                                                                                   gtAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGAt                                                                                                           >  1:63687/1‑95 (MQ=255)
                                                                                         aTGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATgttgtt                                                >  3:62129/1‑148 (MQ=255)
                                                                                           gggACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCaa                                             >  2:119694/1‑149 (MQ=255)
                                                                                           gggACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCACCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAtgt                                                    >  2:120383/1‑142 (MQ=255)
                                                                                           gggACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCACCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCAtgt                                                    <  1:120383/142‑1 (MQ=255)
                                                                                            ggACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGaaaaa                                                                                        <  7:114370/105‑1 (MQ=255)
                                                                                            ggACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGaaaaa                                                                                        >  8:114370/1‑105 (MQ=255)
                                                                                            ggACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAAc                                            >  3:138939/1‑149 (MQ=255)
                                                                                                aaTTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTcgc                                                                                                                         <  5:132981/68‑1 (MQ=255)
                                                                                                aaTTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTcgc                                                                                                                         >  6:132981/1‑68 (MQ=255)
                                                                                                 aTTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCtt                                       >  1:86828/1‑149 (MQ=255)
                                                                                                         cgcAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCt                               >  2:52386/1‑149 (MQ=255)
                                                                                                             atatGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCt                           <  1:52386/149‑1 (MQ=255)
                                                                                                                tGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggt                        >  1:43959/1‑149 (MQ=255)
                                                                                                                 gagaAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTctct                                                            <  4:34160/112‑1 (MQ=255)
                                                                                                                 gagaAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTctct                                                            >  3:34160/1‑112 (MQ=255)
                                                                                                                  agaAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggtg                       <  7:73438/148‑1 (MQ=255)
                                                                                                                  agaAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCtggtg                       >  8:73438/1‑148 (MQ=255)
                                                                                                                                 aCCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGt         <  1:119694/147‑1 (MQ=255)
                                                                                                                                      ttGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAAc                                                                                       >  6:526/1‑64 (MQ=255)
                                                                                                                                      ttGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAAc                                                                                       <  5:526/64‑1 (MQ=255)
                                                                                                                                      ttGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCgg  <  4:62129/149‑1 (MQ=255)
                                                                                                                                        gAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTAt                                                                         <  6:147276/76‑1 (MQ=255)
                                                                                                                                        gAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTAt                                                                         >  5:147276/1‑76 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGG  >  NC_000913/757918‑758200

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
ACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTT  >  NC_000913/757917‑758194
                                                                                                                                                |                                                                                                                                     
ACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAA                                                                                                                                  >  NB501016:228:HFK2FAFXY:3:11402:7622:3644/1‑150 (MQ=60)
 CCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAA                                                                                                                                 >  NB501016:228:HFK2FAFXY:1:11210:6799:8543/1‑150 (MQ=60)
  CCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGAC                                                                                                                                    <  NB501016:228:HFK2FAFXY:1:21210:21791:7171/146‑1 (MQ=60)
          CGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGC                                                                                                                           >  NB501016:228:HFK2FAFXY:1:11102:16111:17338/1‑147 (MQ=60)
                          GGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAG                                                                                                       <  NB501016:228:HFK2FAFXY:1:21105:10819:17072/151‑1 (MQ=60)
                                    ATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGC                                                                                             <  NB501016:228:HFK2FAFXY:1:21211:21225:17670/151‑1 (MQ=60)
                                                     GCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACA                                                                                                        >  NB501016:228:HFK2FAFXY:3:11512:8278:1171/1‑123 (MQ=60)
                                                         GGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGC                                                                        >  NB501016:228:HFK2FAFXY:1:21304:10057:14584/1‑151 (MQ=60)
                                                                CGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAG                                                                                                <  NB501016:228:HFK2FAFXY:4:11504:4002:13561/120‑1 (MQ=60)
                                                                       CCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGA                                                                                                      >  NB501016:228:HFK2FAFXY:2:21109:22802:4445/1‑107 (MQ=60)
                                                                          CGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTG                                                                      >  NB501016:228:HFK2FAFXY:2:21210:15679:11059/1‑136 (MQ=60)
                                                                                   GGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATG                                                                                                    >  NB501016:228:HFK2FAFXY:1:11212:15371:9663/1‑97 (MQ=60)
                                                                                    GTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTT                                             >  NB501016:228:HFK2FAFXY:3:21607:13721:15527/1‑151 (MQ=60)
                                                                                         CATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTC                                         >  NB501016:228:HFK2FAFXY:2:11306:22554:18280/1‑150 (MQ=60)
                                                                                           TGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCACCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTT                                             <  NB501016:228:HFK2FAFXY:1:21109:18439:13339/144‑1 (MQ=60)
                                                                                            GGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACC                                     >  NB501016:228:HFK2FAFXY:2:21308:13228:2789/1‑151 (MQ=60)
                                                                                            GGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAAC                                                                                 <  NB501016:228:HFK2FAFXY:4:21508:3771:7815/107‑1 (MQ=60)
                                                                                                CAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCG                                                                                                                  <  NB501016:228:HFK2FAFXY:3:21508:19626:13616/70‑1 (MQ=60)
                                                                                                 AATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTG                                >  NB501016:228:HFK2FAFXY:1:11309:20028:3181/1‑151 (MQ=60)
                                                                                                             AATATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTG                    <  NB501016:228:HFK2FAFXY:1:11208:19611:4350/151‑1 (MQ=60)
                                                                                                                ATGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTG                 >  NB501016:228:HFK2FAFXY:1:11205:7064:2161/1‑151 (MQ=60)
                                                                                                                 TGAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGCCGGGCTGTATGAATGTATTCTCTG                                                     >  NB501016:228:HFK2FAFXY:2:11205:17005:18887/1‑114 (MQ=60)
                                                                                                                  GAGAAAATTAAGCCTTACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGG                <  NB501016:228:HFK2FAFXY:4:11612:19531:10303/150‑1 (MQ=60)
                                                                                                                                 TACCTGTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTT  <  NB501016:228:HFK2FAFXY:1:21109:3892:7922/149‑1 (MQ=60)
                                                                                                                                      GTTGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACT                                                                                <  NB501016:228:HFK2FAFXY:3:11401:24307:5109/66‑1 (MQ=60)
                                                                                                                                        TGAATAATTGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATG                                                                  >  NB501016:228:HFK2FAFXY:3:21602:17384:18245/1‑78 (MQ=60)
                                                                                                                                                |                                                                                                                                     
ACCCCGATTTCGGCACTCAACCAGCCGGGCAAGAAGATTGTGATTCGCCCGCTGCCAGGTTTACCGGTGATCCGCGATTTGGTGGTAGACATGGGACAATTCTATGCGCAATATGAGAAAATTAAGCCTTACCTGTTGAATAATGGACAAAATCCGCCAGCTCGCGAGCATTTACAGATGCCAGAGCAGCGCGAAAAACTCGACGGGCTGTATGAATGTATTCTCTGCGCATGTTGTTCAACCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTT  >  NC_000913/757917‑758194

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