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
A2 F5 I179 R1 182 40.3 1547118 96.8% 1497610 120.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,015,913 T→G T155P (ACG→CCG)  sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,015,9130TG100.0% 81.6 / NA 24T155P (ACG→CCG) sulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/11);  total (13/11)

AGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACA  >  NZ_CP009273/1015780‑1016037
                                                                                                                                     |                                                                                                                            
aGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGt                                                                                                                         >  1:403837/1‑139 (MQ=255)
      gCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGagag                                                                                                                    >  1:117136/1‑138 (MQ=255)
      gCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGagag                                                                                                                    >  2:570885/1‑138 (MQ=255)
                            aCGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCataat                                                                                             <  2:117136/139‑1 (MQ=255)
                                    aaGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa                                                                                     <  2:403837/139‑1 (MQ=255)
                                          tGGCTGGGCGATAAAAAAAGTTCCAGGCTTAATCCTAAATTTACTTAATGATACACATTAGTGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGTCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGt                                                                               >  2:184152/1‑139 (MQ=255)
                                          tGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGt                                                                               >  1:461373/1‑139 (MQ=255)
                                                 gCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCtt                                                                         >  1:398105/1‑138 (MQ=255)
                                                                        aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACg                                                                                                   <  1:485863/89‑1 (MQ=255)
                                                                        aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACg                                                                                                   >  2:485863/1‑89 (MQ=255)
                                                                         aTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAcccc                                                                                      <  1:208237/101‑1 (MQ=255)
                                                                         aTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAcccc                                                                                      >  2:208237/1‑101 (MQ=255)
                                                                              aaaTTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGct                                           >  1:523737/1‑139 (MQ=255)
                                                                                    aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGc                                                <  1:694352/128‑1 (MQ=255)
                                                                                    aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGc                                                >  2:694352/1‑128 (MQ=255)
                                                                                    aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTcttctt                                     >  2:230508/1‑139 (MQ=255)
                                                                                         aTGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTc                                <  1:575311/139‑1 (MQ=255)
                                                                                                      gagTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCataat                                                                                             <  1:99517/65‑1 (MQ=255)
                                                                                                      gagTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCataat                                                                                             >  2:99517/1‑65 (MQ=255)
                                                                                                      gagTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGcca                   <  2:398105/139‑1 (MQ=255)
                                                                                                           aaTTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGccaaccaa              <  2:507545/139‑1 (MQ=255)
                                                                                                                       ggAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATca       >  1:558548/1‑134 (MQ=255)
                                                                                                                       ggAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATca       <  2:558548/134‑1 (MQ=255)
                                                                                                                       ggAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACcaca  <  2:91735/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                            
AGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACA  >  NZ_CP009273/1015780‑1016037

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAAACCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACACTGTAATTGCCCGTGCGTA  >  NZ_CP009273/1015769‑1016056
                                                                                                                                                |                                                                                                                                               
GAAACCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTG                                                                                                                                           >  SRR3722133.409429/1‑151 (MQ=60)
      GCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGG                                                                                                                                      >  SRR3722133.118702/1‑150 (MQ=60)
                                agtacggtctcgtgggctcggagatgtgtataagagacagGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACC                                                                                                           <  SRR3722133.492422/111‑1 (MQ=60)
                                      agcagaagacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAG                                                                                                     <  SRR3722133.100814/87‑1 (MQ=60)
                                          ACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTT                                                                                                 >  SRR3722133.467640/1‑151 (MQ=60)
                                              gggctcggagatgtgtataagagacAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACC                                                                                              <  SRR3722133.211015/125‑1 (MQ=60)
                                                 GCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTC                                                                                           >  SRR3722133.403636/1‑150 (MQ=60)
                                                                              GGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTC                                                             >  SRR3722133.530779/1‑151 (MQ=60)
                                                                                   GATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTT                                                        <  SRR3722133.703361/151‑1 (MQ=60)
                                                                                                   AATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGC                                        <  SRR3722133.582955/151‑1 (MQ=60)
                                                                                                                       ATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACAC                    >  SRR3722133.565986/1‑151 (MQ=60)
                                                                                                                                          GTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCtgtctcttatacacatctccgagcccacg  >  SRR3722133.338974/1‑121 (MQ=60)
                                                                                                                                                |                                                                                                                                               
GAAACCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACACTGTAATTGCCCGTGCGTA  >  NZ_CP009273/1015769‑1016056

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