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
A3 F1 I1 R1 781 46.9 2231650 93.5% 2086592 109.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,583,454 G→A intergenic (+36/‑29) yrhA → / → insA pseudogene, interrupted by IS1E/IS1 repressor TnpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,583,4540GA100.0% 107.5 / NA 35intergenic (+36/‑29)yrhA/insApseudogene, interrupted by IS1E/IS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/16);  total (19/16)

TAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATC  >  NC_000913/3583324‑3583564
                                                                                                                                  |                                                                                                              
tAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTatgat                                                                                                        >  1:75707/1‑139 (MQ=255)
             ttgAAATAATGATGATTTCATAAACCCTCATCTACATTAACGGTTAGTGATCGGAGATTATAGCACTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATATAATGTATGATGGTGTTTTTGAgg                                                                                           <  2:541280/137‑1 (MQ=37)
                                       cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGc                                                                               >  2:563649/1‑125 (MQ=255)
                                       cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGc                                                                               <  1:563649/125‑1 (MQ=255)
                                          aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCAATTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                                                              >  2:967426/1‑139 (MQ=255)
                                                    aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAg                                                                                   <  2:725753/108‑1 (MQ=255)
                                                    aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAg                                                                                   >  1:313957/1‑108 (MQ=255)
                                                    aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAg                                                                                   <  2:313957/108‑1 (MQ=255)
                                                    aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAg                                                                                   >  1:725753/1‑108 (MQ=255)
                                                     cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAgg                                                                                           >  1:130642/1‑99 (MQ=255)
                                                     cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAgg                                                                                           <  2:130642/99‑1 (MQ=255)
                                                     cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGt                                                   >  1:1110899/1‑139 (MQ=255)
                                                         tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                             >  1:1024430/1‑109 (MQ=255)
                                                         tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                             <  2:1024430/109‑1 (MQ=255)
                                                                cGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                             >  1:1014321/1‑102 (MQ=255)
                                                                cGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                             <  2:1014321/102‑1 (MQ=255)
                                                                  gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTatga                                                                                                         <  1:580142/72‑1 (MQ=255)
                                                                  gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTatga                                                                                                         >  2:580142/1‑72 (MQ=255)
                                                                  gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                      <  2:1110899/139‑1 (MQ=255)
                                                                  gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                      >  1:361232/1‑139 (MQ=255)
                                                                     aTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                                   <  2:361232/139‑1 (MQ=255)
                                                                         taGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtggt                               >  1:203410/1‑139 (MQ=255)
                                                                             aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                                                                 >  1:175871/1‑101 (MQ=255)
                                                                             aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                                                                 <  2:175871/101‑1 (MQ=255)
                                                                              tttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGATGGTGCGTa                          >  1:710646/1‑139 (MQ=255)
                                                                                 cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATgg                                                                                                      >  1:612055/1‑60 (MQ=255)
                                                                                 cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATgg                                                                                                      <  2:612055/60‑1 (MQ=255)
                                                                                   atatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGc                     <  2:492060/139‑1 (MQ=255)
                                                                                    tatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATc                                                                  <  2:350044/93‑1 (MQ=255)
                                                                                    tatTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATc                                                                  >  1:350044/1‑93 (MQ=255)
                                                                                              tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                             <  2:639403/72‑1 (MQ=255)
                                                                                              tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                             >  1:639403/1‑72 (MQ=255)
                                                                                              tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgc          >  2:537098/1‑139 (MQ=37)
                                                                                              tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgc          >  2:309467/1‑139 (MQ=37)
                                                                                                      aaGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATc  <  1:886148/139‑1 (MQ=11)
                                                                                                                                  |                                                                                                              
TAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATC  >  NC_000913/3583324‑3583564

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

N/A