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 I210 R1 226 19.6 1130256 95.4% 1078264 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,670,359 T→C G454G (GGA→GGG yphE ← sugar ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,670,3590TC100.0% 85.9 / NA 27G454G (GGA→GGGyphEsugar ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/7);  total (20/7)

CGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGC  >  NZ_CP009273/2670272‑2670440
                                                                                       |                                                                                 
cGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTccc                                                                                 <  1:26731/90‑1 (MQ=255)
            gtgcTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagccag                                                                     >  2:62653/1‑90 (MQ=255)
                tGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCa                                                                 <  1:19586/90‑1 (MQ=255)
                      aaCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACa                                                           >  2:43162/1‑90 (MQ=255)
                        cAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACCCCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                         <  2:276715/90‑1 (MQ=255)
                           gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACg                                                               <  2:271378/81‑1 (MQ=255)
                           gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACg                                                               >  1:271378/1‑81 (MQ=255)
                           gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                      >  2:414216/1‑90 (MQ=255)
                           gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                      >  2:200399/1‑90 (MQ=255)
                              gCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTa                                                   >  1:249846/1‑90 (MQ=255)
                                ggTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTaaa                                                 >  2:319702/1‑90 (MQ=255)
                                         ccAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGtt                                        >  1:188430/1‑90 (MQ=255)
                                           aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                      >  2:516070/1‑90 (MQ=255)
                                                      ctccACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACgg                                                              <  1:552725/55‑1 (MQ=255)
                                                      ctccACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACgg                                                              >  2:552725/1‑55 (MQ=255)
                                                          aCTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGAcg                       >  1:318984/1‑90 (MQ=255)
                                                           cTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGAcgc                      >  2:320156/1‑90 (MQ=255)
                                                             tCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGAcgccg                    <  1:86215/90‑1 (MQ=255)
                                                               aCTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAACACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCcgcg                  >  2:295066/1‑90 (MQ=255)
                                                                     gaTAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTgg            >  2:94268/1‑90 (MQ=255)
                                                                       tAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCt          >  1:122433/1‑90 (MQ=255)
                                                                          acacCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgt       >  2:134299/1‑90 (MQ=255)
                                                                            accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCtt                                  <  2:263619/61‑1 (MQ=255)
                                                                            accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCtt                                  >  1:263619/1‑61 (MQ=255)
                                                                            accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgtcg     >  2:275438/1‑90 (MQ=255)
                                                                             ccacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGa    >  2:496843/1‑90 (MQ=255)
                                                                               acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc  >  1:256317/1‑90 (MQ=255)
                                                                               acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc  <  1:200399/90‑1 (MQ=255)
                                                                                       |                                                                                 
CGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGC  >  NZ_CP009273/2670272‑2670440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ 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
ACTCCTGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTG  >  NZ_CP009273/2670265‑2670451
                                                                                              |                                                                                            
ACTCCTGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTC                                                                                         <  SRR3722092.399504/100‑1 (MQ=60)
       CGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCC                                                                                  <  SRR3722092.27192/100‑1 (MQ=60)
                       TGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATAC                                                                  <  SRR3722092.19924/100‑1 (MQ=60)
                        GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                 >  SRR3722092.275937/1‑100 (MQ=60)
                          agcgtcagatgtgtataagagacagGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGT                                                               <  SRR3722092.563819/75‑1 (MQ=60)
                           ATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTA                                                              >  SRR3722092.253892/1‑100 (MQ=60)
                                      CGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTT                                                   >  SRR3722092.191563/1‑100 (MQ=60)
                                                       CAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACG                                  >  SRR3722092.324624/1‑100 (MQ=60)
                                                                    TCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCT                     >  SRR3722092.124509/1‑100 (MQ=60)
                                                                    TCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCT                     <  SRR3722092.87692/100‑1 (MQ=60)
                                                                         GGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACctgtctcttatacacatct                >  SRR3722092.267985/1‑81 (MQ=60)
                                                                            GATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGC             >  SRR3722092.260496/1‑100 (MQ=60)
                                                                                      ACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCT   <  SRR3722092.203712/100‑1 (MQ=60)
                                                                                       CGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTG  >  SRR3722092.108710/1‑100 (MQ=60)
                                                                                              |                                                                                            
ACTCCTGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTG  >  NZ_CP009273/2670265‑2670451

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