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 I178 R1 167 28.0 1137820 95.3% 1084342 115.1

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% 57.1 / NA 17T155P (ACG→CCG) sulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/9);  total (8/9)

GCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACA  >  NZ_CP009273/1015829‑1016037
                                                                                    |                                                                                                                            
gCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGc                                                                                                                           <  1:405226/88‑1 (MQ=255)
gCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGc                                                                                                                           >  2:405226/1‑88 (MQ=255)
gCGACAAAAAAAGATCCAGTATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGATGCATAATAAACCCCATAGCGTTACCtt                                                                         >  2:196567/1‑138 (MQ=255)
                       aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCa                                                  <  2:293621/138‑1 (MQ=255)
                        aTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAg                                                 >  1:215275/1‑138 (MQ=255)
                              aaTTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAcccc                                                                                      <  1:233675/95‑1 (MQ=255)
                              aaTTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAcccc                                                                                      >  2:233675/1‑95 (MQ=255)
                                   aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTAcc                                                                                                        >  1:32982/1‑72 (MQ=255)
                                   aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTAcc                                                                                                        <  2:32982/72‑1 (MQ=255)
                                   aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATg                                             <  2:79614/131‑1 (MQ=255)
                                   aCTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATg                                             >  1:79614/1‑131 (MQ=255)
                                              cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCataat                                                                                             >  2:359164/1‑72 (MQ=255)
                                              cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCataat                                                                                             <  1:359164/72‑1 (MQ=255)
                                                          aaTTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGccaaccaa              <  2:140045/139‑1 (MQ=255)
                                                               ttAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGa          <  2:215275/138‑1 (MQ=255)
                                                                    ccGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATcacca    >  2:433222/1‑139 (MQ=255)
                                                                      ggAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACcaca  <  2:405991/139‑1 (MQ=255)
                                                                                    |                                                                                                                            
GCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACA  >  NZ_CP009273/1015829‑1016037

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ 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
ACCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACACTGTAATTGCCCGTGC  >  NZ_CP009273/1015772‑1016053
                                                                                                                                             |                                                                                                                                            
ACCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGA                                                                                                                                     <  SRR3722132.486045/151‑1 (MQ=60)
 CCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAG                                                                                                                                    <  SRR3722132.194660/151‑1 (MQ=60)
   GCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGG                                                                                                                                   >  SRR3722132.390890/1‑150 (MQ=60)
      agtacggtctcgtgggctcggagatgtgtataagagacagGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATG                                                                                                                                <  SRR3722132.412589/110‑1 (MQ=60)
                                   gacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAG                                                                                                  <  SRR3722132.365774/94‑1 (MQ=60)
                                          gtctcgtgggctcggagatgtgtataagagaCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACC                                                                                           <  SRR3722132.238053/120‑1 (MQ=60)
                                                                      TTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGC                                                                >  SRR3722132.219290/1‑150 (MQ=60)
                                                                                 ATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtc                                                    >  SRR3722132.33562/1‑94 (MQ=60)
                                                                                 ATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTC                                                    >  SRR3722132.81057/1‑151 (MQ=60)
                                                                                                                                    GTTGTCTCGGGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaa  >  SRR3722132.221998/1‑82 (MQ=60)
                                                                                                                                             |                                                                                                                                            
ACCGCTTCAGACAAGCCTCCGCAAGGAAAATTAGTCACGACTGAAAGCATTGGCTGGGCGACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCATTTGCCGCATCAACAAGTTCAGCATGCTCTTCTTCAGTCAAATCATCTGCCAACCAACCGATCACCACACTGTAATTGCCCGTGC  >  NZ_CP009273/1015772‑1016053

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