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 F21 I1 R3 778 586.6 19747158 98.0% 19352214 137.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 2,997,620 T→C *144* (TAA→TAG CXP41_RS24400 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252682,997,6200TC100.0% 78.1 / NA 24*144* (TAA→TAGCXP41_RS24400hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/8);  total (16/8)

TGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTA  >  NZ_CP025268/2997488‑2997755
                                                                                                                                    |                                                                                                                                       
tgaAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAATTc                                                                                                                                 >  2:2007087/1‑140 (MQ=255)
tgaAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTc                                                                                                                                 >  2:2007151/1‑140 (MQ=255)
                                                                                                    tataTTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                              <  2:6141972/140‑1 (MQ=255)
                                                                                                      tatTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                      >  2:7074146/1‑98 (MQ=255)
                                                                                                      tatTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                      <  1:7074146/98‑1 (MQ=255)
                                                                                                      tatTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGt                                  >  1:2693017/1‑134 (MQ=255)
                                                                                                      tatTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGt                                  <  2:2693017/134‑1 (MQ=255)
                                                                                                        ttCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                              >  1:410395/1‑136 (MQ=255)
                                                                                                        ttCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                              <  2:410395/136‑1 (MQ=255)
                                                                                                            aaTCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAg                      >  1:7296923/1‑140 (MQ=255)
                                                                                                             aTCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGc                                                   <  1:5455883/110‑1 (MQ=255)
                                                                                                             aTCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGc                                                   >  2:5455883/1‑110 (MQ=255)
                                                                                                             aTCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgg                     >  1:9233871/1‑140 (MQ=255)
                                                                                                             aTCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgg                     >  1:5288534/1‑140 (MQ=255)
                                                                                                                aTTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCg                                            >  2:441668/1‑114 (MQ=255)
                                                                                                                aTTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCg                                            <  1:441668/114‑1 (MQ=255)
                                                                                                                aTTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGga                  >  1:5727552/1‑140 (MQ=255)
                                                                                                                  tAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                              <  1:5990058/126‑1 (MQ=255)
                                                                                                                  tAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                              >  2:5990058/1‑126 (MQ=255)
                                                                                                                         gAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATg         >  2:8787304/1‑140 (MQ=255)
                                                                                                                         gAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATg         >  2:9330578/1‑140 (MQ=255)
                                                                                                                          aGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgggg                   >  1:6968858/1‑129 (MQ=255)
                                                                                                                          aGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgggg                   <  2:6968858/129‑1 (MQ=255)
                                                                                                                                ggAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGAttatta  >  1:8350331/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTA  >  NZ_CP025268/2997488‑2997755

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCA  >  NZ_CP025268/2997470‑2997770
                                                                                                                                                      |                                                                                                                                                      
ATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCT                                                                                                                                                       >  SRR14920172.1098949/1‑151 (MQ=60)
                                                                              tcgtgggctcggagatgtgtataagagacagCAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAA                                                                          <  SRR14920172.7130271/120‑1 (MQ=60)
                                                                                                 agatgtgtataagagacAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGA                                                       <  SRR14920172.5498633/134‑1 (MQ=60)
                                                                                                        gtgtataagagacagATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTT                                                <  SRR14920172.445168/136‑1 (MQ=60)
                                                                                                             CAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGctgtctcttatacacatctccgagcccacgagacgcgtagtaatctcgtatgccgtcttctgctttaaaatggtggggggggggggggggggggggggggggg                                           >  SRR14920172.7425700/1‑48 (MQ=60)
                                                                                                             CAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCT                                           >  SRR14920172.2714160/1‑151 (MQ=60)
                                                                                                               ATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTG                                         >  SRR14920172.413661/1‑151 (MQ=60)
                                                                                                                   GATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAG                                     >  SRR14920172.7354814/1‑151 (MQ=60)
                                                                                                                    ATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGG                                    >  SRR14920172.9307510/1‑151 (MQ=60)
                                                                                                                    ATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGG                                    >  SRR14920172.5329839/1‑151 (MQ=60)
                                                                                                                      cagATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGG                                  <  SRR14920172.6037118/148‑1 (MQ=60)
                                                                                                                       ATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGA                                 >  SRR14920172.5772533/1‑151 (MQ=60)
                                                                                                                                 CATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGA                       >  SRR14920172.7024183/1‑151 (MQ=60)
                                                                                                                                       ACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTA                 >  SRR14920172.8416967/1‑151 (MQ=60)
                                                                                                                                           GAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTA             >  SRR14920172.8273096/1‑151 (MQ=60)
                                                                                                                                            AGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAA            >  SRR14920172.1384844/1‑151 (MQ=60)
                                                                                                                                                     GCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTC   >  SRR14920172.4285624/1‑151 (MQ=60)
                                                                                                                                                     GCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTC   >  SRR14920172.4109665/1‑151 (MQ=60)
                                                                                                                                                     GCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTC   >  SRR14920172.372449/1‑151 (MQ=60)
                                                                                                                                                     GCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTC   >  SRR14920172.8868596/1‑151 (MQ=60)
                                                                                                                                                     GCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTC   >  SRR14920172.2208500/1‑151 (MQ=60)
                                                                                                                                                      CTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCtgtctcttatacacatctccgagcccacgagacgcgtagtaatcgcgtatgccgtcttctgcttgaaaaggggggggggggggggggggggggg  >  SRR14920172.4198680/1‑57 (MQ=60)
                                                                                                                                                      CTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCctgtctcttatacacatctccgagcccacgagacgcgtagta  >  SRR14920172.541370/1‑109 (MQ=60)
                                                                                                                                                      CTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTgtctcttatacacatctccgagcccacgagacgcgtagtaa  >  SRR14920172.6565273/1‑110 (MQ=60)
                                                                                                                                                      CTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCA  >  SRR14920172.4287755/1‑151 (MQ=60)
                                                                                                                                                      CTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGATCA  >  SRR14920172.3614456/1‑151 (MQ=60)
                                                                                                                                                      |                                                                                                                                                      
ATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATTAATTAAATATTGCTCA  >  NZ_CP025268/2997470‑2997770

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