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
A5 F13 I0 R1 213 60.4 2554362 98.6% 2518600 147.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,768,005 T→A 100% intergenic (+112/‑70) 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,768,0050TA98.7% 163.8 / ‑6.2 51intergenic (+112/‑70)fpvA/tpaKTonB‑dependent outer membrane ferripyoverdine receptor FpvA/probable terephthalate transporter, MFS superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (23/27);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAG  >  NC_002947_CJ‑RC/4767862‑4768143
                                                                                                                                               |                                                                                                                                          
aCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAtg                                                                                                                                       >  3:9325/1‑149 (MQ=255)
aCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAtg                                                                                                                                       <  5:47055/149‑1 (MQ=255)
  atCACCATCACCGCCAACTACAGCTTCTGATCGCCTGAAGTAGAACGCAAAAAACCGAACACAGGAGAGGTATTTTGAATACACTAGTAGTCTCTATGGAGGTCAGGTAAGATTACTAATGACAATTAATCATCTGCTCGTATaaa                                                                                                                                        >  6:313291/2‑145 (MQ=21)
       cATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACTCAAAAAACCGCACCCAGGTGCGTTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAg                                                                                                                                <  1:87586/149‑1 (MQ=25)
          cACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAcc                                                                                                                             >  7:126977/1‑149 (MQ=255)
                   gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                    >  8:200286/1‑149 (MQ=255)
                   gTACAGCTTCTGATCGCCTAACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                    >  8:195580/1‑149 (MQ=255)
                    tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                    >  6:150123/1‑148 (MQ=255)
                      cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGa                                                                                                                  <  1:253778/148‑1 (MQ=255)
                      cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAg                                                                                                                 <  4:9325/149‑1 (MQ=255)
                       aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                                                <  3:227299/149‑1 (MQ=255)
                       aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                                                <  5:216585/149‑1 (MQ=255)
                              gATCGCCTGACGTTTAACGCAAAAAACCGCACCCAGGTGCGGATAATTCAAATCACTAGTAGTCTCTAAGGAGGACAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAAGTGATCAGACCTGGAATTGTGAGCGGATat                                                                                                          >  8:127289/1‑147 (MQ=17)
                               aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAc                                                                                                         >  2:226245/1‑148 (MQ=255)
                                   ccTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGa                                                                                                                               <  3:31256/122‑1 (MQ=255)
                                   ccTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGa                                                                                                                               >  4:31256/1‑122 (MQ=255)
                                    cTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTATCTATGGAGGTCAGGTATGATTACTATTGACATTTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAGAAttat                                                                                                   >  2:261792/1‑147 (MQ=18)
                                     tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCt                                                                                                   >  6:65463/1‑148 (MQ=255)
                                       tcGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                                <  1:290662/148‑1 (MQ=255)
                                       aCGTTTAACTCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGATGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCTGATAACAATTCTTaa                                                                                                <  7:144125/149‑1 (MQ=17)
                                       aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                                <  5:234537/149‑1 (MQ=255)
                                             aaCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTaa                                                                                          <  8:15997/149‑1 (MQ=255)
                                               cGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACt                                                                                        <  2:211765/149‑1 (MQ=255)
                                                 cAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                                      <  5:150123/149‑1 (MQ=255)
                                                 cAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTCTGAGCGGATAACAATTCTTAAGATTAACTc                                                                                       >  4:115734/1‑148 (MQ=255)
                                                  aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                                      >  3:42520/1‑148 (MQ=255)
                                                   aaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcaca                                                                                    >  3:297899/1‑149 (MQ=255)
                                                      aaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAg                                                                                 >  8:24965/1‑149 (MQ=255)
                                                        ccGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGga                                                                               <  5:186646/149‑1 (MQ=255)
                                                           cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGat                                                                            <  8:126977/149‑1 (MQ=255)
                                                           cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGat                                                                            <  2:317115/149‑1 (MQ=255)
                                                                  gTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCatat                                                                     <  2:238005/149‑1 (MQ=255)
                                                                   tGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATg                                                                    >  5:288130/1‑149 (MQ=255)
                                                                       gTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCt                                                                >  3:25329/1‑149 (MQ=255)
                                                                                     cTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgcgc                                                  <  2:144325/149‑1 (MQ=255)
                                                                                        gtagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgcgcg                                                 <  3:187102/147‑1 (MQ=255)
                                                                                        gtagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGt                                                <  7:113859/148‑1 (MQ=255)
                                                                                             ctctATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACt                                          >  4:47793/1‑149 (MQ=255)
                                                                                                tATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGcc                                       >  5:32145/1‑149 (MQ=255)
                                                                                                tATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGcc                                       >  2:153403/1‑149 (MQ=255)
                                                                                                        tCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCa                               <  1:226245/149‑1 (MQ=255)
                                                                                                            gTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGa                           <  6:77790/149‑1 (MQ=255)
                                                                                                            gTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGa                           <  6:300215/149‑1 (MQ=255)
                                                                                                             tATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGAt                          >  5:206577/1‑149 (MQ=255)
                                                                                                                        ttgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGt               <  7:135572/146‑1 (MQ=255)
                                                                                                                          gaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGt               >  2:298847/2‑147 (MQ=255)
                                                                                                                           aCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAg            >  6:278487/1‑149 (MQ=255)
                                                                                                                           aCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAg            >  1:295711/1‑149 (MQ=255)
                                                                                                                                  atcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc     <  7:200286/149‑1 (MQ=255)
                                                                                                                                  atcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGAATAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc     <  7:195580/149‑1 (MQ=255)
                                                                                                                                     atcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAg  <  3:200384/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                          
ACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAG  >  NC_002947_CJ‑RC/4767862‑4768143

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
CGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACATCGTGAT  >  NC_002947_CJ‑RC/4767858‑4768154
                                                                                                                                                    |                                                                                                                                                     
CGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATA                                                                                                                                                     >  NB501016:282:H2V7GAFX2:2:11206:14644:6987/1‑150 (MQ=60)
   AACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTTAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGT                                                                                                                                                 <  NB501016:282:H2V7GAFX2:3:11411:25861:16672/151‑1 (MQ=60)
   AACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGT                                                                                                                                                 >  NB501016:282:H2V7GAFX2:2:11103:23789:4613/1‑151 (MQ=60)
          CCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACTCAAAAAACCGCACCCAGGTGCGTTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGA                                                                                                                                          <  NB501016:282:H2V7GAFX2:1:11207:4984:10121/151‑1 (MQ=60)
             TCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCT                                                                                                                                       >  NB501016:282:H2V7GAFX2:4:11604:24038:11131/1‑151 (MQ=60)
                         ACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAG                                                                                                                            <  NB501016:282:H2V7GAFX2:1:21208:2062:17543/150‑1 (MQ=60)
                          CAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCG                                                                                                                          <  NB501016:282:H2V7GAFX2:2:21206:16148:7028/151‑1 (MQ=60)
                          CAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCG                                                                                                                          <  NB501016:282:H2V7GAFX2:3:21501:23899:12013/151‑1 (MQ=60)
                                      GCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAC                                                                                                                                         <  NB501016:282:H2V7GAFX2:2:11108:9734:5573/124‑1 (MQ=60)
                                          GTCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                                          <  NB501016:282:H2V7GAFX2:1:21306:22975:1811/151‑1 (MQ=57)
                                          GACGTTTAACTCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGATGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCTGATAACAATTCTTAAG                                                                                                          <  NB501016:282:H2V7GAFX2:4:11608:8592:10172/151‑1 (MQ=37)
                                          GACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                                          <  NB501016:282:H2V7GAFX2:3:21505:7044:14648/151‑1 (MQ=60)
                                                    GCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCAC                                                                                                <  NB501016:282:H2V7GAFX2:3:11609:5571:16864/151‑1 (MQ=42)
                                                     CAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCAC                                                                                                >  NB501016:282:H2V7GAFX2:2:11110:14809:16982/1‑150 (MQ=40)
                                                      AAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTATCTCACAC                                                                                              >  NB501016:282:H2V7GAFX2:2:21310:2522:14503/1‑151 (MQ=37)
                                                           ACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAG                                                                                         <  NB501016:282:H2V7GAFX2:3:21406:12943:8852/151‑1 (MQ=28)
                                                                      GTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGT                                                                              >  NB501016:282:H2V7GAFX2:3:21605:21395:11024/1‑151 (MQ=5)
                                                                          GGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTT                                                                          >  NB501016:282:H2V7GAFX2:2:11106:4597:18066/1‑151 (MQ=5)
                                                                                           AGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTG                                                          <  NB501016:282:H2V7GAFX2:4:11601:17668:11235/150‑1 (MQ=9)
                                                                                           AGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGT                                                           <  NB501016:282:H2V7GAFX2:2:21106:25681:19233/149‑1 (MQ=7)
                                                                                                   CTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCG                                                 >  NB501016:282:H2V7GAFX2:3:11408:15834:8051/1‑151 (MQ=30)
                                                                                                           GTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCAT                                         <  NB501016:282:H2V7GAFX2:1:21202:7801:14242/151‑1 (MQ=46)
                                                                                                                GTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATT                                    >  NB501016:282:H2V7GAFX2:3:21411:4835:3996/1‑151 (MQ=56)
                                                                                                                            ATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTG                         <  NB501016:282:H2V7GAFX2:4:11606:24272:10820/151‑1 (MQ=60)
                                                                                                                              GACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGT                      >  NB501016:282:H2V7GAFX2:1:21307:21758:6317/1‑151 (MQ=60)
                                                                                                                                      AATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC               <  NB501016:282:H2V7GAFX2:4:21408:1172:18151/151‑1 (MQ=60)
                                                                                                                                      AATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGAATAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC               <  NB501016:282:H2V7GAFX2:4:21407:26750:18156/151‑1 (MQ=60)
                                                                                                                                         CATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGT            <  NB501016:282:H2V7GAFX2:2:21110:7843:4377/151‑1 (MQ=60)
                                                                                                                                                    ATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACATCGTGAT  <  NB501016:282:H2V7GAFX2:1:11310:14141:7026/150‑1 (MQ=60)
                                                                                                                                                    |                                                                                                                                                     
CGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACATCGTGAT  >  NC_002947_CJ‑RC/4767858‑4768154

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: