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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

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% 116.3 / NA 34H143H (CAT→CACrraAribonuclease E activity regulator RraA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/17);  total (17/17)

AGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCG  >  NZ_CP009273/4108749‑4108912
                                                                                 |                                                                                  
aGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGcc                                                                            <  2:1008896/90‑1 (MQ=255)
 ggCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCg                                                                           <  1:1004774/90‑1 (MQ=255)
  gCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCgg                                                                          <  1:881812/90‑1 (MQ=255)
     ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                       <  1:310847/90‑1 (MQ=255)
       ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     >  1:1132021/1‑90 (MQ=255)
       ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  2:589590/90‑1 (MQ=255)
       ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  1:611581/90‑1 (MQ=255)
       ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  1:661146/90‑1 (MQ=255)
              gCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGa                                                              >  1:893707/1‑90 (MQ=255)
                     ttATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGcc                                                       <  1:44000/90‑1 (MQ=255)
                       aTCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCaa                                                     <  2:630092/90‑1 (MQ=255)
                           ttCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAtt                                                 >  1:1166465/1‑90 (MQ=255)
                            tCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                                <  1:1170765/90‑1 (MQ=255)
                               atatCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcg                                             <  2:143820/90‑1 (MQ=255)
                                tatCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                            <  1:105212/90‑1 (MQ=255)
                                    cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                        >  1:733783/1‑90 (MQ=255)
                                    cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                        >  2:1138144/1‑90 (MQ=255)
                                     aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                       >  1:234928/1‑90 (MQ=255)
                                      gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      >  1:986903/1‑90 (MQ=255)
                                           tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                 >  2:699216/1‑90 (MQ=255)
                                            cttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCttt                                <  2:588122/90‑1 (MQ=255)
                                              tctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCg                              >  2:422770/1‑90 (MQ=255)
                                                  aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                          <  2:235707/90‑1 (MQ=255)
                                                               ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  2:390422/1‑90 (MQ=255)
                                                               ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  1:1157671/1‑90 (MQ=255)
                                                               ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             <  2:585953/90‑1 (MQ=255)
                                                               ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  1:1168273/1‑90 (MQ=255)
                                                               ggTATTGTCGGCATAAAGGGGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  1:755800/1‑90 (MQ=255)
                                                                    tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        <  2:986903/90‑1 (MQ=255)
                                                                     gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc       >  2:392038/1‑90 (MQ=255)
                                                                     gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc       >  2:899266/1‑90 (MQ=255)
                                                                     gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc       >  1:809474/1‑90 (MQ=255)
                                                                         gCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAAcc   <  1:217901/90‑1 (MQ=255)
                                                                          cATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCg  >  2:704825/1‑90 (MQ=255)
                                                                                 |                                                                                  
AGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCG  >  NZ_CP009273/4108749‑4108912

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATCGCC  >  NZ_CP009273/4108731‑4108929
                                                                                                    |                                                                                                   
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722089.742860/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722089.722066/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCAGCATAAAGG                                                                                                     <  SRR3722089.259763/100‑1 (MQ=60)
 GAGGTAACGCAGAAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722089.725386/100‑1 (MQ=60)
                AAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAG                                                                                      >  SRR3722089.1149236/1‑100 (MQ=60)
                    GGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGG                                                                                  <  SRR3722089.1019734/100‑1 (MQ=60)
                     GCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGT                                                                                 <  SRR3722089.894513/100‑1 (MQ=60)
                       ACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGA                                                                               >  SRR3722089.906633/1‑100 (MQ=60)
                        CCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGAC                                                                              <  SRR3722089.314782/100‑1 (MQ=60)
                          TTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                            <  SRR3722089.669670/100‑1 (MQ=60)
                          TTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                            <  SRR3722089.619201/100‑1 (MQ=60)
                                    gcTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATT                                                                  >  SRR3722089.1184203/3‑100 (MQ=60)
                                        TTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACG                                                              <  SRR3722089.44552/100‑1 (MQ=60)
                                             ATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAC                                                         >  SRR3722089.743350/1‑100 (MQ=60)
                                              TTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACA                                                        >  SRR3722089.237912/1‑100 (MQ=60)
                                               TCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAT                                                       >  SRR3722089.1001537/1‑100 (MQ=60)
                                               TCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAT                                                       <  SRR3722089.1188577/100‑1 (MQ=60)
                                                   TATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCT                                                   <  SRR3722089.106510/100‑1 (MQ=60)
                                                                        GAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC                              >  SRR3722089.1186046/1‑100 (MQ=60)
                                                                        GAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC                              >  SRR3722089.1175287/1‑100 (MQ=60)
                                                                        GAATAATCCCGGTATTGTCGGCATAAAGGGGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC                              >  SRR3722089.765868/1‑100 (MQ=60)
                                                                              TCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACC                        >  SRR3722089.820766/1‑100 (MQ=60)
                                                                                            GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722089.493809/1‑100 (MQ=60)
                                                                                            GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722089.778019/1‑100 (MQ=60)
                                                                                            GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722089.803155/1‑100 (MQ=60)
                                                                                            GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          <  SRR3722089.220664/100‑1 (MQ=60)
                                                                                                    GTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATCGCC  <  SRR3722089.471584/100‑1 (MQ=60)
                                                                                                    |                                                                                                   
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATCGCC  >  NZ_CP009273/4108731‑4108929

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