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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

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,7770TC100.0% 89.1 / NA 27I137V (ATC→GTC) mntHNramp family divalent metal transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCG  >  NZ_CP009273/2505690‑2505854
                                                                                       |                                                                             
gggCAACAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:83141/89‑1 (MQ=255)
cgGCAACAACCAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1182727/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  2:444298/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             >  2:216146/1‑90 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:855002/90‑1 (MQ=255)
cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:487853/90‑1 (MQ=255)
cgGCAAAAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1085649/90‑1 (MQ=255)
cgGAAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  2:675079/90‑1 (MQ=255)
     acaaacaaCAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTc                                                                        <  1:322071/90‑1 (MQ=255)
          caacaGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAc                                                                   <  1:925200/90‑1 (MQ=255)
          caacaGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAc                                                                   <  2:150315/90‑1 (MQ=255)
             caGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACcgc                                                                <  2:136747/90‑1 (MQ=255)
               gTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              <  1:123783/90‑1 (MQ=255)
               gTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              <  1:257849/90‑1 (MQ=255)
               gTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              >  1:1172311/1‑90 (MQ=255)
                tAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGcc                                                             >  1:329110/1‑90 (MQ=255)
                 aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGccc                                                            <  2:1038660/90‑1 (MQ=255)
                 aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGccc                                                            <  1:216146/90‑1 (MQ=255)
                  aCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             >  2:1011834/1‑72 (MQ=255)
                  aCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1011834/72‑1 (MQ=255)
                  aCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCt                                                           >  2:609279/1‑90 (MQ=255)
                    ccGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGc                                                         >  2:978931/1‑90 (MQ=255)
                    ccGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGc                                                         >  1:560228/1‑90 (MQ=255)
                    ccGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGc                                                         >  2:147215/1‑90 (MQ=255)
                      gCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              >  1:95367/1‑83 (MQ=255)
                      gCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              <  2:95367/83‑1 (MQ=255)
                      gCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGcaa                                                       >  1:401272/1‑90 (MQ=255)
                             aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1215057/61‑1 (MQ=255)
                             aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             >  2:1215057/1‑61 (MQ=255)
                             aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGaa                                                >  1:724657/1‑90 (MQ=255)
                                         ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                    <  1:204357/90‑1 (MQ=255)
                                         ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                    <  2:789401/90‑1 (MQ=255)
                                         ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                    <  1:837209/90‑1 (MQ=255)
                                         ggTTTTCGCCCGCGACGTTGCAGCCTTAAAATCAGGAAAGTAGAGACCCCCGTCAGCCCCGCGCCCTGCCACAACGAAACACCAAGAATg                                    <  1:701100/90‑1 (MQ=255)
                                                    gcgACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAcc                         >  2:999059/1‑90 (MQ=255)
                                                                 aTTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCa            >  1:881130/1‑90 (MQ=255)
                                                                   tAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAAt          <  2:554207/90‑1 (MQ=255)
                                                                           ggAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCg  >  1:834427/1‑90 (MQ=255)
                                                                           ggAAAGTCGCGACCCCCGCCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCg  >  2:698632/1‑90 (MQ=255)
                                                                                       |                                                                             
CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCG  >  NZ_CP009273/2505690‑2505854

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ 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
AGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCA  >  NZ_CP009273/2505685‑2505869
                                                                                            |                                                                                            
AGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTC                                                                                       <  SRR3722089.994596/100‑1 (MQ=60)
     agatgtgtataagagacagCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.1233543/81‑1 (MQ=60)
     agagacagAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.1026905/92‑1 (MQ=60)
     GGGCAACAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.84155/100‑1 (MQ=60)
     CGGCAACAACCAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.1200737/100‑1 (MQ=60)
     CGGCAACAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.867196/100‑1 (MQ=60)
     CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.493644/100‑1 (MQ=60)
     CGGCAAAAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAC                                                                                  <  SRR3722089.1102003/100‑1 (MQ=60)
          ACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGC                                                                             <  SRR3722089.326130/100‑1 (MQ=60)
          ACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGC                                                                             >  SRR3722089.1190155/1‑100 (MQ=60)
           CAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCC                                                                            >  SRR3722089.333235/1‑100 (MQ=60)
               CAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGC                                                                        <  SRR3722089.938700/100‑1 (MQ=60)
               CAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGC                                                                        >  SRR3722089.566761/1‑100 (MQ=60)
                 ACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAA                                                                      >  SRR3722089.406143/1‑100 (MQ=60)
                 ACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAA                                                                      >  SRR3722089.96552/1‑100 (MQ=60)
                    GTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAA                                                                   <  SRR3722089.261162/100‑1 (MQ=60)
                    GTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAA                                                                   <  SRR3722089.125294/100‑1 (MQ=60)
                      AACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACG                                                                 <  SRR3722089.218882/100‑1 (MQ=60)
                        CCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAA                                                               >  SRR3722089.734090/1‑100 (MQ=60)
                                              GGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAC                                         <  SRR3722089.849056/100‑1 (MQ=60)
                                              GGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAC                                         <  SRR3722089.206940/100‑1 (MQ=60)
                                              GGTTTTCGCCCGCGACGTTGCAGCCTTAAAATCAGGAAAGTAGAGACCCCCGTCAGCCCCGCGCCCTGCCACAACGAAACACCAAGAATGAGTTTAAAAC                                         <  SRR3722089.710240/100‑1 (MQ=60)
                                                            ACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCA                           >  SRR3722089.893823/1‑100 (MQ=60)
                                                                      ATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCG                 >  SRR3722089.846222/1‑100 (MQ=60)
                                                                                GGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGT       >  SRR3722089.516700/1‑100 (MQ=60)
                                                                                     GTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCA  >  SRR3722089.1117738/1‑100 (MQ=60)
                                                                                            |                                                                                            
AGCCGCGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCGCCAGGTCGGTTGCCA  >  NZ_CP009273/2505685‑2505869

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