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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,613,222 C→G R238P (CGC→CCC)  lplA ← lipoate‑‑protein ligase LplA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,613,2220CG95.5% 137.7 / ‑0.3 44R238P (CGC→CCC) lplAlipoate‑‑protein ligase LplA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (23/19);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

ACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGC  >  NZ_CP009273/4613141‑4613307
                                                                                 |                                                                                     
aCGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAg                                                                               >  1:936218/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTc                                                                          >  1:1235609/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTc                                                                          >  2:139937/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTc                                                                          >  2:1127300/1‑90 (MQ=255)
     gccCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTc                                                                          >  1:245846/1‑90 (MQ=255)
               aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt                                                                >  1:1362770/1‑90 (MQ=255)
               aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt                                                                >  2:1351272/1‑90 (MQ=255)
                 gCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTgg                                                              >  1:1266758/1‑90 (MQ=255)
                   gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                            <  2:936218/90‑1 (MQ=255)
                   gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                            <  2:911077/90‑1 (MQ=255)
                   gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                            >  2:761390/1‑90 (MQ=255)
                    ttCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCaa                                                           >  2:1009173/1‑90 (MQ=255)
                       aTCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTc                                                        <  2:426984/90‑1 (MQ=255)
                          cagcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTgg                                                     <  2:127850/90‑1 (MQ=255)
                           agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGc                                                    >  2:938642/1‑90 (MQ=255)
                           agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGc                                                    >  2:825799/1‑90 (MQ=255)
                            gcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCg                                                   >  1:841906/1‑90 (MQ=255)
                               gATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGttt                                                <  2:1362770/90‑1 (MQ=255)
                                   cGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGtt                                            <  1:139937/90‑1 (MQ=255)
                                   cGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGtt                                            >  1:191255/1‑90 (MQ=255)
                                      gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg                                         >  2:686704/1‑90 (MQ=255)
                                      gaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCgg                                         >  1:279789/1‑90 (MQ=255)
                                         tGCCGGAGCCTGACCGAAGTTCCATTCTCAGCTACTCTGGGGGGCAAAGGGTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGga                                      >  1:1256032/1‑90 (MQ=255)
                                         tGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGga                                      >  2:1207021/1‑90 (MQ=255)
                                          gCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGgag                                     >  2:1567594/1‑90 (MQ=255)
                                              gAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAgatga                                 <  1:1046612/90‑1 (MQ=255)
                                                     aCCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGc                          <  2:852449/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAGGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:807206/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:955444/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  1:1422750/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:465275/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:180084/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:1109161/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:1067578/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  1:825799/90‑1 (MQ=255)
                                                               cATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc                <  2:245846/90‑1 (MQ=255)
                                                                aTTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCt               <  2:1266758/90‑1 (MQ=255)
                                                                   cccAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGc            >  2:1574248/1‑90 (MQ=255)
                                                                     cAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGtttt                                               >  2:896217/1‑53 (MQ=255)
                                                                     cAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGtttt                                               <  1:896217/53‑1 (MQ=255)
                                                                       gCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATa        >  2:601585/1‑90 (MQ=255)
                                                                             ctGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGc  >  1:1547167/1‑90 (MQ=255)
                                                                             ctGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGc  >  1:742495/1‑90 (MQ=255)
                                                                             ctGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGAAGATTTCCGCTTCCACGCGCTCGCCATAATGGGc  <  2:1534533/90‑1 (MQ=255)
                                                                                 |                                                                                     
ACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGC  >  NZ_CP009273/4613141‑4613307

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ 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
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG  >  NZ_CP009273/4613124‑4613314
                                                                                                  |                                                                                            
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGG                                                                                             <  SRR3722078.399347/100‑1 (MQ=60)
       ATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAG                                                                                      >  SRR3722078.949190/1‑100 (MQ=60)
            GTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTC                                                                                 >  SRR3722078.1255118/1‑100 (MQ=60)
            GTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTC                                                                                 >  SRR3722078.248853/1‑100 (MQ=60)
                      GCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGctgtctcttataca                                                                       >  SRR3722078.1529452/1‑86 (MQ=60)
                      GCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTT                                                                       >  SRR3722078.1384685/1‑100 (MQ=60)
                        CCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGG                                                                     >  SRR3722078.1286859/1‑100 (MQ=60)
                                   CGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCG                                                          >  SRR3722078.853133/1‑100 (MQ=60)
                                          GCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTT                                                   >  SRR3722078.193530/1‑100 (MQ=60)
                                             GCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGG                                                >  SRR3722078.283204/1‑100 (MQ=60)
                                                GATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCTCAGCTACTCTGGGGGGCAAAGGGTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGA                                             >  SRR3722078.1275922/1‑100 (MQ=60)
                                                 gcagcgtcagatgtgtataagagacagGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAG                                            <  SRR3722078.908255/73‑1 (MQ=60)
                                                    CGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATG                                         <  SRR3722078.141569/100‑1 (MQ=60)
                                                               GAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTC                              <  SRR3722078.1062269/100‑1 (MQ=60)
                                                                                aagagaCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722078.446169/94‑1 (MQ=60)
                                                                                CATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722078.1445862/100‑1 (MQ=60)
                                                                                CATCCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722078.836745/100‑1 (MQ=60)
                                                                                    CCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGC         >  SRR3722078.1572370/1‑100 (MQ=60)
                                                                                    CCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGC         >  SRR3722078.751892/1‑100 (MQ=60)
                                                                                      CAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGA       >  SRR3722078.1439730/1‑100 (MQ=60)
                                                                                           ACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG  >  SRR3722078.292862/1‑100 (MQ=60)
                                                                                                  |                                                                                            
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG  >  NZ_CP009273/4613124‑4613314

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