Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,483,721 G→C R156R (CGC→CGG pabA ← aminodeoxychorismate synthase component 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,483,7210GC100.0% 97.1 / NA 30R156R (CGC→CGGpabAaminodeoxychorismate synthase component 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/16);  total (14/16)

CAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACC  >  NZ_CP009273/3483634‑3483799
                                                                                       |                                                                              
cAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACtggcgatgccg                                                                              >  1:474197/1‑90 (MQ=255)
 aGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACtggcgatgccga                                                                             <  1:204653/90‑1 (MQ=255)
  ggAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGaa                                                                            <  2:612607/90‑1 (MQ=255)
   gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  2:301661/90‑1 (MQ=255)
   gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  2:49364/90‑1 (MQ=255)
   gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  2:283502/90‑1 (MQ=255)
   gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  1:286714/90‑1 (MQ=255)
    aaaTTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATc                                                                          <  2:400941/90‑1 (MQ=255)
       ttAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATcccc                                                                       >  1:489555/1‑90 (MQ=255)
             agcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAAtc                                                                 <  2:567429/90‑1 (MQ=255)
             agcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAAtc                                                                 <  2:541970/90‑1 (MQ=255)
              gcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAAtct                                                                >  1:249009/1‑90 (MQ=255)
                   tGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCgg                                                           >  2:497144/1‑90 (MQ=255)
                        tCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTc                                                      <  1:239896/90‑1 (MQ=255)
                        tCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTc                                                      <  1:541140/90‑1 (MQ=255)
                           ttGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCt                                                   >  1:319141/1‑90 (MQ=255)
                                                tggatggaACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCa                              >  2:501878/1‑90 (MQ=255)
                                                   atggaACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGc                           <  2:401808/90‑1 (MQ=255)
                                                         cTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTaa                     >  2:537211/1‑90 (MQ=255)
                                                         cTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTaa                     <  2:489555/90‑1 (MQ=255)
                                                           gCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGTTAATg                   >  2:44086/1‑90 (MQ=255)
                                                           gCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATg                   <  2:603057/90‑1 (MQ=255)
                                                               acCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTc               <  1:78007/90‑1 (MQ=255)
                                                                ccTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGc                               >  2:417428/1‑73 (MQ=255)
                                                                ccTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGc                               <  1:417428/73‑1 (MQ=255)
                                                                ccTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCa              >  1:359142/1‑90 (MQ=255)
                                                                  ttCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAgg            >  1:415342/1‑90 (MQ=255)
                                                                          cccACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccacca    >  2:515105/1‑90 (MQ=255)
                                                                          cccACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccacca    >  2:165162/1‑90 (MQ=255)
                                                                            cACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccaccacc  >  2:287069/1‑90 (MQ=255)
                                                                              cTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTcc                                                 <  2:548021/41‑1 (MQ=255)
                                                                              cTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTcc                                                 >  1:548021/1‑41 (MQ=255)
                                                                                       |                                                                              
CAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACC  >  NZ_CP009273/3483634‑3483799

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCAC  >  NZ_CP009273/3483624‑3483819
                                                                                                 |                                                                                                  
ATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCG                                                                                                  >  SRR3722088.479824/1‑100 (MQ=60)
       ATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCC                                                                                           >  SRR3722088.495389/1‑100 (MQ=60)
           AGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAA                                                                                       <  SRR3722088.206704/100‑1 (MQ=60)
             GAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATC                                                                                     <  SRR3722088.289472/100‑1 (MQ=60)
              AAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCT                                                                                    >  SRR3722088.251443/1‑100 (MQ=60)
                           GTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCT                                                                       >  SRR3722088.322381/1‑100 (MQ=60)
                                  TCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCC                                                                <  SRR3722088.242258/100‑1 (MQ=60)
                                  TCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCC                                                                <  SRR3722088.547734/100‑1 (MQ=60)
                                               gtgggctcggagatgtgtataagagacagTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGC                                                   <  SRR3722088.318754/71‑1 (MQ=60)
                                                         gagacaGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAA                                         <  SRR3722088.422061/94‑1 (MQ=60)
                                                                GAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCA                                  >  SRR3722088.362950/1‑100 (MQ=60)
                                                                  ACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGG                                >  SRR3722088.419939/1‑100 (MQ=60)
                                                                         ACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACC                         <  SRR3722088.78798/100‑1 (MQ=60)
                                                                              CCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACctgtctcttatacacatctccgagcccacgagaccgtac                    >  SRR3722088.554708/1‑61 (MQ=60)
                                                                                               TGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCA   >  SRR3722088.56441/1‑100 (MQ=60)
                                                                                                GCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCAC  >  SRR3722088.519306/1‑100 (MQ=60)
                                                                                                 |                                                                                                  
ATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCAC  >  NZ_CP009273/3483624‑3483819

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: