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 I177 R1 141 27.5 996908 97.7% 973979 124.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 724,034 A→C V52G (GTC→GGC)  kdpA ← potassium‑transporting ATPase subunit KdpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273724,0340AC100.0% 59.1 / NA 17V52G (GTC→GGC) kdpApotassium‑transporting ATPase subunit KdpA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/8);  total (9/8)

CCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGACGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTGGCGATC  >  NZ_CP009273/723901‑724165
                                                                                                                                     |                                                                                                                                   
cccTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAgccgcc                                                                                                                                 >  1:394758/1‑138 (MQ=255)
                           cAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACg                                                                                                     >  1:169469/1‑139 (MQ=255)
                               tAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCg                                                                                                 <  2:209908/139‑1 (MQ=255)
                                    cTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTa                                                                                            >  1:458280/1‑139 (MQ=255)
                                       aCCGAGCAACATAAAAAACAGCACCACCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCg                                                                                         >  1:330927/1‑139 (MQ=255)
                                        ccGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTa                                                                                            >  1:247170/1‑135 (MQ=255)
                                        ccGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTa                                                                                            <  2:247170/135‑1 (MQ=255)
                                                 aTAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGaa                                                                               <  2:82069/139‑1 (MQ=255)
                                                          aGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTaa                                                                      >  2:219817/1‑139 (MQ=255)
                                                          aGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTaa                                                                      >  1:496034/1‑139 (MQ=255)
                                                              ccGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAg                                                                  <  2:496034/139‑1 (MQ=255)
                                                                              gTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTg                                                  >  1:470714/1‑139 (MQ=255)
                                                                                    gCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTa                                              >  2:135625/1‑137 (MQ=255)
                                                                                       gagaATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGa                                         <  2:44286/139‑1 (MQ=255)
                                                                                       gagaATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGa                                         <  2:458280/139‑1 (MQ=255)
                                                                                                            cTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAacac                     <  2:169469/138‑1 (MQ=255)
                                                                                                                        ctcACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTg        >  2:142573/1‑139 (MQ=255)
                                                                                                                              gTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTGGCGATc  <  2:327866/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGACGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTGGCGATC  >  NZ_CP009273/723901‑724165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
GCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGACGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTGGCGATCAGTAAGAACCCT  >  NZ_CP009273/723884‑724177
                                                                                                                                                      |                                                                                                                                               
agcagaagacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGC                                                                                                                                                 <  SRR3722131.469844/87‑1 (MQ=60)
      TCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCA                                                                                                                                            >  SRR3722131.399085/1‑150 (MQ=60)
                                 GGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGC                                                                                                                >  SRR3722131.171420/1‑151 (MQ=60)
                                          GGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTAC                                                                                                       >  SRR3722131.463202/1‑151 (MQ=60)
                                             AGATAGTGCTGACCGAGCAACATAAAAAACAGCACCACCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGG                                                                                                    >  SRR3722131.334646/1‑151 (MQ=60)
                                              GATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGG                                                                                                   >  SRR3722131.250083/1‑151 (MQ=60)
                                                                ACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAAT                                                                                 >  SRR3722131.501322/1‑151 (MQ=60)
                                                                                    AGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGC                                                             >  SRR3722131.475777/1‑151 (MQ=60)
                                                                                                                                               GTCAGAGCCGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTGGCGATCAGTAAGAACCCT  >  SRR3722131.162350/1‑151 (MQ=60)
                                                                                                                                                      |                                                                                                                                               
GCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTCAGAGACGCCAAGTGCGCGAAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAAGAGGAATGTCATTAATCAGCCGCGCCAGCCCGCTGCCTAAAGGACGCGCCAGCACCATTAACACCAGTAAAAACGTGGCGATCAGTAAGAACCCT  >  NZ_CP009273/723884‑724177

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