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,284,519 T→C S181P (TCC→CCC)  narP → nitrate/nitrite response regulator protein NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,284,5190TC100.0% 77.7 / NA 24S181P (TCC→CCC) narPnitrate/nitrite response regulator protein NarP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/13);  total (11/13)

GCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCT  >  NZ_CP009273/2284433‑2284604
                                                                                      |                                                                                     
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                    >  2:254661/1‑90 (MQ=255)
            aGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaa                                                                        >  1:331901/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                       <  1:499701/90‑1 (MQ=255)
                            tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                        >  1:187873/1‑90 (MQ=255)
                                 aCGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctgc                                                   <  1:129640/90‑1 (MQ=255)
                                  cGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctgct                                                  <  1:254661/90‑1 (MQ=255)
                                  cGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctgct                                                  <  2:401103/90‑1 (MQ=255)
                                       tGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTa                                             >  1:259901/1‑90 (MQ=255)
                                        ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                            <  2:275236/90‑1 (MQ=255)
                                        ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGAGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                            <  2:61560/90‑1 (MQ=255)
                                         gCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                            <  1:155852/89‑1 (MQ=255)
                                         gCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                            >  2:155852/1‑89 (MQ=255)
                                          cacaGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAAc                                          <  1:95235/90‑1 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                        >  2:116789/1‑72 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                        <  1:116789/72‑1 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                      >  1:43437/1‑90 (MQ=255)
                                                ggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATg                                    >  1:283243/1‑90 (MQ=255)
                                                   tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                               >  2:242321/1‑44 (MQ=255)
                                                   tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                               <  1:242321/44‑1 (MQ=255)
                                                     tCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGc                               <  1:472260/90‑1 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                    >  1:204977/1‑90 (MQ=255)
                                                                 aTTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                   <  1:506972/90‑1 (MQ=255)
                                                                     ccTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcca               >  2:286166/1‑90 (MQ=255)
                                                                                  tatTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCt  <  2:119060/90‑1 (MQ=255)
                                                                                      |                                                                                     
GCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCT  >  NZ_CP009273/2284433‑2284604

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCA  >  NZ_CP009273/2284435‑2284607
                                                                                    |                                                                                        
GTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAA                                                                           >  SRR3722114.335593/1‑100 (MQ=60)
   cggtctcgtgggctcggagatgtgtataagagacaGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAG                                                                        <  SRR3722114.244656/65‑1 (MQ=60)
           GCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATT                                                                <  SRR3722114.505968/100‑1 (MQ=60)
                AGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAA                                                           >  SRR3722114.189713/1‑100 (MQ=60)
                          agagacaGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGT                                                 <  SRR3722114.118028/93‑1 (MQ=60)
                           CTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTA                                                >  SRR3722114.262520/1‑100 (MQ=60)
                               ACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACT                                            <  SRR3722114.130990/100‑1 (MQ=60)
                                CGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTC                                           <  SRR3722114.257200/100‑1 (MQ=60)
                                  AGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAA                                         >  SRR3722114.43857/1‑100 (MQ=60)
                                    tagtacggtctcgtgggctcggagatgtgtataagagacagGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATG                                       <  SRR3722114.257159/59‑1 (MQ=60)
                                    CTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATG                                       >  SRR3722114.286218/1‑100 (MQ=60)
                                      GGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTC                                     <  SRR3722114.157418/100‑1 (MQ=60)
                                        CACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCC                                   <  SRR3722114.96243/100‑1 (MQ=60)
                                                   TCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGG                        <  SRR3722114.478172/100‑1 (MQ=60)
                                                    CAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGC                       >  SRR3722114.206973/1‑100 (MQ=60)
                                                               ATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTC            <  SRR3722114.513331/100‑1 (MQ=60)
                                                                    agagacaggGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTC       <  SRR3722114.311379/91‑1 (MQ=60)
                                                                         TGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCA  <  SRR3722114.478568/100‑1 (MQ=60)
                                                                                    |                                                                                        
GTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCA  >  NZ_CP009273/2284435‑2284607

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