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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

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

GGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATA  >  NZ_CP009273/2797290‑2797457
                                                                                   |                                                                                    
ggTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAttttcca                                                                                >  2:322798/1‑87 (MQ=255)
   gatgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGcc                                                                            >  1:409074/1‑90 (MQ=255)
     tgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCcgc                                                                          >  2:225196/1‑90 (MQ=255)
       atCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCgg                                                                        <  2:48460/90‑1 (MQ=255)
         ccGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGaa                                                                      >  1:423610/1‑90 (MQ=255)
            cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGTTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                   <  1:471617/90‑1 (MQ=255)
              gAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTg                                                                 <  2:189561/90‑1 (MQ=255)
                    gAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAAt                                                           <  1:322798/90‑1 (MQ=255)
                         tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                      <  2:47741/90‑1 (MQ=255)
                         tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                      <  2:426876/90‑1 (MQ=255)
                                           aaTCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGAt                                    >  1:631345/1‑90 (MQ=255)
                                                  tttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCcgc                             <  1:361609/90‑1 (MQ=255)
                                                       tCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGa                        >  2:753257/1‑90 (MQ=255)
                                                          aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                     >  2:124741/1‑90 (MQ=255)
                                                               ggTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc                >  2:538401/1‑90 (MQ=255)
                                                                gTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCa               <  1:722481/90‑1 (MQ=255)
                                                                 tttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCAc              <  2:213173/90‑1 (MQ=255)
                                                                     tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt          <  1:672717/90‑1 (MQ=255)
                                                                     tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt          <  1:441841/90‑1 (MQ=255)
                                                                     tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt          <  1:388763/90‑1 (MQ=255)
                                                                      ggCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTta         <  1:361614/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  2:107315/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  1:52146/90‑1 (MQ=255)
                                                                            cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt   <  1:304091/90‑1 (MQ=255)
                                                                             gATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATa  <  2:423610/90‑1 (MQ=255)
                                                                                   |                                                                                    
GGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATA  >  NZ_CP009273/2797290‑2797457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
TTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA  >  NZ_CP009273/2797283‑2797466
                                                                                              |                                                                                         
TTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCC                                                                                     >  SRR3722073.417813/1‑100 (MQ=60)
      CGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAA                                                                               >  SRR3722073.432676/1‑100 (MQ=60)
                   CTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGTTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATAC                                                                  <  SRR3722073.481657/100‑1 (MQ=60)
                           GAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACC                                                          <  SRR3722073.329654/100‑1 (MQ=60)
                                        GTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGAT                                             >  SRR3722073.644414/1‑100 (MQ=60)
                                                         TTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAG                            <  SRR3722073.369277/100‑1 (MQ=60)
                                                                gatgtgtataagagacaGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGG                     <  SRR3722073.678455/83‑1 (MQ=60)
                                                                       GTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTAC              <  SRR3722073.737079/100‑1 (MQ=60)
                                                                            TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG         <  SRR3722073.397042/100‑1 (MQ=60)
                                                                            TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG         <  SRR3722073.451290/100‑1 (MQ=60)
                                                                            TGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG         <  SRR3722073.686471/100‑1 (MQ=60)
                                                                             GGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGC        <  SRR3722073.369282/100‑1 (MQ=60)
                                                                                   CGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA  <  SRR3722073.310526/100‑1 (MQ=60)
                                                                                   CGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA  <  SRR3722073.53380/100‑1 (MQ=60)
                                                                                              |                                                                                         
TTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA  >  NZ_CP009273/2797283‑2797466

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