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 F1 I1 R1 315 33.5 1889570 98.7% 1865005 80.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,224,682 T→G D285E (GAT→GAG)  metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,224,6820TG100.0% 149.6 / NA 42D285E (GAT→GAG) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (20/22);  total (20/22)

GGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCT  >  NC_000913/4224604‑4224758
                                                                              |                                                                            
ggAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAgct                                                                            <  1:839940/81‑1 (MQ=255)
 gAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAgctc                                                                           <  1:839941/81‑1 (MQ=255)
 gAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAgctc                                                                           <  1:839942/81‑1 (MQ=255)
  aGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCg                                                                          <  2:839943/81‑1 (MQ=255)
    cTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGAc                                                                        <  1:839944/81‑1 (MQ=255)
         aCGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCg                                                                    <  1:839945/80‑1 (MQ=255)
          cGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGac                                                                  <  2:839946/81‑1 (MQ=255)
          cGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGa                                                                   >  1:840015/1‑80 (MQ=255)
            gATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGacac                                                                <  1:839947/81‑1 (MQ=255)
                gCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGAt                                                             <  1:839948/80‑1 (MQ=255)
                 cGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATgg                                                           >  1:840004/1‑81 (MQ=255)
                  ggAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATgg                                                           >  2:840017/1‑80 (MQ=255)
                  ggAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGc                                                          <  1:839949/81‑1 (MQ=255)
                    aaTGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCa                                                         >  1:840009/1‑80 (MQ=255)
                     aTGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCaaa                                                       <  2:839950/81‑1 (MQ=255)
                     aTGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCaa                                                        >  1:839991/1‑80 (MQ=255)
                              cACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACg                                              >  2:840028/1‑81 (MQ=255)
                                    gcACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                        >  1:840007/1‑81 (MQ=255)
                                    gcACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                        >  1:840008/1‑81 (MQ=255)
                                    gcACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                        <  1:839951/81‑1 (MQ=255)
                                     cACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATgg                                       <  2:839952/81‑1 (MQ=255)
                                          gAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGgcgc                                   >  2:840012/1‑80 (MQ=255)
                                          gAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCa                                  <  1:839953/81‑1 (MQ=255)
                                              gCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCg                              >  2:840002/1‑81 (MQ=255)
                                              gCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCg                              <  1:839954/81‑1 (MQ=255)
                                                 gggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGt                           >  1:840020/1‑81 (MQ=255)
                                                  ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          >  1:840001/1‑81 (MQ=255)
                                                  ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          >  1:840024/1‑81 (MQ=255)
                                                  ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          <  2:839955/81‑1 (MQ=255)
                                                  ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          <  2:839956/81‑1 (MQ=255)
                                                   gCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGttt                         >  2:840016/1‑81 (MQ=255)
                                                    cTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtttt                        <  2:839957/81‑1 (MQ=255)
                                                       cccAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtctc                     <  2:839958/81‑1 (MQ=255)
                                                          aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCaa                   >  2:840032/1‑80 (MQ=255)
                                                          aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                  <  2:839959/81‑1 (MQ=255)
                                                          aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                  <  2:839960/81‑1 (MQ=255)
                                                          aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                  <  2:839961/81‑1 (MQ=255)
                                                            cGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAatat                >  2:840005/1‑81 (MQ=255)
                                                             gCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtc               <  1:839962/81‑1 (MQ=255)
                                                               cTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgt             >  1:840033/1‑81 (MQ=255)
                                                                   ggTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggc         >  1:840031/1‑81 (MQ=255)
                                                                       aGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgct     <  1:839963/81‑1 (MQ=255)
                                                                        gTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGTAAGCGGGTTTTCTCAATATCGTCGGCGgctg    >  2:840025/1‑81 (MQ=255)
                                                                         tACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctgc   <  1:839964/81‑1 (MQ=255)
                                                                          aCGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctgct  >  1:840022/1‑81 (MQ=255)
                                                                              |                                                                            
GGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCT  >  NC_000913/4224604‑4224758

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A