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 I211 R1 223 18.0 1007644 96.3% 970361 85.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,751,237 A→G V241V (GTT→GTC BW25113_RS13685 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,751,2370AG100.0% 73.9 / NA 23V241V (GTT→GTCBW25113_RS13685hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/13);  total (10/13)

CAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTT  >  NZ_CP009273/2751152‑2751315
                                                                                     |                                                                              
cAAACTCTGACATTTGTATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACag                                                                            <  1:356890/90‑1 (MQ=255)
  aaCTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACagag                                                                          >  1:476031/1‑90 (MQ=255)
      ctGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCat                                                                      >  2:43254/1‑90 (MQ=255)
      ctGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCat                                                                      >  1:336832/1‑90 (MQ=255)
       tGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCatt                                                                     >  2:51331/1‑90 (MQ=255)
                   aaaaTCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCACTCGGGACAGAGTCATTAATTAGCTcaac                                                         <  2:122486/90‑1 (MQ=255)
                             aGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGt                                               <  1:460153/90‑1 (MQ=255)
                                ggTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGa                                            >  2:464124/1‑90 (MQ=255)
                                 gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACag                                                                            <  1:161520/57‑1 (MQ=255)
                                 gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACag                                                                            >  2:161520/1‑57 (MQ=255)
                                  ttCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAAt                                          <  1:424266/90‑1 (MQ=255)
                                       tAAGTTCATTAACTGCTACAACCTTATTTCTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGatat                                     <  2:272301/90‑1 (MQ=255)
                                                   cTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCg                         >  1:90921/1‑90 (MQ=255)
                                                    tGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgc                        >  1:312507/1‑90 (MQ=255)
                                                     gCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgct                       >  1:302707/1‑90 (MQ=255)
                                                      cTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctt                      <  2:160074/90‑1 (MQ=255)
                                                          aaCCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctttgct                  <  2:300479/90‑1 (MQ=255)
                                                                  ttttCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCAc          <  1:264289/90‑1 (MQ=255)
                                                                    ttCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCAccc        <  2:339778/90‑1 (MQ=255)
                                                                     tCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCct       <  2:336832/90‑1 (MQ=255)
                                                                       cACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCcttc     <  1:51331/90‑1 (MQ=255)
                                                                       cACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCcttc     <  1:464124/90‑1 (MQ=255)
                                                                          ccATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCttt  >  1:196214/1‑90 (MQ=255)
                                                                                     |                                                                              
CAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTT  >  NZ_CP009273/2751152‑2751315

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTG  >  NZ_CP009273/2751140‑2751330
                                                                                                 |                                                                                             
CATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGAC                                                                                             <  SRR3722094.66763/100‑1 (MQ=60)
    ATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAG                                                                                         >  SRR3722094.484442/1‑100 (MQ=60)
        ATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCAT                                                                                     >  SRR3722094.342361/1‑100 (MQ=60)
            cgtcagatgtgtataagagacagTCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAAT                                                                                 <  SRR3722094.163773/77‑1 (MQ=60)
            CAAACTCTGACATTTGTATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAAT                                                                                 <  SRR3722094.362862/100‑1 (MQ=60)
                                         AGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATAT                                                    <  SRR3722094.468253/100‑1 (MQ=60)
                                              TTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTAC                                               <  SRR3722094.431672/100‑1 (MQ=60)
                                                     AGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCG                                        >  SRR3722094.92242/1‑100 (MQ=60)
                                                      GTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGC                                       >  SRR3722094.317422/1‑100 (MQ=60)
                                                       TTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCT                                      >  SRR3722094.307383/1‑100 (MQ=60)
                                                                            ATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTT                 >  SRR3722094.198898/1‑100 (MQ=60)
                                                                              TTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCC               <  SRR3722094.268219/100‑1 (MQ=60)
                                                                                   CACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCAC          <  SRR3722094.472293/100‑1 (MQ=60)
                                                                                   CACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCAC          <  SRR3722094.52041/100‑1 (MQ=60)
                                                                                           ATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTG  >  SRR3722094.296475/1‑100 (MQ=60)
                                                                                                 |                                                                                             
GATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTG  >  NZ_CP009273/2751140‑2751330

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