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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,189,312 A→G T512A (ACC→GCC)  metG → methionine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,189,3120AG100.0% 83.4 / NA 26T512A (ACC→GCC) metGmethionine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/15);  total (11/15)

CTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAG  >  NZ_CP009273/2189228‑2189395
                                                                                    |                                                                                   
cTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTgg                                                                                >  2:380287/1‑90 (MQ=255)
  gTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGa                                                                              >  1:150908/1‑90 (MQ=255)
     ccGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATgg                                                                           >  2:43047/1‑90 (MQ=255)
       gcgTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTa                                                                         >  2:253278/1‑90 (MQ=255)
                 gaCTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAAcc                                                               <  2:505635/90‑1 (MQ=255)
                     tACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCgctg                                                           >  1:310149/1‑90 (MQ=255)
                              ccGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCcac                                                  <  2:324905/90‑1 (MQ=255)
                              ccGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCcac                                                  <  1:576105/90‑1 (MQ=255)
                              ccGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCcac                                                  <  1:688591/90‑1 (MQ=255)
                                ggTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACaa                                                >  2:195235/1‑90 (MQ=255)
                                        cGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCGCAAAGTGAATc                                        <  1:390965/90‑1 (MQ=255)
                                           aaCTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGt                                     >  2:195236/1‑90 (MQ=255)
                                           aaCTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGt                                     >  1:561135/1‑90 (MQ=255)
                                           aaCTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGAAACCGCTGCTGGGCCACAAAGTGAATCCGt                                     >  2:693673/1‑90 (MQ=255)
                                                aCCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGt                                                                          >  1:399369/1‑48 (MQ=255)
                                                aCCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGt                                                                          <  2:399369/48‑1 (MQ=255)
                                                  cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                              <  1:309982/90‑1 (MQ=255)
                                                  cGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgc                              <  1:270323/90‑1 (MQ=255)
                                                    aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                            <  1:693673/90‑1 (MQ=255)
                                                    aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                            <  1:43047/90‑1 (MQ=255)
                                                    aGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGgcgc                            <  1:383524/90‑1 (MQ=255)
                                                                 aTTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCAt               >  2:275118/1‑90 (MQ=255)
                                                                  ttCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATc              <  1:625917/90‑1 (MQ=255)
                                                                  ttCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATc              <  1:380287/90‑1 (MQ=255)
                                                                            cGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGc    <  1:253278/90‑1 (MQ=255)
                                                                              gAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAg  <  2:561135/90‑1 (MQ=255)
                                                                                    |                                                                                   
CTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAG  >  NZ_CP009273/2189228‑2189395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGAAGC  >  NZ_CP009273/2189220‑2189403
                                                                                            |                                                                                           
GCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGA                                                                                      >  SRR3722112.152769/1‑100 (MQ=60)
                   GCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTG                                                                   >  SRR3722112.313706/1‑100 (MQ=60)
                                      CCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATC                                                <  SRR3722112.584895/100‑1 (MQ=60)
                                      CCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATC                                                <  SRR3722112.699345/100‑1 (MQ=60)
                                         GTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGT                                             >  SRR3722112.569643/1‑100 (MQ=60)
                                              GCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCtgtctcttatacacatctgacgctgccgacg                                        >  SRR3722112.404624/1‑69 (MQ=60)
                                                CGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCGCAAAGTGAATCCGTTCAAGGC                                      <  SRR3722112.396083/100‑1 (MQ=60)
                                                          CGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAAC                            <  SRR3722112.273485/100‑1 (MQ=60)
                                                          CGAGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAAC                            <  SRR3722112.313538/100‑1 (MQ=60)
                                                            AGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCG                          <  SRR3722112.388507/100‑1 (MQ=60)
                                                            AGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCG                          <  SRR3722112.43567/100‑1 (MQ=60)
                                                            AGCGTGCAGAAGCATTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCG                          <  SRR3722112.704508/100‑1 (MQ=60)
                                                                          TTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGC            <  SRR3722112.385224/100‑1 (MQ=60)
                                                                          TTCCTCAATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGC            <  SRR3722112.635579/100‑1 (MQ=60)
                                                                                 ATACGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGA     >  SRR3722112.460221/1‑100 (MQ=60)
                                                                                    CGGAACTGGCCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGAAGC  <  SRR3722112.256267/100‑1 (MQ=60)
                                                                                            |                                                                                           
GCATCAACCTGTTCCGCGTGCTGATGACTTACCTGAAGCCGGTACTGCCGAAACTGACCGAGCGTGCAGAAGCATTCCTCAATACGGAACTGACCTGGGATGGTATCCAGCAACCGCTGCTGGGCCACAAAGTGAATCCGTTCAAGGCGCTGTATAACCGCATCGATATGAGGCAGGTTGAAGC  >  NZ_CP009273/2189220‑2189403

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