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 I230 R1 226 18.8 1048726 96.7% 1014118 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,903,061 A→T intergenic (‑38/‑121) mgrB ← / → yobH PhoP/PhoQ regulator MgrB/YobH family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,903,0610AT100.0% 95.9 / NA 25intergenic (‑38/‑121)mgrB/yobHPhoP/PhoQ regulator MgrB/YobH family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/11);  total (14/11)

AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTA  >  NZ_CP009273/1902979‑1903147
                                                                                  |                                                                                      
agcAAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCg                                                                                 <  2:131636/90‑1 (MQ=255)
     acacCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCgttagt                                                                            >  1:11256/1‑90 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       >  2:82863/1‑90 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:60527/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       >  2:453623/1‑90 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:247892/90‑1 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       >  1:497213/1‑90 (MQ=255)
          acGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAc                                                                       <  2:11256/90‑1 (MQ=255)
            gACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAg                                                                     <  2:112225/90‑1 (MQ=255)
                  cAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTa                                                               >  1:167310/1‑90 (MQ=255)
                      aCGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTa                                                           <  1:519238/90‑1 (MQ=255)
                           ccATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCg                                                      >  2:98130/1‑90 (MQ=255)
                           ccATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCg                                                      >  2:66797/1‑90 (MQ=255)
                              tCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGata                                                   >  1:232238/1‑90 (MQ=255)
                                                      cGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAGACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCg                           >  1:438604/1‑90 (MQ=255)
                                                          tCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgcttg                       >  1:484718/1‑90 (MQ=255)
                                                                acctTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAg                 <  1:16782/90‑1 (MQ=255)
                                                                     aTGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCat            <  2:369575/90‑1 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgctt                        <  2:263441/73‑1 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCgctt                        >  1:263441/1‑73 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc       >  2:255311/1‑90 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc       >  1:406071/1‑90 (MQ=255)
                                                                          aTTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGc       >  1:159130/1‑90 (MQ=255)
                                                                              aaCTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCt   <  1:98130/90‑1 (MQ=255)
                                                                               aCTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTa  <  2:167310/90‑1 (MQ=255)
                                                                                  |                                                                                      
AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTA  >  NZ_CP009273/1902979‑1903147

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCG  >  NZ_CP009273/1902974‑1903157
                                                                                       |                                                                                                
AAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGT                                                                                      >  SRR3722114.11365/1‑100 (MQ=60)
     AGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAAC                                                                                 >  SRR3722114.503445/1‑100 (MQ=60)
             CCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTA                                                                         >  SRR3722114.168962/1‑100 (MQ=60)
                         GAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATA                                                             >  SRR3722114.234441/1‑100 (MQ=60)
                           ACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATT                                                           <  SRR3722114.525747/100‑1 (MQ=60)
                                                 CTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAGACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCG                                     >  SRR3722114.444015/1‑100 (MQ=60)
                                                     AATCTCCGTTTCACCTACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTG                                 >  SRR3722114.490791/1‑100 (MQ=60)
                                                                     ACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCctgtctc                 >  SRR3722114.266127/1‑93 (MQ=60)
                                                                     ACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC                 >  SRR3722114.160721/1‑100 (MQ=60)
                                                                     ACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC                 <  SRR3722114.16945/100‑1 (MQ=60)
                                                                     ACCTTATGTCATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGC                 >  SRR3722114.411032/1‑100 (MQ=60)
                                                                              CATTAAACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGAT        >  SRR3722114.324346/1‑100 (MQ=60)
                                                                                   AACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGC   <  SRR3722114.99167/100‑1 (MQ=60)
                                                                                    ACTTGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCG  >  SRR3722114.507639/1‑100 (MQ=60)
                                                                                       |                                                                                                
AAGCAAGCCAACACCACGACAACCAGAACGACCCATCGAAACTTTTTCACTCCAATCTCCGTTTCACCTACCTTATGTCATTAAACTAGCATGCGTTAGTTAAACAGCGCCTAACTATGTCGATATTCGTGCTTTTTCGGAATGAGTCGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCG  >  NZ_CP009273/1902974‑1903157

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