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 F12 I1 R1 58 76.7 3263832 98.3% 3208346 146.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,767,982 A→T intergenic (+89/‑93) 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,9820AT97.0% 200.7 / ‑4.7 67intergenic (+89/‑93)fpvA/tpaKTonB‑dependent outer membrane ferripyoverdine receptor FpvA/probable terephthalate transporter, MFS superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (30/35);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGG  >  NC_002947_CJ‑RC/4767844‑4768125
                                                                                                                                          |                                                                                                                                               
tgtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTaa                                                                                                                                       <  1:356160/149‑1 (MQ=255)
 gtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTa                                                                                                                                        >  8:361015/1‑147 (MQ=255)
 gtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAat                                                                                                                                      >  3:131619/1‑149 (MQ=255)
  tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                     <  6:326534/149‑1 (MQ=255)
  tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAat                                                                                                                                      >  1:156811/1‑148 (MQ=255)
   aTGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                     >  7:222664/1‑148 (MQ=255)
    tGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatca                                                                                                                                    >  2:206889/1‑148 (MQ=255)
    tGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatca                                                                                                                                    >  8:352495/1‑148 (MQ=255)
      gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatcatc                                                                                                                                  <  3:273724/148‑1 (MQ=255)
      gTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatcatc                                                                                                                                  >  7:121471/1‑148 (MQ=255)
       tGACCCAAGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCgg                                                                                                                                <  6:131213/149‑1 (MQ=255)
          cccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTc                                                                                                                             >  5:374545/1‑149 (MQ=255)
          cccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTc                                                                                                                             <  1:73540/149‑1 (MQ=255)
          cccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTc                                                                                                                             <  1:89285/149‑1 (MQ=255)
           ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCg                                                                                                                            >  3:276861/1‑149 (MQ=255)
           ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCg                                                                                                                            >  6:126557/1‑149 (MQ=255)
             aCGCAACTTCAACTTCACCGAAATGTACAGAATCTGATACCATGACGTTGAACGCAAAAAACAGCACCCAGGTGCGGTTTTTTTAATTCACAACTAGTCTCTAGGGAGGACAGGTATGATTAATATTGACAATTAATCATCGACtataa                                                                                                                          >  6:116591/1‑145 (MQ=255)
                cAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATaa                                                                                                                       <  5:225399/149‑1 (MQ=255)
                  aCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAt                                                                                                                      <  8:35323/148‑1 (MQ=255)
                           tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAc                                                                                                            <  4:131619/149‑1 (MQ=255)
                               cGCCAAGTACAGCATCAGATAGCCTGACGATGAACGCAAAAAACCTCACCCAGGTGCGGTTTTTTGAATTCACTATTAGTCACAATGGACGTAACGTATGAATACAAAAGACAATTAATCATCGGCAAGTATAATGTTATCAGACCt                                                                                                          >  8:178120/1‑147 (MQ=255)
                                 ccAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGa                                                                                                       >  7:243405/1‑148 (MQ=255)
                                     gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGATTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATGACAAATAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                  >  2:324234/1‑149 (MQ=255)
                                      tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                  <  7:409571/148‑1 (MQ=255)
                                      tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                  <  8:106007/148‑1 (MQ=255)
                                        cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGa                                                                                                <  4:32932/148‑1 (MQ=255)
                                         aggttCTGTTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGa                                                                                                <  8:120810/144‑1 (MQ=21)
                                          gCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                              <  8:388253/148‑1 (MQ=255)
                                           cTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                     >  3:273397/1‑108 (MQ=255)
                                           cTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAatc                                                                                                                                     <  4:273397/108‑1 (MQ=255)
                                                 aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATaa                                                                                        <  1:143948/147‑1 (MQ=255)
                                                    gCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAAtt                                                                                   <  1:303987/149‑1 (MQ=255)
                                                      cTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCt                                                                                 >  4:104223/1‑149 (MQ=255)
                                                         aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                              <  3:90639/149‑1 (MQ=255)
                                                          cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTa                                                                               <  7:381644/147‑1 (MQ=255)
                                                           gTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCtt                                                                                >  4:128054/1‑145 (MQ=255)
                                                              gAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGAt                                                                           >  7:306168/1‑147 (MQ=255)
                                                                 cGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTaa                                                                        <  7:178120/147‑1 (MQ=255)
                                                                 cGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCa                                                      >  6:318824/1‑149 (MQ=12)
                                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                              >  1:379498/1‑140 (MQ=255)
                                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                              <  2:379498/140‑1 (MQ=255)
                                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATc                                                       >  7:253636/1‑147 (MQ=12)
                                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCat                                                     <  1:206889/149‑1 (MQ=12)
                                                                    aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                    >  4:88209/1‑148 (MQ=255)
                                                                    aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATAATCGGCTCGTATAATGTGATCAGACATGGAATTGTGAGCGGATAAAAATTCTTAAGATTAACTca                                                                    >  4:272024/1‑148 (MQ=17)
                                                                          ccGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCt                                              >  1:102649/1‑148 (MQ=12)
                                                                          ccGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCt                                              >  1:35775/1‑148 (MQ=12)
                                                                             cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTa                                            >  8:159465/1‑147 (MQ=12)
                                                                             cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAg                                           <  8:193119/148‑1 (MQ=12)
                                                                                ccAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTACGGAGGTCAGGTATGATTACGATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATTTCATATGGCCTTAGCAc                                        >  4:320574/1‑148 (MQ=9)
                                                                                   ggTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATc                                    <  1:33563/149‑1 (MQ=12)
                                                                                    gTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCatat                                                   >  3:207447/1‑149 (MQ=255)
                                                                                     tGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATg                                                  >  4:62624/1‑149 (MQ=255)
                                                                                                gAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACt                        <  5:278129/148‑1 (MQ=12)
                                                                                                  aTTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGcc                     <  5:271518/149‑1 (MQ=12)
                                                                                                    tCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATc                                    <  8:222664/148‑1 (MQ=255)
                                                                                                    tCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCATTAGCACCATCGCGCGTGACACTGCCg                    >  5:90328/1‑148 (MQ=9)
                                                                                                      aCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgc                                  <  4:276861/148‑1 (MQ=255)
                                                                                                      aCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAt                 <  2:136241/149‑1 (MQ=12)
                                                                                                          gtagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGt                              >  6:359795/1‑148 (MQ=255)
                                                                                                             gtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGaca                          >  2:229858/1‑149 (MQ=255)
                                                                                                                  tATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTg                       <  5:406250/147‑1 (MQ=255)
                                                                                                                       aGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAtt                <  5:126557/149‑1 (MQ=255)
                                                                                                                          tCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCATTAGCACCATCGCGCGTGACACTGCCGGATTTCa             <  2:305256/149‑1 (MQ=255)
                                                                                                                                    ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgc    <  2:156811/148‑1 (MQ=255)
                                                                                                                                    ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgc    <  3:62624/148‑1 (MQ=255)
                                                                                                                                      aCTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCgg  <  4:113890/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
TGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGG  >  NC_002947_CJ‑RC/4767844‑4768125

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCG  >  NC_002947_CJ‑RC/4767832‑4768136
                                                                                                                                                      |                                                                                                                                                          
GCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCCCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCTCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAT                                                                                                                                                            >  NB501016:282:H2V7GAFX2:1:21310:19078:16238/1‑151 (MQ=60)
  CCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTT                                                                                                                                                          >  NB501016:282:H2V7GAFX2:2:11106:4078:18269/1‑151 (MQ=60)
      ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACA                                                                                                                                                       >  NB501016:282:H2V7GAFX2:1:11211:10823:16952/1‑150 (MQ=60)
           GTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAAT                                                                                                                                                 <  NB501016:282:H2V7GAFX2:1:21303:15195:11342/151‑1 (MQ=60)
            TGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATC                                                                                                                                                >  NB501016:282:H2V7GAFX2:2:11212:13105:5444/1‑151 (MQ=60)
             GTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATC                                                                                                                                                >  NB501016:282:H2V7GAFX2:1:11303:10913:10123/1‑150 (MQ=60)
              TATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCAT                                                                                                                                              >  NB501016:282:H2V7GAFX2:4:21402:20629:18403/1‑151 (MQ=60)
                 GGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGG                                                                                                                                           >  NB501016:282:H2V7GAFX2:4:11509:11666:5605/1‑151 (MQ=60)
                 GGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCG                                                                                                                                            <  NB501016:282:H2V7GAFX2:2:21201:19645:11763/150‑1 (MQ=60)
                     ACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCG                                                                                                                                       >  NB501016:282:H2V7GAFX2:3:21606:24709:18473/1‑151 (MQ=60)
                     ACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCG                                                                                                                                       <  NB501016:282:H2V7GAFX2:1:11203:10553:13643/151‑1 (MQ=60)
                     ACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCG                                                                                                                                       <  NB501016:282:H2V7GAFX2:1:11112:17189:19460/151‑1 (MQ=60)
                      CCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT                                                                                                                                      >  NB501016:282:H2V7GAFX2:2:21202:14460:11241/1‑151 (MQ=60)
                           GCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAT                                                                                                                                 <  NB501016:282:H2V7GAFX2:3:21404:26727:4901/151‑1 (MQ=60)
                                   ACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTatgaattgtgagc                                                                                                                         <  NB501016:282:H2V7GAFX2:4:21508:16620:13498/151‑14 (MQ=60)
                                            GCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAAT                                                                                                                >  NB501016:282:H2V7GAFX2:4:21406:2756:9503/1‑151 (MQ=60)
                                                AGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTG                                                                                                            <  NB501016:282:H2V7GAFX2:4:21611:21634:5079/151‑1 (MQ=60)
                                                    CAGCTTCTGATCGCCTTACGTTGAACGAAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAG                                                                                        <  NB501016:282:H2V7GAFX2:2:11307:3269:17714/151‑1 (MQ=56)
                                                    CAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAG                                                                                        <  NB501016:282:H2V7GAFX2:3:21508:1326:10264/151‑1 (MQ=60)
                                                      GCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCA                                                                                                                                               >  NB501016:282:H2V7GAFX2:2:21201:17176:10548/1‑110 (MQ=60)
                                                         TCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATT                                                                                     <  NB501016:282:H2V7GAFX2:1:21203:16314:12366/149‑1 (MQ=60)
                                                            GATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAC                                                                                                  <  NB501016:282:H2V7GAFX2:1:11301:12920:4463/149‑1 (MQ=60)
                                                               CGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTC                                                                                             <  NB501016:282:H2V7GAFX2:1:21205:8796:11927/151‑1 (MQ=60)
                                                               CGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACA                                                                             <  NB501016:282:H2V7GAFX2:1:21202:8048:12793/151‑1 (MQ=57)
                                                                    GACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                        <  NB501016:282:H2V7GAFX2:2:11205:5300:4642/151‑1 (MQ=60)
                                                                    GACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAA                                                                                         <  NB501016:282:H2V7GAFX2:4:21606:7238:4750/150‑1 (MQ=60)
                                                                       GTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGAT                                                                     <  NB501016:282:H2V7GAFX2:4:21409:2142:10963/151‑1 (MQ=45)
                                                                         TGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTA                                                                                    >  NB501016:282:H2V7GAFX2:4:21505:11291:1791/1‑150 (MQ=52)
                                                                           AACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAC                                                                                  <  NB501016:282:H2V7GAFX2:4:11607:16581:2776/150‑1 (MQ=46)
                                                                             CGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                        >  NB501016:282:H2V7GAFX2:1:21308:21985:3282/1‑142 (MQ=45)
                                                                             CGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATA                                                               <  NB501016:282:H2V7GAFX2:1:11312:12697:12187/151‑1 (MQ=35)
                                                                             CGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCAT                                                                >  NB501016:282:H2V7GAFX2:4:21408:7060:4257/1‑150 (MQ=35)
                                                                                     ACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTT                                                        >  NB501016:282:H2V7GAFX2:1:11106:7243:14561/1‑150 (MQ=30)
                                                                                     ACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTT                                                        >  NB501016:282:H2V7GAFX2:1:11205:3258:18847/1‑150 (MQ=30)
                                                                                              AGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCG                                              <  NB501016:282:H2V7GAFX2:1:11106:16643:7971/151‑1 (MQ=32)
                                                                                               GGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATT                                                             >  NB501016:282:H2V7GAFX2:2:11312:20809:14394/1‑151 (MQ=5)
                                                                                                           TGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG                                  <  NB501016:282:H2V7GAFX2:3:21502:12796:3286/150‑1 (MQ=30)
                                                                                                             AATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCG                               <  NB501016:282:H2V7GAFX2:3:21412:11342:18858/151‑1 (MQ=32)
                                                                                                               TTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCATTAGCACCATCGCGCGTGACACTGCCGG                              >  NB501016:282:H2V7GAFX2:3:11505:7108:4391/1‑150 (MQ=22)
                                                                                                                         GTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGAT                   >  NB501016:282:H2V7GAFX2:4:11402:4790:4954/1‑151 (MQ=43)
                                                                                                                          TCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATT                  >  NB501016:282:H2V7GAFX2:4:11510:15215:11948/1‑151 (MQ=45)
                                                                                                                             CTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGC                                 <  NB501016:282:H2V7GAFX2:3:21612:14001:15358/149‑1 (MQ=26)
                                                                                                                               ATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGCGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG              <  NB501016:282:H2V7GAFX2:2:11201:25744:20208/150‑1 (MQ=50)
                                                                                                                                  GAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTT                          <  NB501016:282:H2V7GAFX2:3:11511:7225:8626/151‑1 (MQ=40)
                                                                                                                                         GGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTC   <  NB501016:282:H2V7GAFX2:4:21602:20814:9925/151‑1 (MQ=60)
                                                                                                                                          GTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAAGAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCG  <  NB501016:282:H2V7GAFX2:2:21201:5719:15219/151‑1 (MQ=60)
                                                                                                                                               ATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG              <  NB501016:282:H2V7GAFX2:2:11112:6246:5504/150‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                                          
GCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCG  >  NC_002947_CJ‑RC/4767832‑4768136

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