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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,931,448 G→T L265F (TTG→TTT rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,931,4480GT96.5% 104.2 / ‑3.0 29L265F (TTG→TTTrbsKribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (14/14);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTC  >  NZ_CP009273/3931365‑3931534
                                                                                   |                                                                                      
cAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGgaa                                                                                  >  2:450567/1‑90 (MQ=255)
                   gggTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAgg                                                               >  1:313934/1‑90 (MQ=255)
                       gCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGAt                                                           <  2:413420/90‑1 (MQ=255)
                       gCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGAt                                                           <  1:450567/90‑1 (MQ=255)
                                    aCCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAt                                              >  2:105421/1‑90 (MQ=255)
                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                          >  1:204096/1‑90 (MQ=255)
                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                          >  1:209900/1‑90 (MQ=255)
                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                          >  1:362954/1‑90 (MQ=255)
                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                          >  1:162820/1‑90 (MQ=255)
                                        ttGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgccg                                          >  1:165027/1‑88 (MQ=255)
                                                   gaTACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGAtt                               <  2:325758/90‑1 (MQ=255)
                                                     tACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGc                             >  1:336895/1‑90 (MQ=255)
                                                          ttAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTaa                        <  1:49639/90‑1 (MQ=255)
                                                          ttAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTaa                        <  2:209900/90‑1 (MQ=255)
                                                          ttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTaa                        <  2:253627/90‑1 (MQ=255)
                                                                gTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTa                  >  2:564016/1‑90 (MQ=255)
                                                                gTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTa                  >  1:438057/1‑90 (MQ=255)
                                                                  gCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                <  1:507827/90‑1 (MQ=255)
                                                                  gCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                <  2:446124/90‑1 (MQ=255)
                                                                           aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa       <  2:224007/90‑1 (MQ=255)
                                                                           aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa       <  2:394911/90‑1 (MQ=255)
                                                                           aCGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa       <  2:437359/90‑1 (MQ=255)
                                                                            cGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                >  2:451726/1‑80 (MQ=255)
                                                                            cGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                <  1:451726/80‑1 (MQ=255)
                                                                            cGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGca          <  2:577991/86‑1 (MQ=255)
                                                                            cGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGca          >  1:577991/1‑86 (MQ=255)
                                                                            cGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAc      >  2:308611/1‑90 (MQ=255)
                                                                             ggCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAcc     <  1:435815/90‑1 (MQ=255)
                                                                                aTTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTc  <  1:446917/90‑1 (MQ=255)
                                                                                   |                                                                                      
CAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTC  >  NZ_CP009273/3931365‑3931534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ 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
CCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGG  >  NZ_CP009273/3931374‑3931544
                                                                          |                                                                                                
CCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGG                                                                         >  SRR3722111.318061/1‑100 (MQ=60)
              GCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCC                                                           <  SRR3722111.457403/100‑1 (MQ=60)
                     GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTG                                                    >  SRR3722111.164874/1‑100 (MQ=60)
                     GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTG                                                    >  SRR3722111.206616/1‑100 (MQ=60)
                     GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTG                                                    >  SRR3722111.212470/1‑100 (MQ=60)
                     GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTG                                                    >  SRR3722111.368023/1‑100 (MQ=60)
                     GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCCG                                                    >  SRR3722111.167102/1‑100 (MQ=60)
                                  CTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGC                                       >  SRR3722111.341404/1‑100 (MQ=60)
                                             ACCTTTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTA                            >  SRR3722111.444624/1‑100 (MQ=60)
                                                 TTAACGGTGCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGG                        <  SRR3722111.50242/100‑1 (MQ=60)
                                                         GCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAC                <  SRR3722111.458578/100‑1 (MQ=60)
                                                         GCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAC                <  SRR3722111.515687/100‑1 (MQ=60)
                                                         GCGTTAATCACGGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAC                >  SRR3722111.587058/1‑100 (MQ=60)
                                                                    GGCATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCG     <  SRR3722111.442333/100‑1 (MQ=60)
                                                                       ATTTCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGG  <  SRR3722111.453673/100‑1 (MQ=60)
                                                                          |                                                                                                
CCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGG  >  NZ_CP009273/3931374‑3931544

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