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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,276,251 Δ1 bp coding (932/3744 nt) narG → nitrate reductase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,276,2490G.100.0% 80.8 / NA 20coding (930/3744 nt)narGnitrate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (10/10);  total (10/10)

CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGAC  >  NZ_CP009273/1276162‑1276327
                                                                                       |                                                                              
cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                              <  2:834901/90‑2 (MQ=255)
cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                              <  1:400946/90‑2 (MQ=255)
cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                              <  1:970030/90‑2 (MQ=255)
       gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTaa                                                                      >  2:1345479/1‑90 (MQ=255)
           aCGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAAtgct                                                                  <  1:1345479/90‑1 (MQ=255)
              cTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCg                                                               >  1:993800/1‑90 (MQ=255)
                     cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTc                                                        >  2:1244029/1‑90 (MQ=255)
                       ccAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCa                                                      <  2:1081606/90‑1 (MQ=255)
                             tgtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGa                                                <  2:759765/90‑1 (MQ=255)
                              gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                              <  1:631181/60‑2 (MQ=255)
                              gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                              >  2:631181/1‑59 (MQ=255)
                              gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                               >  2:1219698/1‑90 (MQ=255)
                                      tgtgGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                                       >  1:516315/1‑90 (MQ=255)
                                      tgtgGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAg                                       >  1:590232/1‑90 (MQ=255)
                                        tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                     <  1:1219698/90‑1 (MQ=255)
                                        tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                     >  2:579971/1‑90 (MQ=255)
                                        tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                     <  1:558171/90‑1 (MQ=255)
                                            tggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTa                                 <  1:887651/90‑1 (MQ=255)
                                             ggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                                >  1:9882/1‑90 (MQ=255)
                                                aCCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTc                             <  2:1533077/90‑1 (MQ=255)
                                                  cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc                           <  1:103535/90‑1 (MQ=255)
                                                   gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAcc                          >  2:188262/1‑90 (MQ=255)
                                                                           ggCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCCCCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGAc  <  2:1371053/90‑1 (MQ=255)
                                                                              gcTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac        >  1:91894/1‑81 (MQ=255)
                                                                              gcTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac        <  2:91894/81‑1 (MQ=255)
                                                                                       |                                                                              
CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGAC  >  NZ_CP009273/1276162‑1276327

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
CCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGA  >  NZ_CP009273/1276158‑1276350
                                                                                             |                                                                                                   
CCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTA                                                                                              <  SRR3722078.1289353/100‑1 (MQ=60)
    cagatgtgtataagagacagTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGC                                                                                          <  SRR3722078.638443/80‑1 (MQ=60)
    CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGC                                                                                          <  SRR3722078.405810/100‑1 (MQ=60)
    CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGC                                                                                          <  SRR3722078.983893/100‑1 (MQ=60)
        CCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCG                                                                                      >  SRR3722078.1008251/1‑100 (MQ=60)
               ACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTC                                                                               <  SRR3722078.1367074/100‑1 (MQ=60)
                                CTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAG                                                              >  SRR3722078.522424/1‑100 (MQ=60)
                                CTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAG                                                              >  SRR3722078.597092/1‑100 (MQ=60)
                                       ATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAT                                                       >  SRR3722078.9988/1‑100 (MQ=60)
                                            TGGCTGGCACCGAAACAGGGCCCCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAC                                                  <  SRR3722078.1238895/100‑1 (MQ=60)
                                            TGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAC                                                  <  SRR3722078.564745/100‑1 (MQ=60)
                                                TGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGAC                                              <  SRR3722078.899579/100‑1 (MQ=60)
                                                      CGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTG                                        <  SRR3722078.104765/100‑1 (MQ=60)
                                                                  CCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCctgtctcttatacacatctgacgctgccgacgatggcggtcgtgtag                            >  SRR3722078.535866/1‑53 (MQ=60)
                                                                        CGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATG                      >  SRR3722078.92974/1‑100 (MQ=60)
                                                                                GCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCT              >  SRR3722078.5476/1‑100 (MQ=60)
                                                                                             GCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCtgtctcttatacacatctgacgctgccgacgatggcg  >  SRR3722078.1274951/1‑63 (MQ=60)
                                                                                             |                                                                                                   
CCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGA  >  NZ_CP009273/1276158‑1276350

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