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 4,108,830 A→G H143H (CAT→CAC rraA ← ribonuclease E activity regulator RraA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,108,8300AG100.0% 76.0 / NA 24H143H (CAT→CACrraAribonuclease E activity regulator RraA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/11);  total (13/11)

AAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACC  >  NZ_CP009273/4108746‑4108907
                                                                                    |                                                                             
aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTc                                                                          <  2:274654/90‑1 (MQ=255)
     gCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCgg                                                                     <  1:505900/90‑1 (MQ=255)
       aCCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCg                                                                      <  1:460994/87‑1 (MQ=255)
       aCCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCg                                                                      >  2:460994/1‑87 (MQ=255)
        ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                  <  2:39800/90‑1 (MQ=255)
          ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                <  2:51255/90‑1 (MQ=255)
                    tttCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGacac                                                      >  2:28659/1‑90 (MQ=255)
                    tttCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGacac                                                      >  1:122449/1‑90 (MQ=255)
                           tCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                               >  1:17705/1‑90 (MQ=255)
                                  atatCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcg                                        <  1:584811/90‑1 (MQ=255)
                                   tatCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                       <  2:541874/90‑1 (MQ=255)
                                    atCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgcg                                      >  1:575085/1‑90 (MQ=255)
                                     tCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCgg                                     >  1:40178/1‑90 (MQ=255)
                                       cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                   >  1:603786/1‑90 (MQ=255)
                                        aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                <  1:287778/60‑1 (MQ=255)
                                        aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                >  2:287778/1‑60 (MQ=255)
                                         gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                 >  2:97958/1‑90 (MQ=255)
                                         gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                 <  2:173495/90‑1 (MQ=255)
                                         gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                 <  1:589001/90‑1 (MQ=255)
                                            gATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                              >  1:36916/1‑90 (MQ=255)
                                                ttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                          >  2:43938/1‑90 (MQ=255)
                                                     aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                     <  2:631023/90‑1 (MQ=255)
                                                        gAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGc                  >  2:486475/1‑90 (MQ=255)
                                                                        gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc  >  1:589581/1‑90 (MQ=255)
                                                                                    |                                                                             
AAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACC  >  NZ_CP009273/4108746‑4108907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG  >  NZ_CP009273/4108731‑4108921
                                                                                                   |                                                                                           
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                             <  SRR3722088.571686/100‑1 (MQ=60)
                   GGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGG                                                                          <  SRR3722088.466387/100‑1 (MQ=60)
                    GCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGT                                                                         <  SRR3722088.511982/100‑1 (MQ=60)
                         gagatgtgtataagagacagTTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                    <  SRR3722088.290540/80‑1 (MQ=60)
                         TTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                    >  SRR3722088.123721/1‑100 (MQ=60)
                                GCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAA                                                             >  SRR3722088.17898/1‑100 (MQ=60)
                                         ATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCG                                                    >  SRR3722088.582135/1‑100 (MQ=60)
                                          TCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGG                                                   >  SRR3722088.40598/1‑100 (MQ=60)
                                            ATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAC                                                 >  SRR3722088.611222/1‑100 (MQ=60)
                                                 ATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGC                                            >  SRR3722088.37302/1‑100 (MQ=60)
                                                 ATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGC                                            <  SRR3722088.591984/100‑1 (MQ=60)
                                                        GCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCC                                     <  SRR3722088.596231/100‑1 (MQ=60)
                                                                             TCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACC                >  SRR3722088.596817/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCctgtctcttatacacatctccgagcccacgagacc  >  SRR3722088.18473/1‑65 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG  >  SRR3722088.636288/1‑100 (MQ=60)
                                                                                                   |                                                                                           
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG  >  NZ_CP009273/4108731‑4108921

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