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,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% 165.4 / NA 38coding (513/960 nt)nrdFclass 1b ribonucleoside‑diphosphate reductase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/26);  total (12/26)

GGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATA  >  NZ_CP009273/2797290‑2797457
                                                                                   |                                                                                    
ggTGATGATCCGCTGAAAAAGAAAATCGCCAGTGGGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAttttcca                                                                                >  1:283375/1‑87 (MQ=255)
   gatgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGcc                                                                            >  2:194205/1‑90 (MQ=255)
     tgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCcgc                                                                          >  1:87476/1‑90 (MQ=255)
          cGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGaaa                                                                     <  2:114337/90‑1 (MQ=255)
           gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                    >  1:1069207/1‑90 (MQ=255)
           gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                    >  2:473416/1‑90 (MQ=255)
            cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                   <  2:304545/90‑1 (MQ=255)
              gAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTg                                                                 <  2:409210/90‑1 (MQ=255)
                    gAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAAt                                                           <  2:341469/90‑1 (MQ=255)
                         tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                      <  1:751396/90‑1 (MQ=255)
                         tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                      <  2:770298/90‑1 (MQ=255)
                         tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                      <  1:544812/90‑1 (MQ=255)
                         tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                      <  2:676319/90‑1 (MQ=255)
                              aGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGAcc                                                 >  1:998536/1‑90 (MQ=255)
                                           aaTCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGAt                                    >  2:146712/1‑90 (MQ=255)
                                                tttttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                               <  2:1135833/90‑1 (MQ=255)
                                                tttttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc                               >  1:637213/1‑90 (MQ=255)
                                                       tCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAACCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGa                        <  2:512068/90‑1 (MQ=255)
                                                          aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                     >  2:162842/1‑90 (MQ=255)
                                                          aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                     <  2:318397/90‑1 (MQ=255)
                                                               ggTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc                <  2:520320/90‑1 (MQ=255)
                                                               ggTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc                <  1:468842/90‑1 (MQ=255)
                                                               ggTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc                >  1:313333/1‑90 (MQ=255)
                                                               ggTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc                <  2:703281/90‑1 (MQ=255)
                                                                gTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCa               <  1:584616/90‑1 (MQ=255)
                                                                 tttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCAc              <  2:900142/90‑1 (MQ=255)
                                                                  ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGCCCACg             <  2:46757/90‑1 (MQ=255)
                                                                    cTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGt           >  1:207307/1‑90 (MQ=255)
                                                                     tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt          <  2:628525/90‑1 (MQ=255)
                                                                     tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt          <  1:75166/90‑1 (MQ=255)
                                                                     tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt          <  1:1106602/90‑1 (MQ=255)
                                                                         tGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTtacta      >  1:990062/1‑90 (MQ=255)
                                                                          gCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTtactac     >  2:1170481/1‑90 (MQ=255)
                                                                           ccGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACa    <  2:892555/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  2:286698/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  2:1069207/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  2:695403/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  2:72136/90‑1 (MQ=255)
                                                                             gATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATa  <  1:136346/90‑1 (MQ=255)
                                                                                   |                                                                                    
GGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATA  >  NZ_CP009273/2797290‑2797457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
AACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  >  NZ_CP009273/2797279‑2797470
                                                                                                  |                                                                                             
AACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCC                                                                                              <  SRR3722089.704528/100‑1 (MQ=60)
 ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCC                                                                                             >  SRR3722089.968515/1‑100 (MQ=60)
 ACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGGGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCA                                                                                            >  SRR3722089.287001/1‑100 (MQ=60)
      ATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGC                                                                                       >  SRR3722089.88554/1‑100 (MQ=60)
            GTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAG                                                                                 >  SRR3722089.1085259/1‑100 (MQ=60)
                               GAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACC                                                              >  SRR3722089.1013379/1‑100 (MQ=60)
                                    TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                                         <  SRR3722089.551151/100‑1 (MQ=60)
                                    TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC                                                         <  SRR3722089.761378/100‑1 (MQ=60)
                                                 TCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCC                                            >  SRR3722089.645338/1‑100 (MQ=60)
                                                                GTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCC                             >  SRR3722089.317307/1‑100 (MQ=60)
                                                                     ATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGT                        >  SRR3722089.209915/1‑100 (MQ=60)
                                                                          GGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTA                   <  SRR3722089.474445/100‑1 (MQ=60)
                                                                          GGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTA                   >  SRR3722089.1004757/1‑100 (MQ=60)
                                                                           GTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTAC                  <  SRR3722089.591645/100‑1 (MQ=60)
                                                                                TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG             <  SRR3722089.76080/100‑1 (MQ=60)
                                                                                TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG             <  SRR3722089.1123336/100‑1 (MQ=60)
                                                                                      CCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCtgtctcttatacacatctgacgctgccgacgaccatcttagtgtagatc       >  SRR3722089.1136274/1‑51 (MQ=60)
                                                                                        GATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAAT     <  SRR3722089.138006/100‑1 (MQ=60)
                                                                                           GTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTnnnnnnnnnnnnnnnnnnnnn  >  SRR3722089.218833/1‑79 (MQ=60)
                                                                                           GTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  >  SRR3722089.213069/1‑100 (MQ=60)
                                                                                           GTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  >  SRR3722089.783223/1‑100 (MQ=60)
                                                                                                  |                                                                                             
AACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATC  >  NZ_CP009273/2797279‑2797470

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