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 F1 I1 R1 152 24.9 2083062 91.6% 1908084 109.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,056,947 +G intergenic (+211/‑15) h16_A0964 → / → h16_A0965 conserved hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,056,9451.G100.0% 58.8 / NA 16intergenic (+209/‑17)h16_A0964/h16_A0965conserved hypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (9/7);  total (9/7)

CACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACG  >  AM260479/1056835‑1057070
                                                                                                               |                                                                                                                             
cacaAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                                           <  2:866625/132‑1 (MQ=255)
cacaAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcg                                                                                                         >  1:630727/1‑134 (MQ=255)
    aaTGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATAGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcgcg                                                                                                     >  1:932602/1‑134 (MQ=255)
            ccgccgGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCAc                                                                                             >  2:70237/1‑134 (MQ=255)
              gccgGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGa                                                                                           <  2:630727/134‑1 (MQ=255)
                    ggcgCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCAt                                                                                     <  2:356852/134‑1 (MQ=255)
                                               acacGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGc                                                                                 <  1:409341/111‑1 (MQ=255)
                                               acacGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGc                                                                                 >  2:409341/1‑111 (MQ=255)
                                                              gggcAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATcag                                           <  2:932602/134‑1 (MQ=255)
                                                                  aaGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCagcagc                                       >  1:551268/1‑134 (MQ=255)
                                                                                 gTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAg                        <  2:551268/134‑1 (MQ=255)
                                                                                             cccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTg               >  1:666627/1‑131 (MQ=255)
                                                                                             cccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTg               <  2:666627/131‑1 (MQ=255)
                                                                                                     gAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTa    >  1:864566/1‑134 (MQ=255)
                                                                                                     gAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTa    >  2:893444/1‑134 (MQ=255)
                                                                                                       gcgcGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACg  >  1:488038/1‑134 (MQ=255)
                                                                                                               |                                                                                                                             
CACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACG  >  AM260479/1056835‑1057070

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCGGCGCAGTGCATTAGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACGG  >  AM260479/1056800‑1057071
                                                                                                                                                    |                                                                                                                            
gtgggctcggagatgtgtataagagaCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGC                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2213:9291:93344/124‑1 (MQ=60)
                    CTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGC                                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:2108:9050:44715/1‑150 (MQ=60)
                        ttgttttttctcgcgttctttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGC                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1214:15062:88219/51‑1 (MQ=60)
                        TGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATAGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGG                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2210:3532:4346/1‑150 (MQ=60)
                                                          aagagacagGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGAC                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1211:14218:79740/141‑1 (MQ=60)
                                                                                      GTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCC                                       >  GWNJ‑0478:712:GW2002102894th:2:2103:4398:99468/1‑150 (MQ=60)
                                                                                                                 AGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGC            >  GWNJ‑0478:712:GW2002102894th:2:2111:12731:63514/1‑150 (MQ=60)
                                                                                                                         GTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTAC    >  GWNJ‑0478:712:GW2002102894th:2:2209:13481:84105/1‑150 (MQ=60)
                                                                                                                           ggcctcgtgcgcaccGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACGG  >  GWNJ‑0478:712:GW2002102894th:2:1214:19780:7333/16‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                            
GCGGCGCAGTGCATTAGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACGG  >  AM260479/1056800‑1057071

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