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 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 987,566 A→T I497F (ATC→TTC)  pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273987,5660AT100.0% 81.4 / NA 24I497F (ATC→TTC) pepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (11/13);  total (11/13)

ACGCCAGCCACGCCGGATCAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGATCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTA  >  NZ_CP009273/987482‑987651
                                                                                    |                                                                                     
aCGCCAGCCACGCCGGATCAGGCAGAAAAACAGCCGCTGCATATTACGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGt                                                                                 >  1:668936/1‑91 (MQ=255)
           gCCGGATCAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGc                                                                      >  2:511/1‑91 (MQ=255)
               gATCAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTc                                                                  >  2:444855/1‑91 (MQ=255)
                  cAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATc                                                               <  2:83893/91‑1 (MQ=255)
                  cAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATc                                                               <  1:563633/91‑1 (MQ=255)
                        gAAAAACAGCCGCTGCATGTTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGa                                                         >  1:545890/1‑91 (MQ=255)
                              cAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCg                                                   <  1:399262/91‑1 (MQ=255)
                                      gCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAAc                                           <  2:37639/91‑1 (MQ=255)
                                                 ttGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGTAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGc                                <  1:545232/91‑1 (MQ=255)
                                                 ttGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGc                                <  2:133934/91‑1 (MQ=255)
                                                 ttGCCAGCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGc                                <  2:143409/91‑1 (MQ=255)
                                                     cATCGAACTGTATGATAACGAAGGCAAAGCGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGACCGTCACTCAGGCGGaa                            >  2:651516/1‑91 (MQ=255)
                                                             tGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCtt                    <  1:274214/91‑1 (MQ=255)
                                                             tGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCtt                    <  1:821709/91‑1 (MQ=255)
                                                                 tGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTc                >  1:450422/1‑91 (MQ=255)
                                                                    tAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTACGTGCTGACAGTCACTCAGGCGGAACAGACCTTTGTCttt             >  2:650534/1‑91 (MQ=255)
                                                                        gAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATa         >  2:172326/1‑91 (MQ=255)
                                                                        gAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATa         >  2:127293/1‑91 (MQ=255)
                                                                         aaGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATaa        <  2:483823/91‑1 (MQ=255)
                                                                          aGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTc                                        >  2:32065/1‑58 (MQ=255)
                                                                          aGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTc                                        <  1:32065/58‑1 (MQ=255)
                                                                             cAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTc    <  2:852410/91‑1 (MQ=255)
                                                                              aaaGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCt   <  2:31600/91‑1 (MQ=255)
                                                                               aaGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTa  >  1:270226/1‑91 (MQ=255)
                                                                                    |                                                                                     
ACGCCAGCCACGCCGGATCAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGATCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTA  >  NZ_CP009273/987482‑987651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ 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
CAGCCAGCGCACGCCAGCCACGCCGGATCAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGATCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTACTTCCA  >  NZ_CP009273/987472‑987657
                                                                                              |                                                                                           
CAGCCAGCGCACGCCAGCCACGCCGGATCAGGCAGAAAAACAGCCGCTGCATATTACGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGT                                                                                       >  SRR3722118.683430/1‑101 (MQ=60)
                        GGATCAGGCAGAAAAACAGCCGCTGCATGTTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGA                                                               >  SRR3722118.557470/1‑101 (MQ=60)
                            CAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTC                                                           <  SRR3722118.575572/101‑1 (MQ=60)
                                        CAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGT                                               <  SRR3722118.407637/101‑1 (MQ=60)
                                                   tcggagatgtgtataagagacagATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGG                                    <  SRR3722118.32647/78‑1 (MQ=60)
                                                           TTGCCATCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGTAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACC                            <  SRR3722118.556803/101‑1 (MQ=60)
                                                                 TCGAACTGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTC                      >  SRR3722118.459915/1‑101 (MQ=60)
                                                                       TGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGAT                <  SRR3722118.279831/101‑1 (MQ=60)
                                                                       TGTATGATAACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGAT                <  SRR3722118.839794/101‑1 (MQ=60)
                                                                               AACGAAGGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTA        >  SRR3722118.275752/1‑101 (MQ=60)
                                                                                     GGCAAAGTGTTCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTACTTCCA  >  SRR3722118.90843/1‑101 (MQ=60)
                                                                                              |                                                                                           
CAGCCAGCGCACGCCAGCCACGCCGGATCAGGCAGAAAAACAGCCGCTGCATATTCCGTTTGCCATCGAACTGTATGATAACGAAGGCAAAGTGATCCCGTTGCAGAAAGGCGGTCATCCGGTGAATTCCGTGCTGAACGTCACTCAGGCGGAACAGACCTTTGTCTTTGATAATGTCTACTTCCA  >  NZ_CP009273/987472‑987657

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