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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 987,365 A→G N430D (AAT→GAT)  pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273987,3650AG100.0% 142.5 / NA 42N430D (AAT→GAT) pepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

GGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCG  >  NZ_CP009273/987284‑987445
                                                                                 |                                                                                
gggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAt                                                                          >  2:758748/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  2:1175738/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  2:386639/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:138744/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:643427/90‑1 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:6118/90‑1 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:928633/90‑1 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:443565/1‑90 (MQ=255)
                       tCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                   >  1:1051782/1‑90 (MQ=255)
                            atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  2:492043/90‑1 (MQ=255)
                            atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  1:930602/90‑1 (MQ=255)
                            atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  1:248274/90‑1 (MQ=255)
                              ggTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            <  1:347011/90‑1 (MQ=255)
                                tAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTc                                          <  2:575988/90‑1 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTcc                                                             >  1:665715/1‑70 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTcc                                                             <  2:665715/70‑1 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  2:927208/75‑1 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:927208/1‑75 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         >  2:797156/1‑90 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         <  1:300272/90‑1 (MQ=255)
                                  gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        >  1:1082597/1‑90 (MQ=255)
                                  gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        >  2:1213787/1‑90 (MQ=255)
                                    gCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                      <  1:1241307/90‑1 (MQ=255)
                                        cGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTacac                                  <  1:386639/90‑1 (MQ=255)
                                        cGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTacac                                  >  2:350230/1‑90 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        >  1:1183634/1‑83 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        <  2:1183634/83‑1 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGCACAcc                                 >  2:1013418/1‑90 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAACTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                 >  2:943179/1‑90 (MQ=255)
                                            cTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGAt                              >  2:824413/1‑90 (MQ=255)
                                                 acgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                         <  2:517878/90‑1 (MQ=255)
                                                   gacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGAcc                       <  1:874402/90‑1 (MQ=255)
                                                        tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                  <  1:453799/90‑1 (MQ=255)
                                                          tgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAg                <  2:138744/90‑1 (MQ=255)
                                                             aGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacg             <  1:824413/90‑1 (MQ=255)
                                                              ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga            >  2:13141/1‑90 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            >  2:655478/1‑53 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            <  1:655478/53‑1 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         >  2:1157697/1‑56 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         <  1:1157697/56‑1 (MQ=255)
                                                                        gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg  >  2:168124/1‑90 (MQ=255)
                                                                        gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg  <  2:1114960/90‑1 (MQ=255)
                                                                                 |                                                                                
GGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCG  >  NZ_CP009273/987284‑987445

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCC  >  NZ_CP009273/987275‑987464
                                                                                          |                                                                                                   
TTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATC                                                                                            <  SRR3722089.368906/100‑1 (MQ=60)
  CCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTC                                                                                          <  SRR3722089.984778/100‑1 (MQ=60)
                 GCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGT                                                                           >  SRR3722089.140433/1‑100 (MQ=60)
                 GCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGT                                                                           >  SRR3722089.448888/1‑100 (MQ=60)
                      ATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTA                                                                      >  SRR3722089.1067550/1‑100 (MQ=60)
                           GAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCC                                                                 <  SRR3722089.942196/100‑1 (MQ=60)
                           GAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCC                                                                 <  SRR3722089.651648/100‑1 (MQ=60)
                           GAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCC                                                                 <  SRR3722089.6182/100‑1 (MQ=60)
                                TCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACtgtctctta                                                            >  SRR3722089.674316/1‑91 (MQ=60)
                                TCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCctgtc                                                            >  SRR3722089.940741/1‑95 (MQ=60)
                                 CATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCG                                                           >  SRR3722089.1098906/1‑100 (MQ=60)
                                     ATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTAC                                                       <  SRR3722089.251448/100‑1 (MQ=60)
                                     ATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTAC                                                       <  SRR3722089.944194/100‑1 (MQ=60)
                                       gcagcgtcagatgtgtataagagacagGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAC                                                     <  SRR3722089.663912/73‑1 (MQ=60)
                                       GGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAC                                                     <  SRR3722089.351324/100‑1 (MQ=60)
                                        GTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACC                                                    >  SRR3722089.1201655/1‑100 (MQ=60)
                                          gcgtcagatgtgtataagagacagGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGA                                                  <  SRR3722089.1175313/76‑1 (MQ=60)
                                          AGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGA                                                  <  SRR3722089.304087/100‑1 (MQ=60)
                                             GCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTG                                               <  SRR3722089.1260189/100‑1 (MQ=60)
                                                 CGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGAC                                           <  SRR3722089.391356/100‑1 (MQ=60)
                                                            GACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACG                                <  SRR3722089.886957/100‑1 (MQ=60)
                                                                 TGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTAC                           <  SRR3722089.459236/100‑1 (MQ=60)
                                                                      AGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCC                      <  SRR3722089.836030/100‑1 (MQ=60)
                                                                                         GGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACC   <  SRR3722089.627890/100‑1 (MQ=60)
                                                                                          GATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACGCCC  <  SRR3722089.974507/100‑1 (MQ=60)
                                                                                          GATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCC  <  SRR3722089.857761/100‑1 (MQ=60)
                                                                                          GATGTCGATCTCCCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCC  <  SRR3722089.170221/100‑1 (MQ=60)
                                                                                          |                                                                                                   
TTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCC  >  NZ_CP009273/987275‑987464

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