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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 579,943 (T)8→7 intergenic (+57/+193) appY → / ← ompT DNA‑binding transcriptional activator AppY/omptin family outer membrane protease OmpT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273579,9360T.92.3% 86.1 / 2.8 26intergenic (+50/+200)appY/ompTDNA‑binding transcriptional activator AppY/omptin family outer membrane protease OmpT
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (11/13);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTT  >  NZ_CP009273/579862‑580018
                                                                          |                                                                                  
tGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGata                                                                    >  2:784859/1‑90 (MQ=255)
  gAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATAtt                                                                  >  1:1193759/1‑90 (MQ=255)
            aacaaTTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCAt                                                        >  2:615355/1‑90 (MQ=255)
                  tGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAAttt                                                  <  1:295729/90‑1 (MQ=255)
                   gACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAAttt                                                  <  1:245096/90‑1 (MQ=255)
                        aTAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGt                                                                                    >  1:14074/1‑51 (MQ=255)
                        aTAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGt                                                                                    <  2:14074/51‑1 (MQ=255)
                        aTAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAAtt                                            >  1:237829/1‑90 (MQ=255)
                                 attaCAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTg                                   <  2:1271109/90‑1 (MQ=255)
                                 attaCAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTg                                   <  1:567436/90‑1 (MQ=255)
                                      aaGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGcc                              <  1:1310641/90‑1 (MQ=255)
                                      aaGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGcc                              <  2:521474/90‑1 (MQ=255)
                                      aaGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGcc                              <  2:694733/90‑1 (MQ=255)
                                      aaGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGcc                              <  1:280754/90‑1 (MQ=255)
                                       aGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAAGTTGTAATTGGGTTGCTGTCGcc                              <  2:1569195/89‑1 (MQ=255)
                                       aGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGCGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGcc                              >  1:1569195/1‑89 (MQ=255)
                                           gTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTAttt                         >  1:255528/1‑90 (MQ=255)
                                            tCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTAtttt                        <  1:1523735/90‑1 (MQ=255)
                                                ccagTTAATTATAATATTATTCCGCGTTTTTTTTATTTTATATTTTTGAACCATTAAATTTTAATTGGGTTGGTGGCCCCTTTTTTatac                    <  1:1549865/87‑1 (MQ=255)
                                                cATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTatac                    <  1:318579/90‑1 (MQ=255)
                                                cATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTatac                    <  1:1260864/90‑1 (MQ=255)
                                                  tGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTatacat                  <  1:1082957/90‑1 (MQ=255)
                                                  tGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTatacat                  <  1:677852/90‑1 (MQ=255)
                                                     tAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACt               >  2:1497021/1‑90 (MQ=255)
                                                         tataATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATaa           >  1:787064/1‑90 (MQ=255)
                                                                  ataCAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGtt  >  2:1570674/1‑90 (MQ=255)
                                                                  ataCAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGtt  >  2:381210/1‑90 (MQ=255)
                                                                  ataCAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGtt  >  2:1183234/1‑90 (MQ=255)
                                                                          |                                                                                  
TGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTT  >  NZ_CP009273/579862‑580018

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ 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
GTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAG  >  NZ_CP009273/579854‑580033
                                                                                         |                                                                                          
GTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATT                                                                                 >  SRR3722078.1212456/1‑100 (MQ=60)
                      CAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATT                                                           >  SRR3722078.240734/1‑100 (MQ=60)
                          TGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGT                                                       <  SRR3722078.299337/100‑1 (MQ=60)
                           GACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGT                                                       <  SRR3722078.248092/100‑1 (MQ=60)
                                     GTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGCGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCT                                            >  SRR3722078.1594739/1‑100 (MQ=60)
                                         ATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTT                                        <  SRR3722078.574085/100‑1 (MQ=60)
                                         ATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTT                                        >  SRR3722078.258656/1‑100 (MQ=60)
                                              cgtcagatgtgtataagagacagATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATAC                                   <  SRR3722078.1331197/77‑1 (MQ=60)
                                              agcgtcagatgtgtataagagacagATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATAC                                   <  SRR3722078.1117006/75‑1 (MQ=60)
                                              AAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATAC                                   <  SRR3722078.1331565/100‑1 (MQ=60)
                                              AAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATAC                                   <  SRR3722078.284180/100‑1 (MQ=60)
                                                    TCTACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTA                             <  SRR3722078.1548571/100‑1 (MQ=60)
                                                       ACATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAA                          >  SRR3722078.797344/1‑100 (MQ=60)
                                                        CATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAAT                         <  SRR3722078.1280849/100‑1 (MQ=60)
                                                        CATGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAAT                         <  SRR3722078.322482/100‑1 (MQ=60)
                                                          TGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTG                       <  SRR3722078.1099342/100‑1 (MQ=60)
                                                          TGTTAATTATAATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTG                       <  SRR3722078.685691/100‑1 (MQ=60)
                                                                     ATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGctgtctcttat            >  SRR3722078.1320930/1‑89 (MQ=60)
                                                                     ATATTATACAGCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTA            >  SRR3722078.843792/1‑100 (MQ=60)
                                                                               GCGTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAG  <  SRR3722078.385863/100‑1 (MQ=60)
                                                                                         |                                                                                          
GTCACAGATGGAATAAACAAAACAATTGACTGATAATGTTTATTACAAGTTGTCTACATGTTAATTATAATATTATACAGCGTTTTTTTTGATGTGATATTCTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAG  >  NZ_CP009273/579854‑580033

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