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 I230 R1 226 18.8 1048726 96.7% 1014118 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,797,377 (T)5→4 coding (517/960 nt) nrdF → class 1b ribonucleoside‑diphosphate reductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,797,3730T.100.0% 95.4 / NA 23coding (513/960 nt)nrdFclass 1b ribonucleoside‑diphosphate reductase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (8/15);  total (8/15)

GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAT  >  NZ_CP009273/2797291‑2797456
                                                                                  |                                                                                   
gTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAttttccag                                                                              >  2:258150/1‑86 (MQ=255)
  gatgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGcc                                                                           >  1:279020/1‑90 (MQ=255)
          gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                   >  1:25155/1‑90 (MQ=255)
          gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                   >  2:347627/1‑90 (MQ=255)
          gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                   >  1:331996/1‑90 (MQ=255)
           cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                  <  1:195285/90‑1 (MQ=255)
           cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                  <  2:183347/90‑1 (MQ=255)
                aaaGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGAcc                                                             <  2:390947/90‑1 (MQ=255)
                        tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                     <  2:348092/90‑1 (MQ=255)
                        tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                     <  2:331996/90‑1 (MQ=255)
                        tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                     <  1:398023/90‑1 (MQ=255)
                        tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                     <  1:84997/90‑1 (MQ=255)
                                  gtTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATc                                           >  2:24777/1‑90 (MQ=255)
                                     tCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGt                                        >  1:166573/1‑90 (MQ=255)
                                      cTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTc                                       <  1:520307/90‑1 (MQ=255)
                                                              ggTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc               >  1:123320/1‑90 (MQ=255)
                                                                tttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCAc             <  1:24777/90‑1 (MQ=255)
                                                                 ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg            <  2:166573/90‑1 (MQ=255)
                                                                 ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg            >  2:361611/1‑90 (MQ=255)
                                                                    tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt         <  1:59980/90‑1 (MQ=255)
                                                                           cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt  <  1:5631/90‑1 (MQ=255)
                                                                           cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt  <  1:371235/90‑1 (MQ=255)
                                                                           cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt  <  1:347627/90‑1 (MQ=255)
                                                                           cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt  <  1:258150/90‑1 (MQ=255)
                                                                                  |                                                                                   
GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAT  >  NZ_CP009273/2797291‑2797456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  >  NZ_CP009273/2797280‑2797470
                                                                                                 |                                                                                             
ANATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCA                                                                                            >  SRR3722114.367961/1‑100 (MQ=60)
   TTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCC                                                                                         >  SRR3722114.281950/1‑100 (MQ=60)
           GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAG                                                                                 >  SRR3722114.25400/1‑100 (MQ=60)
           GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAG                                                                                 >  SRR3722114.335690/1‑100 (MQ=60)
                      CTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATAC                                                                      <  SRR3722114.197187/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                                         <  SRR3722114.85903/100‑1 (MQ=60)
                                   TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                                         <  SRR3722114.402855/100‑1 (MQ=60)
                                      CCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGT                                                      >  SRR3722114.168222/1‑100 (MQ=60)
                                                 CTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCG                                           <  SRR3722114.526832/100‑1 (MQ=60)
                                                               GTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCC                             >  SRR3722114.124611/1‑100 (MQ=60)
                                                                           TTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACA                 <  SRR3722114.25019/100‑1 (MQ=60)
                                                                               TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG             <  SRR3722114.60573/100‑1 (MQ=60)
                                                                                     CCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAA       >  SRR3722114.503472/1‑100 (MQ=60)
                                                                                      CGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA      <  SRR3722114.260745/100‑1 (MQ=60)
                                                                                      CGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA      <  SRR3722114.351607/100‑1 (MQ=60)
                                                                                      CGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA      <  SRR3722114.375673/100‑1 (MQ=60)
                                                                                      CGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA      <  SRR3722114.5688/100‑1 (MQ=60)
                                                                                          GTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  <  SRR3722114.517346/100‑1 (MQ=60)
                                                                                                 |                                                                                             
ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  >  NZ_CP009273/2797280‑2797470

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: