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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,505,777 T→C I137V (ATC→GTC)  mntH ← Nramp family divalent metal transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,505,7770TC90.0% 44.2 / 0.1 20I137V (ATC→GTC) mntHNramp family divalent metal transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (6/12);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGG  >  NZ_CP009273/2505690‑2505859
                                                                                       |                                                                                  
gggCAAAAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  <  2:701133/89‑1 (MQ=255)
cgGCAACAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  <  2:339573/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  <  1:323865/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  <  1:578441/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  >  1:604655/1‑90 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  <  1:693602/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  >  2:319682/1‑90 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  >  2:251817/1‑90 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATcc                                                                                  <  1:810299/90‑1 (MQ=255)
     acaaacaaCAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTc                                                                             <  2:663036/90‑1 (MQ=255)
          caacaGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAc                                                                        >  2:657835/1‑90 (MQ=255)
             caGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACcgc                                                                     <  2:646936/90‑1 (MQ=255)
                 aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  >  2:795179/1‑73 (MQ=255)
                 aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                                  <  1:795179/73‑1 (MQ=255)
                         aaTCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGcaacaa                                                         >  1:779096/1‑90 (MQ=255)
                           tCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACg                                                       <  1:748705/90‑1 (MQ=255)
                             aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGaa                                                     <  2:193467/90‑1 (MQ=255)
                             aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGaa                                                     <  2:382146/90‑1 (MQ=255)
                                     cAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAg                                             >  2:98701/1‑90 (MQ=255)
                                         ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                         <  1:98701/90‑1 (MQ=255)
                                         ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                         <  1:251817/90‑1 (MQ=255)
                                         ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                         <  1:238226/90‑1 (MQ=255)
                                                           tGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATcgcc                       <  2:410497/90‑1 (MQ=255)
                                                             cagcaTTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATcgccgc                     <  2:3929/90‑1 (MQ=255)
                                                                 aTTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCa                 >  2:572412/1‑90 (MQ=255)
                                                                 aTTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCa                 >  2:210041/1‑90 (MQ=255)
                                                                   tAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAAt               <  1:466636/90‑1 (MQ=255)
                                                                       aTCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAAt           <  2:791894/90‑1 (MQ=255)
                                                                         cAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTc         >  1:314919/1‑90 (MQ=255)
                                                                                gTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAgg  >  2:820768/1‑90 (MQ=255)
                                                                                       |                                                                                  
CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGG  >  NZ_CP009273/2505690‑2505859

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ 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
CAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCC  >  NZ_CP009273/2505678‑2505856
                                                                                                   |                                                                               
CAATGTAAGCCGCGGAAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAC                                                                                 <  SRR3722079.759500/100‑1 (MQ=60)
  ATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCC                                                                               >  SRR3722079.610753/1‑100 (MQ=60)
       AGCCGCGGCAAAAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTC                                                                          <  SRR3722079.777326/100‑1 (MQ=60)
            gagacagAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                     <  SRR3722079.803711/93‑1 (MQ=60)
            CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCAC                                                                     <  SRR3722079.818991/100‑1 (MQ=60)
            CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                     <  SRR3722079.326464/100‑1 (MQ=60)
            CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                     <  SRR3722079.584165/100‑1 (MQ=60)
            CGGCAAAAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                     <  SRR3722079.700902/100‑1 (MQ=60)
                           GTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAA                                                      >  SRR3722079.787445/1‑100 (MQ=60)
                                       TCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAG                                          <  SRR3722079.756705/100‑1 (MQ=60)
                                                     GGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAC                            <  SRR3722079.240232/100‑1 (MQ=60)
                                                     GGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAC                            <  SRR3722079.253927/100‑1 (MQ=60)
                                                     GGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAC                            <  SRR3722079.99481/100‑1 (MQ=60)
                                                                           GCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTC      >  SRR3722079.317464/1‑100 (MQ=60)
                                                                               TAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCC  <  SRR3722079.470716/100‑1 (MQ=60)
                                                                                                   |                                                                               
CAATGTAAGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCC  >  NZ_CP009273/2505678‑2505856

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