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 557,194 A→G Q681R (CAG→CGG)  sfmD → fimbrial biogenesis usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273557,1940AG100.0% 61.5 / NA 18Q681R (CAG→CGG) sfmDfimbrial biogenesis usher protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/5);  total (13/5)

AGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGA  >  NZ_CP009273/557061‑557326
                                                                                                                                     |                                                                                                                                    
aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCt                                                                                                                                   <  2:234453/137‑1 (MQ=255)
    aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTgg                                                                                                                             >  2:351497/1‑139 (MQ=255)
                            cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg                                                                                                     >  2:27631/1‑139 (MQ=255)
                                         cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc                                                                                        >  1:525285/1‑139 (MQ=255)
                                          acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcc                                                                                       <  2:525285/139‑1 (MQ=255)
                                                        ggCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGa                                                                          >  1:27562/1‑138 (MQ=255)
                                                                         cTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag                                                                                         <  2:102659/106‑1 (MQ=255)
                                                                         cTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag                                                                                         >  1:102659/1‑106 (MQ=255)
                                                                                       gggggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAg                                          <  2:208283/139‑1 (MQ=255)
                                                                                          ggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAACCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCg                                       >  2:397620/1‑139 (MQ=255)
                                                                                              ggTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgtg                                    >  1:384945/1‑138 (MQ=255)
                                                                                                         gACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTaa                                                        >  1:307822/1‑107 (MQ=255)
                                                                                                         gACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTaa                                                        <  2:307822/107‑1 (MQ=255)
                                                                                                         gACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaa                         >  1:143217/1‑138 (MQ=255)
                                                                                                          aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa                       >  1:308959/1‑139 (MQ=255)
                                                                                                          aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa                       >  1:190018/1‑139 (MQ=255)
                                                                                                              cgATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa                    >  1:35414/1‑138 (MQ=255)
                                                                                                                               aCTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGa  >  2:239974/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                    
AGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGA  >  NZ_CP009273/557061‑557326

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCG  >  NZ_CP009273/557091‑557332
                                                                                                       |                                                                                                                                          
AATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCC                                                                                             >  SRR3722132.534582/1‑151 (MQ=60)
               GGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGAT                                                                               >  SRR3722132.28042/1‑150 (MQ=60)
                                CTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGctgtctcttatacacatctccg                                                              >  SRR3722132.104574/1‑128 (MQ=60)
                                                     GCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGT                                         >  SRR3722132.391961/1‑150 (MQ=60)
                                                                GGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCtgtctcttatacacatctccg                             >  SRR3722132.313610/1‑130 (MQ=60)
                                                                GGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAA                              >  SRR3722132.145861/1‑150 (MQ=60)
                                                                 GTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAAT                            >  SRR3722132.193581/1‑151 (MQ=60)
                                                                 GTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAAT                            >  SRR3722132.314762/1‑151 (MQ=60)
                                                                     AACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAA                         >  SRR3722132.36030/1‑150 (MQ=60)
                                                                                           ACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCG  >  SRR3722132.187442/1‑151 (MQ=60)
                                                                                                       |                                                                                                                                          
AATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCG  >  NZ_CP009273/557091‑557332

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