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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,796,543 G→C G95G (GGC→GGG thrS ← threonine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,796,5430GC76.7% 61.8 / 13.6 30G95G (GGC→GGGthrSthreonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base C (7/16);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCT  >  NZ_CP009273/1796483‑1796630
                                                            |                                                                                       
ttCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAg                                                            >  1:153551/1‑90 (MQ=255)
          aaCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATc                                                  >  2:837591/1‑90 (MQ=255)
          aaCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATc                                                  >  1:135188/1‑90 (MQ=255)
          aaCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATc                                                  >  1:520323/1‑90 (MQ=255)
           aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcg                                                 >  2:146137/1‑90 (MQ=255)
           aCGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcg                                                 >  1:370069/1‑90 (MQ=255)
            cGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcgc                                                <  2:143983/90‑1 (MQ=255)
            cGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcgc                                                <  1:582257/90‑1 (MQ=255)
             gTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATcgcg                                               <  1:520401/90‑1 (MQ=255)
              tGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGt                                              <  1:61774/90‑1 (MQ=255)
               gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCGAAGTTGTTTAATCGCGTg                                             <  2:252581/90‑1 (MQ=255)
               gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  1:478104/90‑1 (MQ=255)
               gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  2:667048/90‑1 (MQ=255)
               gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  2:627845/90‑1 (MQ=255)
               gCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTg                                             <  2:436858/90‑1 (MQ=255)
                 ggTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGcc                                           >  1:319409/1‑90 (MQ=255)
                         tCAACGTCGTAATAAAAACCGTTGTCAATACCCGGGCCGATTGCCGTTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAgg                                   <  2:690106/90‑1 (MQ=255)
                            aCGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTgc                                >  1:696988/1‑90 (MQ=255)
                            aCGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTg                                 >  1:722067/1‑89 (MQ=255)
                            aCGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTg                                 <  2:722067/89‑1 (MQ=255)
                              gtcgtAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTgcgc                              <  1:408561/90‑1 (MQ=255)
                                  taataaAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAg                          <  2:319409/90‑1 (MQ=255)
                                  taataaAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAg                          <  2:422801/90‑1 (MQ=255)
                                          aCCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACg                  <  2:54855/90‑1 (MQ=255)
                                             gTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAAt               >  1:657803/1‑90 (MQ=255)
                                              ttGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATg              <  1:835268/90‑1 (MQ=255)
                                              ttGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATg              <  1:640316/90‑1 (MQ=255)
                                               tGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGa             <  2:426653/90‑1 (MQ=255)
                                                        ccGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAc    <  2:455673/90‑1 (MQ=255)
                                                        ccGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAc    <  2:716063/90‑1 (MQ=255)
                                                          gggCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACct  >  1:246412/1‑90 (MQ=255)
                                                            |                                                                                       
TTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCT  >  NZ_CP009273/1796483‑1796630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ 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
CTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTC  >  NZ_CP009273/1796473‑1796635
                                                                      |                                                                                            
CTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAG                                                                 >  SRR3722079.154882/1‑100 (MQ=60)
          TTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATC                                                       >  SRR3722079.525123/1‑100 (MQ=60)
          TTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATC                                                       >  SRR3722079.136330/1‑100 (MQ=60)
           TCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCG                                                      >  SRR3722079.373032/1‑100 (MQ=60)
                 GTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCC                                                >  SRR3722079.321979/1‑100 (MQ=60)
                      CGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAAC                                           <  SRR3722079.588019/100‑1 (MQ=60)
                       GTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACA                                          <  SRR3722079.525201/100‑1 (MQ=60)
                        TGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAG                                         <  SRR3722079.62275/100‑1 (MQ=60)
                         GCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGG                                        <  SRR3722079.482294/100‑1 (MQ=60)
                            GTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGC                                     >  SRR3722079.704337/1‑100 (MQ=60)
                            GTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGC                                     >  SRR3722079.729781/1‑100 (MQ=60)
                                        GTCGTAATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGA                         <  SRR3722079.412043/100‑1 (MQ=60)
                                             AATAAAAACCGTTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAAT                    >  SRR3722079.664627/1‑100 (MQ=60)
                                                        TTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAC         <  SRR3722079.646914/100‑1 (MQ=60)
                                                        TTGTCAATAACCGGCCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGAC         <  SRR3722079.844208/100‑1 (MQ=60)
                                                          GTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCT       >  SRR3722079.248483/1‑100 (MQ=60)
                                                               TAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTC  >  SRR3722079.2520/1‑100 (MQ=60)
                                                                      |                                                                                            
CTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAAAACCGTTGTCAATAACCGGGCCGATTGCCATTTTGGTATGCGGCCAAAGTTGTTTAATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTC  >  NZ_CP009273/1796473‑1796635

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