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 I175 R1 221 18.5 1047116 95.1% 995807 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,595,668 A→G G16G (GGA→GGG yjjZ → DUF1435 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,595,6680AG100.0% 101.2 / NA 28G16G (GGA→GGGyjjZDUF1435 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/10);  total (18/10)

ATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGAT  >  NZ_CP009273/4595582‑4595754
                                                                                      |                                                                                      
atGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTg                                                                                     >  2:41750/1‑90 (MQ=255)
   atatTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTg                                                                                  <  1:38019/90‑1 (MQ=255)
            atcattatcaATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAggg                                                                         >  1:300225/1‑90 (MQ=255)
             tcattatcaATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGc                                                                        >  1:282350/1‑90 (MQ=255)
                  atcaATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATg                                                                   >  2:355816/1‑90 (MQ=255)
                                      cATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGGGTTGCTGCCGGGGTCGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGa                                               >  1:431613/1‑90 (MQ=255)
                                        tGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGGTGTATGCGGATTTATCGTCAGATc                                             >  2:279609/1‑90 (MQ=255)
                                            gCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTg                                         <  1:260993/90‑1 (MQ=255)
                                              aaCGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGc                                       >  2:238002/1‑90 (MQ=255)
                                              aaCGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGc                                       >  2:387336/1‑90 (MQ=255)
                                                 gTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGa                                    >  2:104308/1‑90 (MQ=255)
                                                  tACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGAt                                   <  1:387336/90‑1 (MQ=255)
                                                  tACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGAt                                   <  1:238002/90‑1 (MQ=255)
                                                   aCGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGAtt                                  >  1:135004/1‑90 (MQ=255)
                                                     gCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGt                                <  2:300225/90‑1 (MQ=255)
                                                        gggCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCAt                             >  2:459720/1‑90 (MQ=255)
                                                        gggCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCAt                             >  2:46982/1‑90 (MQ=255)
                                                        gggCAGTGGCTGGGGAGTGTTGCTGCCGGGGGTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGAATGGCAt                             >  2:237843/1‑90 (MQ=255)
                                                           cAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCt                          >  1:377251/1‑90 (MQ=255)
                                                             gTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGa                        >  1:158278/1‑90 (MQ=255)
                                                               ggCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATg                      <  2:142972/90‑1 (MQ=255)
                                                               ggCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATg                      <  1:459720/90‑1 (MQ=255)
                                                                    gggAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGAtt                 <  2:325636/90‑1 (MQ=255)
                                                                     ggAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTa                <  2:331409/90‑1 (MQ=255)
                                                                       aGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGt              >  1:338849/1‑90 (MQ=255)
                                                                          gtTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTAtt           >  1:454340/1‑90 (MQ=255)
                                                                             gctgcCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCgg                                     <  2:203699/61‑1 (MQ=255)
                                                                             gctgcCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCgg                                     >  1:203699/1‑61 (MQ=255)
                                                                                   ggggTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCgat  <  2:509863/90‑1 (MQ=255)
                                                                                      |                                                                                      
ATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGAT  >  NZ_CP009273/4595582‑4595754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ 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
GATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCAC  >  NZ_CP009273/4595584‑4595767
                                                                                    |                                                                                                   
GATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGGGTTGCTGCCGGGGTTGCTGATTGCAGGG                                                                                      >  SRR3722048.304720/1‑100 (MQ=60)
 ATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGC                                                                                     >  SRR3722048.286546/1‑100 (MQ=60)
 ATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGC                                                                                     <  SRR3722048.38547/100‑1 (MQ=60)
                          GAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGGGTTGCTGCCGGGGTCGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGA                                                            >  SRR3722048.438671/1‑100 (MQ=60)
                                       GTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATT                                               >  SRR3722048.136956/1‑100 (MQ=60)
                                          GCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTC                                            <  SRR3722048.264794/100‑1 (MQ=60)
                                               GTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCT                                       >  SRR3722048.383350/1‑100 (MQ=60)
                                                TACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTG                                      <  SRR3722048.241343/100‑1 (MQ=60)
                                                TACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTG                                      <  SRR3722048.393577/100‑1 (MQ=60)
                                                 ACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGA                                     >  SRR3722048.160542/1‑100 (MQ=60)
                                                           GTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGG                           >  SRR3722048.344107/1‑100 (MQ=60)
                                                             GGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTAT                         <  SRR3722048.467299/100‑1 (MQ=60)
                                                              GCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATT                        >  SRR3722048.461815/1‑100 (MQ=60)
                                                                 GGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCctgtctcttatacacatct                     >  SRR3722048.206542/1‑81 (MQ=60)
                                                                                    GTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCtgtctcttatacacatctccgagcccacgagacctctctacatctcgtat  >  SRR3722048.461569/1‑50 (MQ=60)
                                                                                    |                                                                                                   
GATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCAC  >  NZ_CP009273/4595584‑4595767

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