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 F14 I0 R1 248 61.6 2601574 98.7% 2567753 146.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,767,984 T→G 100% 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% 189.5 / NA 55intergenic (+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 (25/30);  total (25/30)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG  >  NC_002947_CJ‑RC/4767841‑4768124
                                                                                                                                               |                                                                                                                                            
acgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAAt                                                                                                                                         <  2:93593/149‑1 (MQ=255)
   tgtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTaa                                                                                                                                      <  8:171438/149‑1 (MQ=255)
   tgtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTaa                                                                                                                                      <  8:182206/149‑1 (MQ=255)
        ggTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatcatc                                                                                                                                 >  4:111835/1‑149 (MQ=255)
            aCCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGc                                                                                                                              <  4:171917/148‑1 (MQ=255)
                       ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAtgt                                                                                                                   >  6:307520/1‑148 (MQ=255)
                          aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGAt                                                                                                                <  6:146807/148‑1 (MQ=255)
                              tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAc                                                                                                           >  4:235712/1‑149 (MQ=255)
                                  cGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTg                                                                                                        >  5:77453/1‑148 (MQ=255)
                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAAtt                                                                                                   <  6:579/149‑1 (MQ=255)
                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAAtt                                                                                                   >  7:34670/1‑149 (MQ=255)
                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAAtt                                                                                                   >  3:82431/1‑149 (MQ=255)
                                        gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                 <  5:82246/149‑1 (MQ=255)
                                        gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                 <  8:34670/149‑1 (MQ=255)
                                        gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                 <  1:93891/149‑1 (MQ=255)
                                          aCAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                    >  2:305639/1‑112 (MQ=255)
                                          aCAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                    <  1:305639/112‑1 (MQ=255)
                                           cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAg                                                                                              <  3:244672/149‑1 (MQ=255)
                                            aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                             <  7:163440/149‑1 (MQ=255)
                                               ttCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGa                                                                                          >  2:191510/1‑149 (MQ=255)
                                                    aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATaa                                                                                       <  7:321145/147‑1 (MQ=255)
                                                    aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATaa                                                                                       <  8:302482/147‑1 (MQ=255)
                                                    aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAc                                                                                      >  8:277720/1‑148 (MQ=255)
                                                       gCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAAtt                                                                                  >  5:12662/1‑149 (MQ=255)
                                                        ccTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTc                                                                                 >  5:34701/1‑149 (MQ=255)
                                                          tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCt                                                                                >  8:331117/1‑148 (MQ=255)
                                                            aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCtt                                                                               <  4:236083/147‑1 (MQ=255)
                                                             cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAg                                                                            <  6:77453/149‑1 (MQ=255)
                                                                  aaCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTaa                                                                       <  3:257237/149‑1 (MQ=255)
                                                                   aCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAc                                                                      >  4:87195/1‑149 (MQ=255)
                                                                   aCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAc                                                                      <  2:153747/149‑1 (MQ=255)
                                                                       aaaaaaCCGCACCTGGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                   >  3:143404/1‑148 (MQ=14)
                                                                       aaaaaaCCGCACCTGGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                   >  6:154943/1‑148 (MQ=14)
                                                                        aaaaaCCGCACCTGGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGTATAACAATTCTTAAGATTAAATaaca                                                                 >  3:299188/1‑149 (MQ=14)
                                                                               gCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGgaga                                                          <  6:158917/149‑1 (MQ=255)
                                                                                cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGat                                                         <  7:331117/149‑1 (MQ=255)
                                                                                   ccAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGatat                                                       >  3:168805/1‑148 (MQ=255)
                                                                                    cAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCa                                                     >  8:42104/1‑149 (MQ=255)
                                                                                        tGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGAt                                                                                                                <  4:311490/86‑1 (MQ=14)
                                                                                        tGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGAt                                                                                                                >  3:311490/1‑86 (MQ=14)
                                                                                           ggTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTcc                                              >  4:160684/1‑149 (MQ=255)
                                                                                             ttttttGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCtt                                            <  5:307520/149‑1 (MQ=255)
                                                                                                   gAATTCACTAGTAGTCTCTATTGATGTCAGGTATTATTACTATTGACAATTAATCATCGTCTCTTATAATGTGATGAGACCTGGATTTGTGAGCGGATAACAATTCTTAAGATTAACTCACAGAGGATATATCATATGTCCTTAGCAc                                       <  4:299188/148‑1 (MQ=11)
                                                                                                     aTTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCAt                                    <  1:191510/149‑1 (MQ=255)
                                                                                                             gtagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGt                             <  4:194579/148‑1 (MQ=255)
                                                                                                                gtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAAATCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGaca                         >  1:190102/1‑149 (MQ=255)
                                                                                                                 tctcTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGacac                        >  3:193959/1‑149 (MQ=255)
                                                                                                                  ctctATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACt                       >  8:322599/1‑149 (MQ=255)
                                                                                                                     tATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACt                       <  3:235712/146‑1 (MQ=255)
                                                                                                                         gAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAt                >  1:278067/1‑149 (MQ=255)
                                                                                                                          aGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAGAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAtt               <  3:162592/149‑1 (MQ=255)
                                                                                                                               aGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATc          >  6:6331/1‑149 (MQ=255)
                                                                                                                               aGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATc          <  4:143404/149‑1 (MQ=255)
                                                                                                                                       ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgcg  <  5:154943/149‑1 (MQ=255)
                                                                                                                                       ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgc   <  7:64267/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
ACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG  >  NC_002947_CJ‑RC/4767841‑4768124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ 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
CCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGT  >  NC_002947_CJ‑RC/4767874‑4768134
                                                                                                              |                                                                                                                                                       
CCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGG                                                                                                                 >  NB501016:282:H2V7GAFX2:3:11506:7515:3391/1‑150 (MQ=60)
    CAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTG                                                                                                            >  NB501016:282:H2V7GAFX2:2:11207:23251:12000/1‑151 (MQ=60)
    CAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTG                                                                                                            >  NB501016:282:H2V7GAFX2:4:11408:2023:5413/1‑151 (MQ=60)
      AGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTG                                                                                                          <  NB501016:282:H2V7GAFX2:1:11208:6929:12448/151‑1 (MQ=60)
      AGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTG                                                                                                          <  NB501016:282:H2V7GAFX2:3:11507:13434:3812/151‑1 (MQ=60)
        TACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCA                                                                                                                                             <  NB501016:282:H2V7GAFX2:1:21308:12252:6525/114‑1 (MQ=60)
         ACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGC                                                                                                       <  NB501016:282:H2V7GAFX2:2:21209:20000:14306/151‑1 (MQ=60)
          CAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCG                                                                                                      <  NB501016:282:H2V7GAFX2:4:11611:18119:19172/151‑1 (MQ=60)
                  GATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAC                                                                                                <  NB501016:282:H2V7GAFX2:4:21609:13891:13715/149‑1 (MQ=60)
                     CGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTC                                                                                           >  NB501016:282:H2V7GAFX2:3:11403:25703:16494/1‑151 (MQ=60)
                      GCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCT                                                                                          >  NB501016:282:H2V7GAFX2:3:11408:9808:15471/1‑151 (MQ=60)
                                GAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAC                                                                                <  NB501016:282:H2V7GAFX2:2:21212:25096:14312/151‑1 (MQ=50)
                                     CAAAAAACCGCACCTGGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCAC                                                                            >  NB501016:282:H2V7GAFX2:2:11308:11299:14199/1‑150 (MQ=19)
                                      AAAAAACCGCACCTGGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGTATAACAATTCTTAAGATTAAATAACAC                                                                          >  NB501016:282:H2V7GAFX2:2:21310:16074:7166/1‑151 (MQ=16)
                                              GCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATA                                                                  <  NB501016:282:H2V7GAFX2:4:21611:14604:17891/151‑1 (MQ=22)
                                                 CCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATC                                                                >  NB501016:282:H2V7GAFX2:2:21102:22976:4106/1‑150 (MQ=16)
                                                      GTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATC                                                                                                                         >  NB501016:282:H2V7GAFX2:2:21312:9781:18331/1‑88 (MQ=4)
                                                           GTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTA                                                     <  NB501016:282:H2V7GAFX2:3:21608:16370:9450/151‑1 (MQ=5)
                                                                   AATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATC                                             <  NB501016:282:H2V7GAFX2:1:21106:17134:2766/151‑1 (MQ=5)
                                                                              AGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAAATCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACAC                                  >  NB501016:282:H2V7GAFX2:1:21105:8873:16097/1‑151 (MQ=19)
                                                                               GTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACT                                 >  NB501016:282:H2V7GAFX2:2:21108:3700:5167/1‑151 (MQ=22)
                                                                                   CTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG                                <  NB501016:282:H2V7GAFX2:2:21207:15129:14303/148‑1 (MQ=24)
                                                                                       GGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATT                         >  NB501016:282:H2V7GAFX2:1:21301:18424:4519/1‑151 (MQ=38)
                                                                                        GAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAGAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTT                        <  NB501016:282:H2V7GAFX2:2:11312:13279:14723/151‑1 (MQ=39)
                                                                                                     ATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGG           <  NB501016:282:H2V7GAFX2:3:11610:8374:7561/151‑1 (MQ=60)
                                                                                                     ATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG            <  NB501016:282:H2V7GAFX2:4:11502:14323:6926/150‑1 (MQ=60)
                                                                                                            ATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTG     <  NB501016:282:H2V7GAFX2:4:21511:17306:18524/151‑1 (MQ=60)
                                                                                                            ATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAG       <  NB501016:282:H2V7GAFX2:2:21110:11945:9375/149‑1 (MQ=60)
                                                                                                             TGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAG   <  NB501016:282:H2V7GAFX2:2:11301:19621:20030/150‑1 (MQ=60)
                                                                                                              GACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGAACCATCGCGCGTGACACTGCAGGATTTCATCGATTCGCGGCCAGTGAGT  >  NB501016:282:H2V7GAFX2:4:21401:1318:4339/1‑151 (MQ=60)
                                                                                                              |                                                                                                                                                       
CCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGT  >  NC_002947_CJ‑RC/4767874‑4768134

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