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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,064,564 A→G Y185H (TAT→CAT)  yihS ← sulfoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,064,5640AG100.0% 84.7 / NA 26Y185H (TAT→CAT) yihSsulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/15);  total (11/15)

TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTC  >  NZ_CP009273/4064479‑4064642
                                                                                     |                                                                              
tGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACa                                                                            <  1:196401/90‑1 (MQ=255)
       ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                     <  2:216730/90‑1 (MQ=255)
           ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                                 <  1:250057/90‑1 (MQ=255)
           ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                                 >  2:280914/1‑90 (MQ=255)
           ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                                 >  1:530551/1‑90 (MQ=255)
               aCGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTcc                                                             >  1:148434/1‑90 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                      <  2:530551/90‑1 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                      <  2:465502/90‑1 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                      <  1:722553/90‑1 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                      <  2:104592/90‑1 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                      <  2:148434/90‑1 (MQ=255)
                          aaTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCat                                                  >  2:96043/1‑90 (MQ=255)
                           aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                                 >  1:135138/1‑90 (MQ=255)
                            tCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCatat                                                <  2:614165/90‑1 (MQ=255)
                                      cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                      <  2:1763/90‑1 (MQ=255)
                                        cgcgAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgccgc                                    <  1:437464/90‑1 (MQ=255)
                                              tCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGt                              <  1:117266/90‑1 (MQ=255)
                                                gcgcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAc                            >  2:567798/1‑90 (MQ=255)
                                                  gcgATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGGActc                          >  1:149933/1‑90 (MQ=255)
                                                               tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg             >  1:523632/1‑90 (MQ=255)
                                                               tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg             >  2:31650/1‑90 (MQ=255)
                                                                      cGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCtt      <  2:282283/90‑1 (MQ=255)
                                                                        tgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCcgccgc                                 >  1:642239/1‑61 (MQ=255)
                                                                        tgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCcgccgc                                 <  2:642239/61‑1 (MQ=255)
                                                                         gagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGt   >  2:38737/1‑90 (MQ=255)
                                                                          agtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTc  <  2:674137/90‑1 (MQ=255)
                                                                                     |                                                                              
TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTC  >  NZ_CP009273/4064479‑4064642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCACATC  >  NZ_CP009273/4064479‑4064661
                                                                                     |                                                                                                 
TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAG                                                                                     <  SRR3722112.198796/100‑1 (MQ=60)
 GATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGC                                                                                    >  SRR3722112.538414/1‑100 (MQ=60)
     ATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCC                                                                                >  SRR3722112.150273/1‑100 (MQ=60)
           GGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCG                                                                          <  SRR3722112.253013/100‑1 (MQ=60)
                 GTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATA                                                                    >  SRR3722112.136825/1‑100 (MQ=60)
                      GGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                               <  SRR3722112.733832/100‑1 (MQ=60)
                                        CGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTC                                             <  SRR3722112.443481/100‑1 (MQ=60)
                                        CGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGGACTC                                             >  SRR3722112.151787/1‑100 (MQ=60)
                                              TCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGT                                       <  SRR3722112.118727/100‑1 (MQ=60)
                                                     ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                                >  SRR3722112.531353/1‑100 (MQ=60)
                                                              TTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGCACTCTTCctgtctcttatacacatct                       >  SRR3722112.652230/1‑81 (MQ=60)
                                                                                GTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC     >  SRR3722112.272559/1‑100 (MQ=60)
                                                                                   ATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCACATC  <  SRR3722112.382805/100‑1 (MQ=60)
                                                                                     |                                                                                                 
TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCACATC  >  NZ_CP009273/4064479‑4064661

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