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 F27 I1 R1 56 84.2 3463790 98.6% 3415296 146.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,767,984 T→G intergenic (+91/‑91) fpvA → / → tpaK TonB‑dependent outer membrane ferripyoverdine receptor FpvA/probable terephthalate transporter, MFS superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,767,9840TG100.0% 145.0 / NA 42intergenic (+91/‑91)fpvA/tpaKTonB‑dependent outer membrane ferripyoverdine receptor FpvA/probable terephthalate transporter, MFS superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (26/16);  total (26/16)

TGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG  >  NC_002947_CJ‑RC/4767844‑4768124
                                                                                                                                            |                                                                                                                                            
cgaaTGGTGACCCACGCAACTACCCCATCACCGCCAAGTCCAGCTTCTCATCGCCTGACGTTGAACGCAACAAACCCCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTa                                                                                                                                       <  4:403349/145‑1 (MQ=37)
   aTGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                    <  6:403511/148‑1 (MQ=255)
       tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatcatc                                                                                                                                 <  5:358229/147‑1 (MQ=255)
        gACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCgg                                                                                                                               >  6:51761/1‑148 (MQ=255)
             aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGt                                                                                                                          >  6:97739/1‑148 (MQ=255)
             aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGt                                                                                                                          >  2:131594/1‑148 (MQ=255)
             aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGt                                                                                                                          >  8:184136/1‑148 (MQ=255)
             aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGt                                                                                                                          >  4:274780/1‑148 (MQ=255)
                       aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGa                                                                                                                 <  4:139880/147‑1 (MQ=255)
                          aTCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAg                                                                                                             <  1:181264/148‑1 (MQ=255)
                            cACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAc                                                                                                           >  8:26698/1‑148 (MQ=255)
                             aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAcc                                                                                                          >  1:44296/1‑148 (MQ=255)
                                               tGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATa                                                                                        >  5:244036/1‑148 (MQ=255)
                                                          cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                             <  4:66172/148‑1 (MQ=255)
                                                          cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                             <  5:315902/148‑1 (MQ=255)
                                                               aaagCAAAAAACCGCAACCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGAtt                                                                         <  1:269320/144‑1 (MQ=17)
                                                                 cGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACATTTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAc                                                                      >  1:140366/1‑148 (MQ=255)
                                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTCAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGCATAACAATTCTAAAGATTAAc                                                                      >  8:147716/1‑147 (MQ=17)
                                                                   cAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGATTTGTGAGCGGATAACAATTCTTAAGATTAACt                                                                     >  4:253551/1‑147 (MQ=255)
                                                                    aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAGTc                                                                    >  7:303118/1‑147 (MQ=255)
                                                                     aaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGCCCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcac                                                                  >  8:357036/1‑148 (MQ=255)
                                                                      aaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcaca                                                                 <  3:283308/148‑1 (MQ=255)
                                                                          ccGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAgg                                                             <  6:221000/148‑1 (MQ=255)
                                                                           cGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcacaca                                                               >  5:229575/1‑145 (MQ=255)
                                                                            gCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcacaca                                                               >  8:56162/1‑144 (MQ=255)
                                                                              aCCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGat                                                         <  5:97739/148‑1 (MQ=255)
                                                                                                      aCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgc                                 >  6:375648/1‑148 (MQ=255)
                                                                                                      aCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgc                                 >  7:181468/1‑148 (MQ=255)
                                                                                                        tagtagTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgcgc                               >  4:147151/1‑148 (MQ=255)
                                                                                                           tagtCTCTATGGATGTCAGGTATGATTACTATTGAGAATTAATAATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGTTAACAATTCTTAAGATTAACTCACACAGGAGATTTCATATGTCCTTAGCACCATCGCGCGt                             <  6:353219/147‑1 (MQ=11)
                                                                                                           tagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGTCCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTg                            <  7:357036/148‑1 (MQ=255)
                                                                                                           tagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGt                             <  2:44296/147‑1 (MQ=255)
                                                                                                           tagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGt                             <  3:233589/147‑1 (MQ=255)
                                                                                                              tctcTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGaca                         >  5:248228/1‑148 (MQ=255)
                                                                                                               ctctATGGAGGTCAGGTATGGTTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGacac                        >  5:259480/1‑148 (MQ=255)
                                                                                                               ctctATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGacac                        >  7:174856/1‑148 (MQ=255)
                                                                                                               ctctATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGacac                        <  5:67747/148‑1 (MQ=255)
                                                                                                                    tGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTAAGCACCATCGCGCGTGACACTGcc                    >  4:331231/1‑147 (MQ=255)
                                                                                                                              gTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCg         >  6:248118/1‑148 (MQ=255)
                                                                                                                                   aTTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcg    >  6:14738/1‑148 (MQ=255)
                                                                                                                                   aTTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcg    >  6:20533/1‑148 (MQ=255)
                                                                                                                                     tACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAAGTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgcg  >  8:182305/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG  >  NC_002947_CJ‑RC/4767844‑4768124

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
GGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG  >  NC_002947_CJ‑RC/4767849‑4768104
                                                                                                                                       |                                                                                                                        
GGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCG                                                                                                            <  NB501016:282:H2V7GAFX2:3:21601:2211:3254/150‑1 (MQ=60)
                   CCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGA                                                                                        <  NB501016:282:H2V7GAFX2:1:11305:14378:13171/151‑1 (MQ=60)
                       CACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTG                                                                                    >  NB501016:282:H2V7GAFX2:1:11107:4480:12006/1‑151 (MQ=60)
                                         CTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAC                                                                  >  NB501016:282:H2V7GAFX2:3:21405:10254:2631/1‑151 (MQ=60)
                                                   GACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                        <  NB501016:282:H2V7GAFX2:3:21505:6235:20000/151‑1 (MQ=60)
                                                        TGAAAGCAAAAAACCGCAACCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTA                                                    <  NB501016:282:H2V7GAFX2:1:21108:3371:9706/150‑1 (MQ=30)
                                                           ACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACATTTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTC                                                >  NB501016:282:H2V7GAFX2:1:11210:5234:16588/1‑151 (MQ=44)
                                                              CAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAGTCAC                                              >  NB501016:282:H2V7GAFX2:4:21501:13104:6099/1‑150 (MQ=39)
                                                               AAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACAC                                            <  NB501016:282:H2V7GAFX2:2:21111:16693:15163/151‑1 (MQ=38)
                                                                     CCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGG                                         >  NB501016:282:H2V7GAFX2:3:21402:16371:8260/1‑148 (MQ=26)
                                                                       GCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATA                                    <  NB501016:282:H2V7GAFX2:3:11505:17286:13832/151‑1 (MQ=22)
                                                                                                CACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGC           >  NB501016:282:H2V7GAFX2:4:11605:11824:15091/1‑151 (MQ=5)
                                                                                                    AGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGTCCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGA       <  NB501016:282:H2V7GAFX2:4:21510:25013:2569/151‑1 (MQ=11)
                                                                                                    AGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTG        <  NB501016:282:H2V7GAFX2:2:21102:5785:16563/150‑1 (MQ=9)
                                                                                                        GTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACT   <  NB501016:282:H2V7GAFX2:3:11412:10595:18874/151‑1 (MQ=22)
                                                                                                        GTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACT   >  NB501016:282:H2V7GAFX2:3:21405:4336:17169/1‑151 (MQ=22)
                                                                                                         TCTCTATGGAGGTCAGGTATGGTTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG  >  NB501016:282:H2V7GAFX2:3:21407:22788:18617/1‑151 (MQ=23)
                                                                                                         TCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG  >  NB501016:282:H2V7GAFX2:4:11604:11463:13342/1‑151 (MQ=24)
                                                                                                                                       |                                                                                                                        
GGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG  >  NC_002947_CJ‑RC/4767849‑4768104

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: