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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,165,187 (A)5→4 intergenic (‑233/‑129) coaA ← / → BW25113_RS20630 type I pantothenate kinase/tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,165,1830A.100.0% 148.4 / NA 36intergenic (‑229/‑133)coaA/BW25113_RS20630type I pantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (23/13);  total (23/13)

AAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCA  >  NZ_CP009273/4165105‑4165263
                                                                              |                                                                                
aaGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTa                                                                      <  1:196469/90‑1 (MQ=255)
    aaTGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCt                                                                  >  2:65951/1‑90 (MQ=255)
            aGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCg                                                          >  2:730772/1‑90 (MQ=255)
              atatTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGtt                                                        >  1:719778/1‑90 (MQ=255)
               tatTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTc                                                       >  2:509154/1‑90 (MQ=255)
               tatTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTc                                                       >  1:422044/1‑90 (MQ=255)
               tatTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTc                                                       >  1:671982/1‑90 (MQ=255)
                 ttCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCgg                                                     >  2:278453/1‑90 (MQ=255)
                     aTCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt                                                 >  1:285482/1‑90 (MQ=255)
                       cTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGttt                                               >  1:286608/1‑90 (MQ=255)
                       cTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGttt                                               >  1:499425/1‑90 (MQ=255)
                          aTAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTttctt                                            >  1:404877/1‑90 (MQ=255)
                                   aTTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTc                                   >  1:517941/1‑90 (MQ=255)
                                   aTTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTc                                   >  1:893131/1‑90 (MQ=255)
                                   aTTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTc                                   >  1:849181/1‑90 (MQ=255)
                                   aTTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATAGTGTc                                   <  2:600773/90‑1 (MQ=255)
                                    ttATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCg                                  >  2:55252/1‑90 (MQ=255)
                                       tGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCAt                               <  1:65951/90‑1 (MQ=255)
                                       tGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCAt                               <  1:342271/90‑1 (MQ=255)
                                        ggTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATa                              <  1:39400/90‑1 (MQ=255)
                                        ggTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATa                              <  2:404877/90‑1 (MQ=255)
                                           ttAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATaaaa                           >  1:60872/1‑90 (MQ=255)
                                           ttAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATaaaa                           >  2:304781/1‑90 (MQ=255)
                                           ttAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATaaaa                           >  2:460572/1‑90 (MQ=255)
                                                aGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGa                      <  1:144918/90‑1 (MQ=255)
                                                aGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGa                      <  2:539869/90‑1 (MQ=255)
                                                  gAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGAcc                    >  1:830621/1‑90 (MQ=255)
                                                   aGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTc                                   >  1:672542/1‑74 (MQ=255)
                                                   aGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTc                                   <  2:672542/74‑1 (MQ=255)
                                                      gcATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGccc                                                              >  1:688796/1‑44 (MQ=38)
                                                      gcATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGccc                                                              <  2:688796/44‑1 (MQ=38)
                                                      gcATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCt                                                   <  1:178259/55‑1 (MQ=255)
                                                      gcATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCt                                                   >  2:178259/1‑55 (MQ=255)
                                                         ttgttgAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATaaaa             <  2:35503/90‑1 (MQ=255)
                                                                   caATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGc   >  1:504887/1‑90 (MQ=255)
                                                                    aaTGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCa  <  1:304781/90‑1 (MQ=255)
                                                                              |                                                                                
AAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCA  >  NZ_CP009273/4165105‑4165263

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
AGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATG  >  NZ_CP009273/4165088‑4165279
                                                                                                   |                                                                                            
AGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAG                                                                                              <  SRR3722072.423725/100‑1 (MQ=60)
   AAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGT                                                                                          <  SRR3722072.625011/100‑1 (MQ=60)
                 AAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTC                                                                            <  SRR3722072.199861/100‑1 (MQ=60)
                     AATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTT                                                                        >  SRR3722072.732582/1‑100 (MQ=60)
                      ATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTC                                                                       >  SRR3722072.684013/1‑100 (MQ=60)
                      ATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTC                                                                       >  SRR3722072.429707/1‑100 (MQ=60)
                            CAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGT                                                                 >  SRR3722072.290394/1‑100 (MQ=60)
                              GATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTT                                                               >  SRR3722072.291541/1‑100 (MQ=60)
                              GATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTT                                                               >  SRR3722072.508448/1‑100 (MQ=60)
                                 ATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTT                                                            >  SRR3722072.412184/1‑100 (MQ=60)
                                    gctcggagatgtgtataagagacagTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCAT                                                         <  SRR3722072.181342/75‑1 (MQ=60)
                                          GATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTC                                                   >  SRR3722072.527264/1‑100 (MQ=60)
                                          GATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTC                                                   >  SRR3722072.863861/1‑100 (MQ=60)
                                          GATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTC                                                   >  SRR3722072.908506/1‑100 (MQ=60)
                                                  GAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAA                                           >  SRR3722072.61983/1‑100 (MQ=60)
                                                        TGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGAC                                     <  SRR3722072.348301/100‑1 (MQ=60)
                                                        TGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGAC                                     <  SRR3722072.67168/100‑1 (MQ=60)
                                                         GGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACC                                    <  SRR3722072.40122/100‑1 (MQ=60)
                                                         GGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACC                                    >  SRR3722072.845037/1‑100 (MQ=60)
                                                          GTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGctgtct                                   >  SRR3722072.684586/1‑94 (MQ=60)
                                                             TAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCctgtctcttatacacatctccgagcccacgagacgt                                >  SRR3722072.701122/1‑64 (MQ=60)
                                                                 AGGAGCGCATTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAAC                            <  SRR3722072.147500/100‑1 (MQ=60)
                                                                          TTGTTGAGCACAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGC                   >  SRR3722072.513993/1‑100 (MQ=60)
                                                                                    CAATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAG         >  SRR3722072.466377/1‑100 (MQ=60)
                                                                                     AATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGT        <  SRR3722072.310062/100‑1 (MQ=60)
                                                                                      ATGATGTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCtgtctcttatacacatctccgagcccacgagacgtccgattatc       >  SRR3722072.712527/1‑56 (MQ=60)
                                                                                           GTTGAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCctgtctcttatacacatctccgagcccacgagac  >  SRR3722072.322302/1‑66 (MQ=60)
                                                                                                   |                                                                                            
AGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATG  >  NZ_CP009273/4165088‑4165279

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