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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,965,782 T→G E389D (GAA→GAC tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,965,7820TG100.0% 93.8 / NA 28E389D (GAA→GACtarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/15);  total (13/15)

AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTC  >  NZ_CP009273/1965696‑1965868
                                                                                      |                                                                                      
aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                                     >  2:580959/1‑90 (MQ=255)
           ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                          >  2:918494/1‑90 (MQ=255)
           ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                          >  2:546573/1‑90 (MQ=255)
              ttttGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGc                                                                       >  1:531601/1‑90 (MQ=255)
                    cgccGGGGCGCTGCGACTGGCAAGATTCCGCCCTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAc                                                                 <  2:305678/90‑1 (MQ=25)
                       cTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACcgc                                                              >  1:256865/1‑90 (MQ=255)
                                gCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAg                                                     >  1:795118/1‑90 (MQ=255)
                                  gccTGGCAAGATTCCGCTCTTCTCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcg                                                   <  1:201025/88‑1 (MQ=21)
                                   aCTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcgc                                                  >  1:723110/1‑90 (MQ=255)
                                    cTGGCAAGATTCCCCCCTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                                 <  2:119586/90‑1 (MQ=255)
                                    cTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                                 <  1:498820/90‑1 (MQ=255)
                                              tACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAg                                       <  1:124254/90‑1 (MQ=255)
                                                   aCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                                  <  1:723266/90‑1 (MQ=255)
                                                     ttCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGaa                                >  2:468732/1‑90 (MQ=255)
                                                      tCACCCGCCACCACGGCAAAACCACGGGCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                               >  1:480469/1‑90 (MQ=255)
                                                      tCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg                                               >  1:806149/1‑74 (MQ=255)
                                                      tCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg                                               <  2:806149/74‑1 (MQ=255)
                                                          ccGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCaa                           <  2:531601/90‑1 (MQ=255)
                                                              caccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                               <  2:642101/82‑1 (MQ=255)
                                                              caccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                               >  1:642101/1‑82 (MQ=255)
                                                              caccacGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATAcc                       >  1:50182/1‑90 (MQ=255)
                                                               accacGGCAAAACCACGGCCCTGGTCCCCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCg                      <  1:920932/90‑1 (MQ=255)
                                                                     gCAAAACCACGGCCCTGGTCGCCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa                <  1:23785/90‑1 (MQ=255)
                                                                     gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa                <  1:468732/90‑1 (MQ=255)
                                                                       aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc              <  1:424066/90‑1 (MQ=255)
                                                                       aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc              <  1:92515/90‑1 (MQ=255)
                                                                       aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc              <  1:374473/90‑1 (MQ=255)
                                                                                   cTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTc  >  1:897680/1‑90 (MQ=255)
                                                                                      |                                                                                      
AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTC  >  NZ_CP009273/1965696‑1965868

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 26 ≤ 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
GAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTT  >  NZ_CP009273/1965684‑1965880
                                                                                                  |                                                                                                  
GAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGT                                                                                                   >  SRR3722072.686043/1‑100 (MQ=60)
 AGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTC                                                                                                  <  SRR3722072.299757/100‑1 (MQ=60)
                CTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGC                                                                                   >  SRR3722072.541185/1‑100 (MQ=60)
                         CTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGC                                                                          >  SRR3722072.261220/1‑100 (MQ=60)
                                  CCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAG                                                                 >  SRR3722072.809020/1‑100 (MQ=60)
                                     GGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGC                                                              >  SRR3722072.735966/1‑100 (MQ=60)
                                             CGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCctgtctcttatac                                                      >  SRR3722072.545260/1‑87 (MQ=60)
                                              GCCTGGCAAGATTCCGCTCTTCTCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATT                                                     <  SRR3722072.204483/100‑1 (MQ=60)
                                                CTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAG                                                   <  SRR3722072.507836/100‑1 (MQ=60)
                                                        ATTACGCACTTCACCCGCCACCACGGCAAAACCACGGGCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAG                                           >  SRR3722072.489186/1‑100 (MQ=60)
                                                        ATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCctgtct                                           >  SRR3722072.820208/1‑94 (MQ=60)
                                                          TACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGC                                         <  SRR3722072.126517/100‑1 (MQ=60)
                                                               ACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATAC                                    <  SRR3722072.736122/100‑1 (MQ=60)
                                                                CTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACC                                   >  SRR3722072.51101/1‑100 (MQ=60)
                                                                CTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACC                                   >  SRR3722072.653553/1‑100 (MQ=60)
                                                                           ACCACGGCAAAACCACGGCCCTGGTCCCCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGC                        <  SRR3722072.936736/100‑1 (MQ=60)
                                                                                 GCAAAACCACGGCCCTGGTCGCCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA                  <  SRR3722072.24206/100‑1 (MQ=60)
                                                                                 GCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA                  <  SRR3722072.477278/100‑1 (MQ=60)
                                                                                   AAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATG                <  SRR3722072.431779/100‑1 (MQ=60)
                                                                                   AAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATG                <  SRR3722072.381145/100‑1 (MQ=60)
                                                                                   AAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATG                <  SRR3722072.94217/100‑1 (MQ=60)
                                                                                     AACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTC              >  SRR3722072.913129/1‑100 (MQ=60)
                                                                                                 GGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTT  >  SRR3722072.744431/1‑100 (MQ=60)
                                                                                                  |                                                                                                  
GAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTT  >  NZ_CP009273/1965684‑1965880

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