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
A3 F14 I1 R1 57 58.2 2652386 98.6% 2615252 146.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,768,005 T→A 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,0050TA100.0% 200.1 / NA 59intergenic (+112/‑70)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 A (23/36);  total (23/36)

AACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACA  >  NC_002947_CJ‑RC/4767861‑4768147
                                                                                                                                                |                                                                                                                                              
aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAt                                                                                                                                            <  6:167345/149‑1 (MQ=255)
aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAt                                                                                                                                            <  3:150589/149‑1 (MQ=255)
 aCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAtg                                                                                                                                           >  6:229970/1‑149 (MQ=255)
   ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAtgt                                                                                                                                          >  7:182519/1‑148 (MQ=255)
          tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAc                                                                                                                                  >  1:288666/1‑149 (MQ=255)
          tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAc                                                                                                                                  <  7:164188/149‑1 (MQ=255)
                    gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACTGCACCCAGGTGCTGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                        <  4:244790/149‑1 (MQ=21)
                    gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                        <  8:14776/149‑1 (MQ=255)
                    gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                        <  3:280471/149‑1 (MQ=255)
                        aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                                                    <  4:64505/149‑1 (MQ=255)
                        aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                                                    <  3:87033/149‑1 (MQ=255)
                             cTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGAt                                                                                                                >  2:201477/1‑148 (MQ=255)
                               gATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATaa                                                                                                              >  7:116838/1‑148 (MQ=255)
                                 tCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACaa                                                                                                           <  4:278137/149‑1 (MQ=255)
                                      tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCt                                                                                                       >  6:330660/1‑148 (MQ=255)
                                      tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCt                                                                                                       >  6:311738/1‑148 (MQ=255)
                                        aCGTTTAACTCAAAAAACCGCTCCCATTTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGATGTCAGGTATTATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                                    <  5:171906/149‑1 (MQ=14)
                                        aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                                    <  8:209798/149‑1 (MQ=255)
                                        aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                                    <  8:8999/149‑1 (MQ=255)
                                         cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAg                                                                                                   <  4:248378/149‑1 (MQ=255)
                                            tGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGAt                                                                                                 >  4:257006/1‑148 (MQ=255)
                                             gAACGCTAAAAACCGCACCCAGTTGCGTTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGTTAACAATTCTTAAGAtt                                                                                                <  3:272567/148‑1 (MQ=14)
                                              aaCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTaa                                                                                              <  6:11441/149‑1 (MQ=255)
                                               aCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTaa                                                                                              <  7:230449/148‑1 (MQ=255)
                                               aCGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAc                                                                                             <  5:311738/149‑1 (MQ=255)
                                                   aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                                          >  2:82398/1‑148 (MQ=255)
                                                    aaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcaca                                                                                        <  8:210995/149‑1 (MQ=255)
                                                     aaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcaca                                                                                        >  7:312336/1‑148 (MQ=255)
                                                      aaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTcacaca                                                                                      <  7:827/149‑1 (MQ=255)
                                                                cAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCa                                                                            >  6:133095/1‑149 (MQ=255)
                                                                 aGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCa                                                                            >  6:38573/1‑148 (MQ=255)
                                                                    tGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATg                                                                        >  1:304664/1‑149 (MQ=255)
                                                                      cGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTc                                                                      <  5:214791/149‑1 (MQ=255)
                                                                       ggTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTcc                                                                     <  6:201834/149‑1 (MQ=255)
                                                                         ttttttGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCtt                                                                   <  8:326281/149‑1 (MQ=255)
                                                                           ttttGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTa                                                                  <  8:182519/148‑1 (MQ=255)
                                                                                           agtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGa                                                  >  7:250953/1‑148 (MQ=255)
                                                                                            gtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGaca                                                <  5:330660/149‑1 (MQ=255)
                                                                                            gtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGac                                                 <  5:229970/148‑1 (MQ=255)
                                                                                               tctATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTg                                             >  4:306711/1‑149 (MQ=255)
                                                                                                ctATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGc                                            >  1:204182/1‑149 (MQ=255)
                                                                                                   tGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCgg                                         <  5:243143/149‑1 (MQ=255)
                                                                                                      aGGTTAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAtt                                      >  6:38509/1‑149 (MQ=255)
                                                                                                           aGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCAt                                  <  8:115226/148‑1 (MQ=255)
                                                                                                                   ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgcg                         <  1:199003/149‑1 (MQ=255)
                                                                                                                   ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgcg                         <  8:77398/149‑1 (MQ=255)
                                                                                                                      ctattgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGc                       <  4:27174/142‑1 (MQ=255)
                                                                                                                      ctattgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGc                       <  3:257006/142‑1 (MQ=255)
                                                                                                                         ttgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGt                   >  8:150532/4‑149 (MQ=255)
                                                                                                                         ttgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGt                   >  1:169500/4‑149 (MQ=255)
                                                                                                                         ttgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGt                   >  4:107278/4‑149 (MQ=255)
                                                                                                                             cAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGt               <  4:156239/149‑1 (MQ=255)
                                                                                                                              aaTTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTc              >  6:47119/1‑149 (MQ=255)
                                                                                                                                  aaTCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGAt           <  7:80994/148‑1 (MQ=255)
                                                                                                                                   atcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc         <  1:95840/149‑1 (MQ=255)
                                                                                                                                   atcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc         <  2:288666/149‑1 (MQ=255)
                                                                                                                                   atcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc         >  5:291781/1‑149 (MQ=255)
                                                                                                                                      atcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCa       <  8:116838/148‑1 (MQ=255)
                                                                                                                                          gCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACa  >  4:32296/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
AACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACA  >  NC_002947_CJ‑RC/4767861‑4768147

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
CGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC  >  NC_002947_CJ‑RC/4767858‑4768141
                                                                                                                                                    |                                                                                                                                        
CGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATA                                                                                                                                        >  NB501016:282:H2V7GAFX2:2:11310:8146:3688/1‑150 (MQ=60)
  CAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATG                                                                                                                                     <  NB501016:282:H2V7GAFX2:2:11309:15940:1213/151‑1 (MQ=60)
     CTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTG                                                                                                                                   >  NB501016:282:H2V7GAFX2:4:21403:17342:3892/1‑150 (MQ=60)
            ATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACC                                                                                                                           >  NB501016:282:H2V7GAFX2:1:21303:13128:8005/1‑151 (MQ=60)
            ATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACC                                                                                                                           <  NB501016:282:H2V7GAFX2:4:11611:6540:5038/151‑1 (MQ=60)
                      AGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTG                                                                                                                 <  NB501016:282:H2V7GAFX2:2:21304:1604:1484/151‑1 (MQ=60)
                          CAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCG                                                                                                             <  NB501016:282:H2V7GAFX2:2:11207:17413:17938/151‑1 (MQ=60)
                                 TGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAC                                                                                                       >  NB501016:282:H2V7GAFX2:4:11512:4946:17267/1‑150 (MQ=60)
                                      GCCTGACGTTGAAGGCAAAATCCCGCACCCTGGTGCGGTTTTTTGATTTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCT                                                                                                 <  NB501016:282:H2V7GAFX2:1:21302:9841:2961/151‑1 (MQ=17)
                                          GACGTTTAACTCAAAAAACCGCTCCCATTTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGATGTCAGGTATTATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                             <  NB501016:282:H2V7GAFX2:3:21401:20314:8283/151‑1 (MQ=7)
                                               TGAACGCTAAAAACCGCACCCAGTTGCGTTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGTTAACAATTCTTAAGATTA                                                                                         <  NB501016:282:H2V7GAFX2:2:21302:7273:5243/150‑1 (MQ=19)
                                                 AACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACT                                                                                      <  NB501016:282:H2V7GAFX2:3:21608:25298:10082/151‑1 (MQ=48)
                                                 AACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAC                                                                                       <  NB501016:282:H2V7GAFX2:4:21412:2155:19538/150‑1 (MQ=48)
                                                       AAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACAC                                                                                 >  NB501016:282:H2V7GAFX2:4:21606:22572:10194/1‑150 (MQ=36)
                                                        AAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAG                                                                               <  NB501016:282:H2V7GAFX2:4:11401:12910:4390/151‑1 (MQ=34)
                                                                      GTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGT                                                                 >  NB501016:282:H2V7GAFX2:1:21307:6090:8227/1‑151 (MQ=5)
                                                                        GCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCC                                                               <  NB501016:282:H2V7GAFX2:3:21411:4131:11573/151‑1 (MQ=5)
                                                                                             TAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGAC                                           >  NB501016:282:H2V7GAFX2:4:21505:10834:6478/1‑150 (MQ=16)
                                                                                              AGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACAC                                         <  NB501016:282:H2V7GAFX2:3:21612:22170:17378/151‑1 (MQ=20)
                                                                                              AGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACA                                          <  NB501016:282:H2V7GAFX2:3:21502:15250:18772/150‑1 (MQ=18)
                                                                                                  TCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCC                                     >  NB501016:282:H2V7GAFX2:1:21108:17468:3829/1‑151 (MQ=28)
                                                                                                     ATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGA                                  <  NB501016:282:H2V7GAFX2:3:21505:7443:18160/151‑1 (MQ=34)
                                                                                                                     ATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGG                  <  NB501016:282:H2V7GAFX2:1:21106:14864:19985/151‑1 (MQ=60)
                                                                                                                        ACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCC                <  NB501016:282:H2V7GAFX2:2:21210:10441:8238/150‑1 (MQ=60)
                                                                                                                            ATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTG            >  NB501016:282:H2V7GAFX2:1:11312:3201:18698/1‑151 (MQ=60)
                                                                                                                                     TAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATA    <  NB501016:282:H2V7GAFX2:4:11505:14967:10145/150‑1 (MQ=60)
                                                                                                                                      AATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC  >  NB501016:282:H2V7GAFX2:3:21604:9527:7981/1‑151 (MQ=60)
                                                                                                                                      AATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC  <  NB501016:282:H2V7GAFX2:1:11208:19976:11956/151‑1 (MQ=60)
                                                                                                                                                    |                                                                                                                                        
CGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC  >  NC_002947_CJ‑RC/4767858‑4768141

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: