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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 635,449 T→C S81P (TCC→CCC)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273635,4490TC100.0% 77.7 / NA 23S81P (TCC→CCC) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/13);  total (10/13)

ACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTAT  >  NZ_CP009273/635343‑635580
                                                                                                          |                                                                                                                                   
aCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTa                                                                                                     <  2:501256/139‑1 (MQ=255)
      aGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCg                                                                                               <  2:845684/139‑1 (MQ=255)
                           agaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCAt                                                                          <  2:751393/139‑1 (MQ=255)
                           agaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCAt                                                                          <  1:483276/139‑1 (MQ=255)
                            gaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATc                                                                         >  2:372963/1‑139 (MQ=255)
                            gaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTCCAGGCCCCCCGCTGGGCAACGAGTTCACCTCCCTGGTCCTGGCGTTGCTGTGGCCCGGGGGTCATc                                                                         >  1:143325/1‑139 (MQ=255)
                                 cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGt                                                                                              >  1:813307/1‑113 (MQ=255)
                                 cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGt                                                                                              <  2:813307/113‑1 (MQ=255)
                                    cTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGa                                                                   >  2:718099/1‑137 (MQ=255)
                                      tGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAaag                                                                >  2:672301/1‑138 (MQ=255)
                                              cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGt                                                       <  2:751390/139‑1 (MQ=255)
                                                   gCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTctg                                                  <  2:794904/139‑1 (MQ=255)
                                                      gTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTctgctg                                               >  1:21412/1‑139 (MQ=255)
                                                                          cAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGa                           <  2:43115/139‑1 (MQ=255)
                                                                            ggTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGc                                                                                                                               >  1:661711/1‑37 (MQ=255)
                                                                            ggTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGc                                                                                                                               <  2:661711/37‑1 (MQ=255)
                                                                            ggTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACg                         <  1:836729/139‑1 (MQ=255)
                                                                                   aCCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGAttt                  >  1:836494/1‑139 (MQ=255)
                                                                                   aCCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGAttt                  >  1:656998/1‑139 (MQ=255)
                                                                                            cACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGa         <  1:186112/139‑1 (MQ=255)
                                                                                                   tGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTa   <  2:497152/138‑1 (MQ=255)
                                                                                                   tGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTAt  >  1:426964/1‑139 (MQ=255)
                                                                                                    gCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTAt  <  2:135369/138‑1 (MQ=255)
                                                                                                          |                                                                                                                                   
ACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTAT  >  NZ_CP009273/635343‑635580

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 35 ≤ 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
AGAAATCAAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCAC  >  NZ_CP009273/635310‑635598
                                                                                                                                           |                                                                                                                                                     
AGAAATCAAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGC                                                                                                                                            <  SRR3722120.579642/151‑1 (MQ=60)
                                                  CCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTCCAGGCCCCCCGCTGGGCAACGAGTTCACCTCCCTGGTCCTGGCGTTGCTGTGGCCCGGGGGTCATCC                                                                                          >  SRR3722120.145258/1‑151 (MQ=60)
                                                       AAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCctgtctcttatacac                                                                                      >  SRR3722120.824259/1‑135 (MQ=60)
                                                           AAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGA                                                                                 <  SRR3722120.490234/151‑1 (MQ=60)
                                                                            GTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGG                                                                >  SRR3722120.21653/1‑151 (MQ=60)
                                                                                                  TCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaacccacgagacact                                          >  SRR3722120.670785/1‑59 (MQ=60)
                                                                                                         CCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTT                                   >  SRR3722120.666027/1‑151 (MQ=60)
                                                                                                         CCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTT                                   >  SRR3722120.847754/1‑151 (MQ=60)
                                                                                                            AGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAA                                <  SRR3722120.847991/151‑1 (MQ=60)
                                                                                                                         GGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATT                   >  SRR3722120.433193/1‑151 (MQ=60)
                                                                                                                            CCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACT                <  SRR3722120.188702/151‑1 (MQ=60)
                                                                                                                               CGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGC             >  SRR3722120.234375/1‑151 (MQ=60)
                                                                                                                                         GCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCC    >  SRR3722120.345871/1‑150 (MQ=60)
                                                                                                                                          CCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaacacatacacgaaacaaagaacc  >  SRR3722120.260831/1‑52 (MQ=60)
                                                                                                                                           |                                                                                                                                                     
AGAAATCAAGGAACTGTTGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCAC  >  NZ_CP009273/635310‑635598

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